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Response Functions for Types of Vision According 
to the Miller Theory 


By Deane B. Judd 





According to the Miiller theory of vision there are three stages in the visual process, 
an initial photochemical stage, an intermediate chemical stage relating to the chromatic 
aspect, and a final stage of excitations of the optic-nerve fibers. By taking advantage of 
recent precise information regarding the metamers characteristic of normal, protanopic, and 
deuteranopic vision there have been derived the spectral variations of the responses for each 
stage as functions of wavelength. These response functions account precisely for the same 


normal metamers as the ICI standard observer, and closely for the same confusions by 











color-blind observers as the simpler Kénig theory. 


Furthermore these functions describe 


chromatic thresholds of the normal eye (Abney, Priest-Brickwedde) as a gradual approach 


to tritanopic vision as field size and luminance are decreased. 


I. Introduction 


In the ninteenth century, two rival theories of 
vision monopolized most of the interest of in- 
vestigators. One of these is the Young-Helmholtz 
three-components formulation; the other is the 
Hering opponent-colors theory. The discovery of 
the facts of red-green blindness dealt fatal blows 
to the then current forms of both of these simple 
theories, though proponents of the respective 
theories continued to pump a semblance of life 
into them with wordy battles. The opponent- 
colors theory in its original simple form can be 
made to yield but a single form of red-green 
blindness, that known as deuteranopia. It must 
overlook the established fact of a second type of 
red-green blindness, protanopia, in which the 
luminosity function is deficient in the long-wave 
portion of the spectrum and in which the chrom- 
aticity confusions are consistently different from 
those of deuteranopia. The three-components 
theory explains the confusions made by both 
types perfectly but in its original simple form has 
to predict that deuteranopic vision consists of 
mixtures of red and violet and protanopic vision 
consists of mixtures of green and violet. When 
cases of unilateral red-green blindness showed 
consistently that the perceptions of red-green- 





blind observers have the hues blue and yellow and 
no others, the original simple three-components 
formulation became obsolete. Some advocates 
of this simple theory took refuge in a suggestion 
by Fick [1]! that red-green confusion is the result, 
not of the nonfunctioning of either the red or the 
green receptor system, but rather of the two 
receptor systems having identical photosensitive 
substances, either that for red (deuteranopia) or 
that for green (protanopia). By this suggestion, 
the responses from the red cones combine with 
those from the green, regardless of the photo- 
sensitive substance in either, to give yellow. This 
combination can take place in the postretinal 
portion of the nervous system, as emphasized by 
Hecht [2], for binocular fusion of colors, and it 
is permissible to assume that it always takes place 
there even in binary stimulation of one eye alone. 
From this view it is only a step to the theoretical 
position originally proposed by Donders [3] and 
later espoused by Kénig [4], von Kries [5], and 
Adams [6] that the three-components formula- 
tion holds for processes in one stage of the visual 
mechanism (perhaps the photosensitive-substance 
stage), while the opponent-colors theory holds for 


1 Figures in brackets indicate the literature references at the end of this 
paper 

















processes in a later stage (perhaps the optic nerve). 
This view may be called the stage or “zone” 
theory of vision. Furthermore, a very able ad- 
vocate of the opponent-colors theory, G. E. 
Miiller, adopted a theoretical view [7] that, al- 
though divergent in detail and elaborated to 
include an additional stage, was essentially in 
agreement with the stage theories favored by 
Donders, Kénig, von Kries, and Adams. 
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Figure 1. Absolute chromatic thresholds for homogeneous 


light as a function of waveiength. 


The dotted curve is based on the Miller theory; see section VII. @, Ob- 
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It was shown by Abney [8] in 1910, by Priest 
and Brickwedde [9] in 1926, by Guild [10] in 1928, 
by Holmes [11] in 1941, by MacAdam [12] in 1942, 
and probably by others, that the nearly achro- 
matic color of noon sunlight is more confusible with 
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Ficure.2. Least fraction of spectrum light detectable as a 
e ’ ¢ 
chromatic difference in a mixture with sunlight, 2° observing 
field. 


i The solid line is an empirical representation of these data based on the 
uniform-chromaticity-seale triangle (Judd); the dotted curve is based on the 
Miller theory; see section VIT. 
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the greenish yellow color of spectrum light at 570 
mu by a normal observer than with any other 
nearby portion of the spectrum, definitely more 
than with the vellow portion (575 to 585 my). 
Figure | shows the brightness in foot lamberts 
found by Abney to be required to produce the 
perception of a chromatic color noticeably dis- 
tinct from the achromatic color of light from the 
carbon are. Figure 2 shows the Priest-Brick- 
wedde determination of minimum perceptible 
colorimetric purity. In both of these figures the 
maximum near 570 my is outstanding. Similar 
results were found by Guild, Holmes, and Mac- 
Adam. This outstanding maximum might sug- 
gest that pigmentation of the eye media of the 
normal eye absorbs a large fraction of the short- 
wave (violet) portion of the spectrum, or for small 
fields it could mean that because of the chromatic 
aberration of the eye, the short-wave portion of 
the sunlight spectrum is out of focus and largely 
lost. But the most likely explanation is that the 
normal eye, at least in the fovea, has some of the 
characteristics of a tritanopic eye; a tritanope 
has a neutral point in the spectrum near 570 mg 
where the normal observer has this quasi-neutral 
point. Furthermore, for very small fields sub- 
tending 20’ or less, it has been shown by Kénig 
[4], Willmer [13], Hartridge [15], and Wright [14] 
that the fovea is tritanopiec. 

In an attempt to describe the chromaticity 
sensibility of the normal observer in terms of an 
approach to tritanopia, Judd [16] derived a trans- 
formation of the OSA “excitations” corresponding 
to an 80-percent dilution of the violet excitation 
with red and green. Figure 3 shows the resulting 
excitation curves and Maxwell triangle, and Figure 
4 shows how this formulation corresponds with 
Priest's data on minimum perceptible colorimetric 
purity. This formulation corresponds to a the- 
oretical suggestion similar to Fick’s proposal to 
account for red-green confusion; it suggests that 
in the fovea the red and green substance from the 
red and green cones has leaked into the violet 
cones to a serious degree (SO0% leakage). A sim- 
ilar degree of success was demonstrated by Hecht 
[17] in another development of the Young-Helm- 
holtz theory. Both explantions suffer, however, 
from a failure to permit an account of dichromatic 
vision, as do various coordinate systems empiri- 
cally derived to represent in a simple way the 
facts of chromaticity sensibility [18, 19, 20). 
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Ficure 3. The response functions and Maxwell triangle 


corresponding to a three-components explanation (Judd) 


of the data shown in figure 2. 


The suggested explanation is that segregation of photosensitive substances 
specific for long-wave energy is poor, violet cones having nearly as much 
(80°) as the red and green cones. Chromaticity coordinates of the spectrum 
referred to a set of stimulus primaries that have been found useful in deriving 
the “minimum purity perceptible.” (See BS J. Research 4, 515 (1930) R P163 
also J. Opt. Soe. Am. and Rev. Sci. Instr. 16, 115 (1928)). A, The three 
“distribution” curves which give the mixture diagram shown. 


There are, however, two accounts of chroma- 
ticity sensibility that do seem also to permit good 
explanations of dichromatice vision, that by Adams 
21] and the recent excellent treatment of chro- 
maticity sensibility by Stiles [22]. An account of 
protanopia and tritanopia by the Adams theory 
has not yet been worked out in detail. 

An outstanding defect of the three-component 
accounts of chromaticity sensibility is that there 
is no satisfactory explanation of the primary char- 
acter of the spectrum in the neighborhood of 475 
mu. Most normal observers (though not all) see 
this portion of the spectrum as blue, and they see 
the short-wave extreme as binary in character, a 
mixture of red and blue. In commenting on this 
difficulty, it was remarked by Judd [23] in 1932, 
“The most satisfactory solution vet offered is 
Miiller’s theory which ascribes primacy to both 
blue and violet, the latter in the retinal processes, 
and the former in the optic nerve.” 
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10.-$——_—_—_ ——-—— —-—- 
a | 
aa 08 \ | 
iS 4 
7 oO 
Be 
® 3 06 : . 
2x 
r > 05 i _ i 
25.04 ! 
so y | 
2 § 03 ‘ 
H | 
5 02 \ 
. ' 
Ol oe ° ~ . 4 
00 Andante @ 2 Se 
400 450 500 550 600 650 700 


WAVELENGTH IN MILLIMICRONS 
Ficure 4. Least fraction of spectrum light detectable as a 
chromatic difference in a mixture with sunlight, 2° observ- 
ing field. 


The dotted line represents a part of the data shown in figure 2; the solid 
line is based on a three-components explanation (Judd, see fig. 3) of the tend- 
ency of the normal eye under these conditions to make tritanopie confusions. 
© 8KS (c—1/3); K=0.0024, 0.8 blue deficient; +, (dp/dE),-%0. Irwin G Priest, 
observer. 


As a prerequisite to a quantitative explanation 
in terms of the Miiller theory for the confusibility 
of sunlight with the spectrum at 570 my and for 
chromaticity sensibility generally, there must be 
derived the colorimetric coordinate systems corre- 
sponding to the two additional stages, the retinal 
and the optic nerve stages, of the Miller theory. 
They have so far been described only qualitatively, 
or at least semiquantitatively. Since the theory 
has been adjusted to correspond qualitatively with 
the facts of colorblindness, and since those facts 
have recently become known quantitatively (chief- 
ly through the work of Pitt [33}), it is now possible 
to evaluate these coordinate systems, and so lay 
the ground work for a possible explanation of 
chromaticity sensibility based on the Miler 
theory. 


II. Formulation for Normal Vision 


According to the Miller theory? light stimuli 
can elicit: three different primary sensitizing proc- 
esses (P-processes) in the cone mechanism, whose 
strengths are determined according to wave- 
length of the incident radiant energy according to 
functions similar to those defining the three com- 
ponents of the Young-Helmholtz theory. The 
P,-process is aroused by the spectral region 475 

2? Acknowledgment is made to Michael J. Zigler, Department of 


Psychology, Wellesley College, who kindly supplied a very helpful transla- 
tion into English of these parts of Miiller’s discussion of color-blindness [7]. 
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my up to the long-wave visible extreme. The 
P,-process is aroused by the spectral region be- 
tween the long-wave end stretch (770 my to a 
wavelength greater than 655 mz) and the short- 
wave end stretch (380 mug to a wavelength less 
than 450 my). The P;-process is aroused by the 
spectral region between 540 mp and the short- 
wave visible extreme. From this description, the 
distribution curves of the P-processes are seen to 
resemble closely the OSA excitation curves [24]. 
The best modern evaluation of these distribution 
curves based upon the ICI standard observer [25] 
is to be obtained [26] by the following transfor- 
mation: 


P,=3.1956.X +-2.4478Y —0.6434Z 
P= —2.5455.X +-7.0492Y+0.4963Z>), (1) 
P;=0.0000X +-0.0000 Y+5.0000Z 


which are graphed in figure 5. The reverse trans- 
formation [27] is given by 
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WAVELENGTH IN MILLIMICRONS 


Fiat RE 5. Rates of decom position of the photose nsilive 
substances of the three-components theory as functions of 
wavele ngth. 


These functions are suited to the first stage of the Miller theory. 
N=0.24513P,;—0.08512P,+-0.03999P, 
Y=0.08852P,+0.11112P,+-0.00036P,>. (la) 


Z=0.00000P,, +-0.00000P, +-0.20000P, 


Except for the small secondary maximum of the 
P, curve in the neighborhood of 430 mg, these 
functions conform essentially to the description 
by Miller of the primary sensitizing processes. 
These processes contribute immediately to excita- 
tion of a black-white “substance” of the optic- 
nerve fibers in the sense of producing white. They 
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also act upon certain assumed chromatic sensory 
substances in the cones of the normal eye as 
follows: the P,-process acts to produce a major 
yellowish red (yR) process in a yellowish red- 
bluish green (yR-bG) substance and a minor 
greenish yellow (gY) process in a greenish yellow- 
reddish blue (gY-rB) substance. The P.-process 
acts to produce a major bluish green (6G) process 
in the yR-bG@ substance and a minor greenish yel- 
low (gY) process in the gY-rB substance. The 
P;-process acts to produce a reddish blue (rB) 
process in the gY-rB substance. It is further as- 
sumed that the processes within each of these 
chromatic sensory substances are antagonistic so 
that a yR-process cancels completely a bG-process 
of equal strength; and a gY-process may cancel 
completely an rB-process. The wavelengths at 
which stimulation by homogeneous radiant energy 
would produce these cancellations are between 560 
and 570 mu for yR to cancel 6G, and near 495 my 
for gY to cancel rB. 

From this description it would seem that the 
transformation from the amounts of the primary 
sensitizing processes to the amounts of the chro- 
matic sensory processes might take on the simple 
form: 


yR bG=a,P,—aP, 
aY =—rB=b,P,+6.P.—},P,) 


where a, a, 6;, bo, and 6; are constants greater 
than zero, with a, greater than 6,, and a, greater 
than bdo. 

In addition to the primary sensitizing processes 
(P;, Ps, Ps) and the chromatic sensory processes 
(yR-bG, gY-rB), there are six different excita- 
tions of the optic nerve (w, s, r, g, y, 6), which 
correlate with the introspectively pure white, 
black, red, green, yellow, and blue sensations, 
respectively. The chromatic excitations are 
assumed to arise from the chromatic sensory 
processes alone, the yf-process arousing a major 
r-excitation and a minor y-excitation as indicated 
by the notation yR. Similar major and minor 
excitations are aroused by the bG, gY, and rB 
The white-excitation of the optic 
nerve comes from the immediate effect of the 
primary sensitizing processes (P;, P:, P3) to which 
secondary contributions from the yR- and gY- 
processes are added. The black-excitation comes 
chiefly by induction [7a, p. 85], from a white 


processes. 
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surrounding field or from a white preexposure 
field, but a secondary contribution comes from 
the bG and rB processes. Like the chromatic 
sensory processes, the four chromatic excitations 
of the optic nerve make antagonistic pairs, an 
amount of r-excitation cancelling a like amount 
of g-excitation, and the same cancellations for 
y- and 6-excitation. For stimulation by homo- 
geneous radiant energy, the one cancellation 
occurs at the wavelength for arousing unitary 
yellow, which is probably between 575 and 582 my 
under usual observing conditions for an average 
normal observer. A second (r,g)-cancellation 
occurs at the wavelength for arousing unitary 
blue, and a (y,b)-cancellation occurs at the 
wavelength for unitary green. 

For self-luminous areas with a neutral surround- 
ing field the s-excitation acts as negative w-excita- 
tion; however, they do not cancel, but combine to 
give gray. From this description it would appear 
that the excitations of the optic nerve could be 
found for the normal observer by the following 
transformations: 


r=—g=eyR+erB=—e,bG—egY ) 
y= —b=d,gY+dyR=—d,rB—dbG 


w= eP,+e.P2+e3P3+eyR+egY 





s=ebG+-e,rB 


where the luminance of the area is given by the 
difference, w—s, between the white and black 
excitation, and the symbols with subscripts 
represent constants evaluated so far only by the 
conditions that ¢, is greater than qd), and d, is 
greater than ¢>. 

Equation 3 is similar to those set up by Schré- 
dinger [28] in accord with the theoretical views of 
von Kries [5], and by Adams [21] in accord with 
his own theory [6]. Adams has, moreover, pointed 
out the advantages and theoretical plausibility of 
the view that the various stages are not linearly 
connected. For simplicity in the present deriva- 
tion, attention will be confined to the assumption 
represented in eq 2 and 3 that the connection is 
linear and homogeneous. 


Ill. Dichromatic Vision 


According to the Miller theory, protanopia cor- 
responds to the failure of the (yR, 6G)-chromatic 
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substance; and on this account, it is called by him 
outer red-green blindness. Thus, for protanopia, 
yR=bG=0, and we may write from eq 3: 


rp=—G,—cr B= —egY , 


Up = —b, dgY —d,rB 


Wp =@P\+6,P24 é3P3+e9gY 








sp—erB 
From eq 4 one might think that the Miiller theory 
predicts for protanopia simply the sensations of 
black and white plus the two chromatic sensations 
greenish yellow and reddish blue, since a given 
amount of blue excitation is always bound up in- 
extricably with the same minor red excitation. 
This is indeed the simplest prediction from the 
formulation and corresponds fairly well with 
reports of unilaterally protanopic observers. 
Miller, however, points out that although failure 
of the (yR, bG@) chromatic process is sufficient to 
produce the symptoms of protanopia completely, 
this failure could be accompanied by a failure of 
some of the chromatic processes in the optic nerve, 
which combination of circumstances could give 
rise to an observer having protanopic vision by all 
tests actually sensing only the hues yellow and 
blue, or even only the hues red and green. Such 
observers could be distinguished from each other 
only if one eye had trichromatic, and the other, 
protanopic vision. 

Deuteranopia, on the other hand, is ascribed by 
Miiller to failure of the (r,q)-sense of the optic 
nerve and is called inner red-green blindness. 
Thus, for deuteranopia, re=—ga=0, and we may 
write from eq 3: 


Ya= —bg=d gY +-dyR | 
Wa =¢,P,; t él. -e3P3 1 eyR t egy . (5) 
Sa —=€ bG+-e7rB | 


From eq 5 it is plain that the predicted sensations 
of deuteranopes must be black, white, vellow, and 
blue; there are no alternatives. 

Tritanopia, like protanopia, is ascribed to a 
retinal defect. It corresponds to failure of the 
(gY,rB)-chromatic substance and is called outer 
yellow-blue blindness. Thus, for tritanopia, g Y= 
rB—0, and we may write from eq 3: 








, g.=—eyR=—e,bG,  ) 








yy b,— dyR=—d,bG, 

. (6) 
w= P+ eP2+e,P3+eyR, 
a é,bG ) 


The sensations of tritanopes are seen to be pre- 
dicted as black, white, and either yellow and blue, 
or red and green, or some fixed combination such 
The latter 


hues correspond to the simplest prediction, since 


as yellowish red and bluish green. 


the (yR,6G) chromatic sensory process is un- 
affected. These hues agree well with the reports 
of tritanopes who have acquired the defect through 
a disease of the retina. 


IV. Evaluation of the Constants 


Response functions for normal and dichromatic 
vision according to the Miiller theory can be 
evaluated from eq 1, 2, and 3, for all three stages 
of excitation, provided the 14 constants of eq 2 and 
3 be evaluated. For the normal mechanism 
adapted to a stimulus yielding an achromatic 
color, the stimuli for the unitary hues, red, yellow, 
green, and blue, are known [29, 30) within limits; 
these stimuli must excite only the respective r-, 
y-, g-, and b-processes of the optic nerve. The 
stimulus yielding the achromatic color, itself, must 
cause the chromatic sensory processes yR, g), 
bG, and rB to vanish, and also reduce to zero 
the chromatic r-, 
optic nerve. 


y-, g-, and b-processes of the 
The colors confused with gray by 
the typical protanope, deuteranope, and tritanope 
must conform to the yR-process, the g-process, 
the rB-process, and their complements, respec- 
tively. The difference between the spectral lu- 
minosity function of the typical protanope and the 
same function for the normal observer must be a 
constant fraction of the (yR, 6G@)-process in order 
to conform to Miiller’s proposal; and similarly 
the difference normal and tritanopic 
luminosity must be a constant fraction of the 
(¢Y, rB)-process. 


between 


There are many more than 14 conditions to be 
satisfied, including several that are set down by 
Miller only in qualitative terms, and some that 
relate to sensibility to chromaticity differences. 
In evaluating these constants, the best determined 
conditions (marked by asterisks in the next sec- 
tions) have been satisfied perfectly; and from the 
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resulting excitation curves it may be seen to what 
extent the less well-determined conditions are 
satisfied. For example, of the data for the normal 
stimuli for the unitary hues, only those for unitary 
yellow have been used, since they were specifi- 
cally mentioned by Miiller as indicating the stimu- 
lus to be between 575 and 582 u. Unitary blue 
has necessarily to be taken as the complement of 
unitary yellow relative to the stimulus for white 
or gray; and unitary red and green were taken as 
the confusion colors for typical deuteranopia. 


1. Hueless Point and Unitary Yellow 


For an achromatic color, both the chromatic 
processes of the optic nerve and the chromatic 
cancel to zero. It is 
known [31] that stimulation of the normal eye by 
a source of equal energy results under ordinary 
conditions of observation in a closely achromatic, 
or hueless color. For such a source, Y= Y=Z, 
and from eq 1, P,=P,:—P;; hence we may write 
with sufficient accuracy from eq 2 


sensory processes must 


yR=a,P,—a,P,—0, 
whence: 


a,/d,—= 1.000; (7)* 


and similarly 


gy bP; T bP» bP. 0), 


whence we find 
b, +b. by. (8)* 


By setting r=g=y=—6=0 in eq 3, and substi- 
tuting eq 7 and 8, the expressions vanish, and no 
further relation is found. 

For the r-g cancellation point between 575 and 
582 my, we may take somewhat arbitrarily for 
simplicity the crossing point of VY and Y¥ at 578.1 
mu. For stimulation by homogeneous energy of 
this wavelength we may write 


r g=—cyR+erB=0 
¢:(a,P,—a_-P,) —¢2(b,P,+-b,P2—b,P3) =0 
(a,¢;—},e2) (3.1956.N + 2.4478 Y — 0.64342) 


(age, + bee) (— 2.5455.N +-7.0492 Y +-0.49637) +- 
b,c.(5.00007) =0 


By substituting X= )=0.4996, Z 
refer closely to 578.1 mu, we obtain 


0.0008, which 
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2.8189 (a,c, —_ bye2) — 2.2505 (are, + boc.) + 
0.0040b,c,=0. (9)* 


2. Dichromatic Copunctal Points 


The chromaticity confusions of dichromatic 
vision may be represented on the Maxwell triangle 
by families of straight lines, all lines of one family 
intersecting at a common point, known as the 
copunctal point [32]. These copunctal points have 
been evaluated from recent determinations, chiefly 
by Pitt [83], and expressed in terms of the ICI 
standard coordinate system. They embody the 
essential information regarding the chromaticity 
confusions of dichromats and lead to a convenient 
expression of four conditions affecting the unknown 
constants. 

According to the Miiller theory, protanopia 
corresponds to the failure of the (yR,bG@)-chromatic 
substance. For certain stimuli the (g¥,rB)-process 
is also reduced to zero. These stimuli are the 
equal-energy stimulus and all of those confused 
with it by the protanope. These chromaticity con- 
fusions are indicated on the Maxwell triangle by 
a straight line passing through the equal-energy 
point, and for every point on this line, gY’=—rB 
0. In particular, since all of the chromaticity con- 
fusion lines pass through a single point, these 
conditions must hold for the protanopic copunctal 
point defined by 2=X/(X+ Y¥+-Z)=0.747, y=¥ 
(V+ ¥+Z)=0.253, -=Z/(N+Y¥+Z)=0. Hence 
we may write from eq 2 

bP, +b.P,—bsP; 
b,(3.1956.N + 2.4478 Y¥—0.6434Z) 4 
b.(—2.5455 N +-7.0492 Y +0.4963Z) — 
h,(5.00007) =—0; 
whence we find 
bb. —0.0393. (10)* 
Similarly, from the less well-determined tritano- 


pic copunctal point, r=0.18, y=0.00, 2=0.82, 
we may write 


yR ae bG a,P,; _ a.P.= 0, 


a,(3.1956X + 2.4478 Y—0.6434Z) — 
a2(—2.5455.N +-7.0492 ) +-0.49637) —0; 


whence we find a,/a;=—0.0512/0.0476 1.07, 
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a value contradictory to eq 7 based upon the 
achromatic stimulus for the normal observer, 
which is much more reliably established than the 
tritanopic copunctal point. By setting a,—4a, 
in conformity with eq 7, we find that a condition 
for the tritanopic copunctal point by the Miler 
formulation is that P?,=P., which from eq 1 is 
equivalent to y=1.9872—0.327. This line inter- 
sects the z-axis of the chromaticity diagram at 
r=0.165, y=0.000, which is in as good agreement 
with the facts as the approximate evaluation of 
the copunctal point (r=0.18, y=0.00) estimated 
with the help of the Kénig theory [82]. The latter 
accords with Miller’s view that tritanopia is 
characterized by a single neutral point in the 
spectrum (near 570 my); the former places a 
second neutral point near the violet extreme of 
the spectrum (430 my) and arises from applica- 
tion of the Miller theory to the properties of the 
standard observer. As already noted, P; evalu- 
ated by this means has a small secondary maxi- 
mum in the neighborhood of 430 mu, causing it 
to cross the P,-curve at this point, a result quite 
unanticipated by Miiller. Actual reports of 
tritanopic vision are fairly well divided in this 
respect; for example, the cases reported by 
Konig [4] and Kéllner |34) vielded a single neutral 
point; those by Collin and Nagel [35] and Piper 
[36] had neutral points or areas in the neighbor- 
hood of 430 my. Willmer and Wright {14} 
found an indication of such a neutral region for 
small fields in the normal fovea, and Pitt [37] 
considers this to be typical of tritanopia. The 
difference in chromaticity between 430 my and 
the short-wave end of the spectrum is small, 
and it is possible that individual variations among 
tritanopes can account for the slight discrepancy 
in the report between no short-wave neutral point 
and one near 430 mg. On the other hand, it 
seems to be fairly frequent that tritanopes have 
ocular media pigmented heavily with brown pig- 
ment, and it is also possible that this pigmenta- 
tion would cause the spectrum to become in- 
visible to many tritanopes at a wavelength greater 
than 430 mu as in a case reported by Farnsworth 
[38]. The formulation could be made to accord 
strictly with the Miller view on this point by 
choosing constants in eq 1, so that the repre- 
sentation of P, is everywhere higher than P, 
for wavelengths less than 500 my, such as that 
found by Stiles [22], but such a wavelength dis- 
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tribution of P, is itself contrary to Miiller’s view. 
It seems more useful to proceed with an account 
of tritanopia involving a second neutral point 
near 430 my. Take therefore a;—a, as in eq 7. 

By the same argument given for protanopia, 
we may set y=—b=0, for the deuteranopic 
copunctal point defined by Y=Z=0. For eq 1, 
2, and 3 we may write 


y=—b=d,gY +dayR=d,(b,P,+b2P;—b,Ps) + 
d,(a,P,—a,P:) =(, 


= (bd +a,d2) (3.1956. X +2.4478 Y —0.6434Z) 
+ (byd, — ad) (—2.5455X 47.0492 Y+ 
0.4963Z) —byd;(5.0000Z) ; 


whence we find, since Y=Z=0, 
3.1956 (b,d, + a,d,) — 2.5455 (bd, —al2) =0. 
(11)* 


Two other conditions may be derived from the 
protanopic and deuteranopic copunctal points. 
It has been shown [32] that the deuteranopic 
copunctal point lies on the deuteranopic alychne, 
that is, the line on the Maxwell triangle associated 
with zero deuteranopic luminosity; and _ the 
protanopic copunctal point lies on the protanopic 
alychne. Hence for the respective copunctal 
points we may set (w—s)q and (w—s), equal to 
zero, and from eq 1, 2, and 3 there are found: 


€,—0.0393¢e2+-2e;(6;— 0.0393b;) =0, (12)* 
and 
3.1956 (e, + 2ega; + 2e56,) — 2.5455 
(€2—2eqa2+ 2e5b2) =0. (13)* 


3. Dichromatic Luminosity Functions 


It has been shown that the luminosity functions 
of red-green-blind observers can be expressed as 
functions of X, Y, and Z [39] by eq 14 and 15, in 
which W, is the protanopic luminosity and W,, 


the deuteranopic 
W,=0.460X+-1,359)'+-0.101Z, (14) 
W.=Y. (15) 


From eq la, these functions may be written in 
terms of the primary processes, P;, P,, and P;: 


W,=0.00754P, +-0.19017P,+-0.00229P3, 
(14a) 


W,=0.08852P, +-0.11112P,+0.00036P,. 
(15a) 


But from eq 2, 4, and 5, we may write 


W,= (w—s) p= (€:+-2b1e5) Pit (¢2+2b2¢5) P2+ 
(€;—2bses) Ps, (14b) 


We=(w—s)a= (w—8) = (e,:+-2a,¢,+2b,¢5) P+ 
(€2—2az€4+ 2b.5) P2+ (€3—2bse5) Ps. (15b) 


By equating the coefficients of P,, P2, and P; in 
eq 14a and 14b, we obtain three additional con- 
ditions to be satisfied by the constants: 


€, +26,e,=0.00754, (16) 
€2+2b,e,=0.19017, (17) 
€3— 2b,¢;=0.00229. (18) 


Similarly, by equating the coefficients in eq 15a 
and 15b we obtain three more conditions: 


€; + 2a,¢4+2b,e;=0.08852, (19) 
€3—2a.¢,+ 2b.e,=0.11112, (20)* 
eé3— 2b,¢5= 0.00036. (2 1 ) ° 


Equations 16 to 21 are not entirely independent 
of the conditions previously found, nor are they 
all congruent. From eq 16 and 17 there may be 
derived eq 12; and from eq 19 and 20 there may 
be derived eq 13. Furthermore, eq 7, 10, 12, 19, 
and 20 combine to give eq 16 and 17. Equations 
18 and 21 are contradictory. Since it is an essen- 
tial part of the Miiller theory that the deuteranopic 
luminosity function be the same as the normal, 
we must accept eq 21 and reject eq 18. This 
choice will prevent (w—s), in eq 14b from being 
exactly equal to W, in eq 14. It remains to be 
found whether the resulting evaluation of (w—s), 
is as representative of available data on the 
protanopic luminosity function as is W,. Thus, 
we have obtained only two additional independ- 
ent conditions from protanopic and deuteranopic 
luminosity functions, eq 20 and 21. 

If it be assumed for the moment that tritanopic 
luminosity is the same as deuteranopic and nor- 
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mal luminosity, as it may well be judging from 
the available information [7a, pp 53 and 102; 34, 
36, 40, 41, 42, 43], then from eq 2 and 6 we 
obtain eq 22: 


(w—s)=(w—s) = 6,P)+e.P2+63P34+-2e, (a,Pi:— 
a,P,). (22) 


By comparing the coefficients of P;, P:, and Ps 
in eq 15b and 22 we see from each of the three 
comparisons that e; must be zero; that is, there 
can be no darkening effect from the gY-process, 
such as implied by eq 3. The Miiller theory 
thus cannot abide having equality between tri- 
tanopic and normal luminosity. Miller was 
well aware that his theory required the tritanopic 
luminosity function to be different from normal 
and remarks [7a, p. 63] “In regard to spectral 
luminosity distribution in tritanopia, there must 
be, if no complications exist, because of the absence 
of the w value of the gY process, a decrease in 
the luminosity of yellow and yellowish lights in 
comparison to normal. Unfortunately 
there have been up to now no investigations of 
the spectral luminosity characteristic of tritano- 
pia.” Since the present purpose is to find the 
coordinate systems implied by the Miller theory, 
we must disregard the rather inconclusive indi- 
cations that there is no difference between tri- 
tanopic and normal luminosity; hence no atten- 
tion can be paid to eq 22 in evaluating the con- 
stants, and e, must be given a positive, though 
small value. Take arbitrarily, then: 


é5=0.03 €4. (23)* 


4. Chromaticity Sensibility and Theory 


We may now take stock of the conditions that 
must be satisfied by the 14 constants: 


Criterion Resulting condition 


Hueless point | Eq 7 and 8 
Unitary vellow- acne] BD 
Protanopic copunctal point-._ Eq 10 and 12 
Deuteranopic copunctal point Eq 11 and 13 


Dichromatie liminosity -- - -- Eq 20, 21, and 23. 


These 10 equations have been marked with 
asterisks to show that they were used in the deriva- 
tion of the constants. 
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There remain four conditions to be set up before 
the 14 constants can be evaluated. Three of the 
four conditions refer to the relative sizes of a, and 
bs, c, and d,, and a, and ¢,. 

The first ratio, a,/b;, has to do with the relative 
sensibility of the eye to yellowish red-bluish green 
differences on the one hand and greenish yellow- 
reddish blue differences on the other. Empirical 
studies on large fields [20, 21] indicate that a,/), is 
about 2.5, that is, the normal eye detects (yR, 
bG)-differences more readily by about a factor of 
2.5 than would be judged from the relatively great 
overlap of P; and P,. By setting ¢,/d;=1.0, this 
greater sensibility to (yR, bG@)-differences is pre- 
served in the optic-nerve stage. 

The ratio of a,/e, has to do with the comparative 
amounts of the chromatic sensory process and the 
chromatic excitations of the optic nerve. There 
does not seem to be any fundamental meaning to 
this comparison. It has merely to do with a rela- 
tion between the units expressing the rate of a 
chemical process in the retinal receptors and those 
expressing the frequency of the resulting impulses 
in the fibers of the optic nerve. This ratio may be 
set arbitrarily, and for simplicity we set a,/e,;—1. 

The final condition refers to the size of the 
arbitrary units in which the chromatic responses 
are to be expressed ; for simplicity take a,=1. 

Solution of these 14 equations simultaneously 
gives the values of the constants: 


a,=1.0000,a,= 1.0000, 

b, =0.0151,b.=0.3849,b,;=0.40000, 
¢,= 1.0000,c.=0.6265, 

d,=1.0000, d,=0.1622, 


€,=0.0075,¢.=0.1912,¢e,=0.0013,e,=0.0405, 
é;=0.0012 


It will be noted that, as required by eq 2, a, is 
greater than 6,, and a, is greater than b,. Further- 
more, as required by eq 3, ¢, is greater than d), 
and d; is greater than c,. This correspondence 
with the Miiller description indicates how thor- 
ough was his grasp of the facts from purely 
qualitative data, though probably c.=0.6265 is 
not as small compared to d;=1.0000 as would be 
expected from Miiller’s designation of dg)" and 




















egY as a major y-excitation and a minor g- 
excitation, respectively, resulting from the g) 
process. The value of c, would be reduced some- 
what by taking a higher value for the wavelength 
of the spectrum stimulus for unitary yellow, 
say 582 mu instead of 578 my; see eq 9. 


V. Definition of the Coordinate Systems 


We may now insert these constants into eq 
2 and 3, and so give explicit definitions of the two 
new coordinate systems implied in the Miler 
theory. The coordinate system applying to the 
chromatic sensory processes of the retinal recep- 
tors is defined by eq 2a: 


yR=—bG=P,—P,. ) 
gY = —rB=0.0151P, +0.3849P,—0.4000P, 5 


2a) 


The coordinate system applying to the precesses in 
the optic nerve fibers is defined by eq 3a; 
r g—yR + 0.6265rB bG—0.62659gY, ) 


y b=gY +0.1622yR rB—0.162266, 





w—0.0075P, +0.1912P,4 0.0013P,+ 0.0405 r (3a) 
yR+0.0012gY, 
8 = 0.040566 + 0.0012rB. : 


These two coordinate systems may also be de- 
fined in terms of the standard 1931 ICL coordinate 
system for colorimetry from eq 1. Equation 2b 
gives the definition of the colorimetric coordinate 
system for normal observers corresponding to the 
chromatic sensory processes combined with the 
luminosity function Y; and eq 2¢ gives the re- 
verse transformation from this coordinate system 
to the standard ICI system: 


yR bG 5.741.N— 4.601 Y— 1.1402, 
g)¥=—-rB 0.932.V + 2.750 —1.8192Z,) (2b) 
y 1.000); / 


N—O0.1581yR—0.0991gY i 
y ff 
Z 0.0O810yR—04991gYV + 


| (2¢) 
| 2 


Equation 3b gives the definition of the colorimet- 


’ 


ric coordinate system for normal observers corres- 
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ponding to the excitations of the optic nerve, and 
eq 3c gives the reverse transformation from this 
coordinate system to the standard ICI system: 


r= —g=—6.325X—6.325Y 
y= —b- 2.004 Y— 2.0047 (2b) 
(w—s) = 1.000" 


NV=0.1581Ir+ (w—s) 
Y= (w—s) | (3e) 
Z=—O0A9DLy + (w—s8) 


Table 1 gives the response functions of the normal, 
protanopic, deuteranopic, and tritanopic types of 
vision derived from the ICI standard observer 
according to the Miller theory (eq 1, 2b, 3b, 4, 5, 
and 6). 

The very simple transformation equations be- 
tween the chromatic excitations of the optic nerve 
according to the Miller theory and the standard 
1931 ICL coordinate system for colorimetry arise, 
of course, from the fact that the Y-primary of the 
ICI system corresponds to a stimulus for unitary 
red, and the Z-primary corresponds to unitary 
blue. The r,y (w-s) system corresponds to the 
central stage of the Adams theory [*], and coordi- 
nate systems closely resembling that described 
by eq 3b have been used by Adams with consider- 
able success to explain chromaticity spacing for 
large fields, chiefly studies of the spacing of the 
Munsell colors |*'). 

The coordinate system set up by Schrédinger [44] 
in 1925 resembles closely that defined by eq 3b 
except that it was not adjusted to correspond to 
the same balance between (7,b)-excitation and the 
(r.g)-excitation (a)/b;=2.5, e,/d)=1.0). Schouten 
[45] made use of conditions derived from the hueless 
point (eqs 7 and 8), the deuteranopic neutral point 
(eq 11) and the unitary vellow point (eq 9) to 
compute response functions for assumed central 
r-, y-, g-, and b-processes. These functions bear a 
considerable resemblance to r and y evaluated 
from eq 3b. Thus it is seen that the essence of eq 
3b is neither new nor confined to the Miller 
theory; it arises from the opponent-colors theory 
of Hering and has been used in at least three 
theoretical studies since 1925. 

Figure 6 shows as functions of wavelength 
P,, Py, Py; yR and g¥;r and y; and finally in the 
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Response functions of the normal, protanopic, deuteranopic, and tritanopic types of vision derived from the ICT 
standard observer according to the Muller theory 


TABLE 1. 


[Spectrum: equienergy] 








- ! 
Primary sensitizing Chromatic ; 
ete cam a Optic-nerve excitation 
Wave . - Gees penta =e oe a ~ | 
length Normal vision yROG | oY+R r-9 y-h w-8 
P, P» P, Protanopic vision a. | (g¥-rB)»p (*) (a) (w-8) p 
Deuteranopic vision.) (yROG)a |) @Y-rBia (y-b)a (w-s) a 
Tritanopic vision (yR-bG): (*) (*) (Be | 
- - _ — _—-+ ’ 
Di 
380 32 , -13 +8 —13 ; | 
390 l 100 +1 —40 +26 —40 0.1 0.2 | 
400 3 340 +3 —136 +88 —135 4 0.1 6 
410 8 1 1, 087 +8 —414 +267 —413 1.2 5 2.1 
420 a4 6 3, 228 +17 —1, 288 +24 —1, 286 4.0 2.5 7.0 
430 44 7 6, 928 -3 —2, 752 +1, 722 —2, 753 11.6 12 18.1 
440 45 143 8, 736 —o —3, 439 +2, 057 ~—3, 454 23 31 31.1 
410 27 292 8, 860 —24 —3, 432 +1, 886 —3,475 38 bY 46 
40 2 Sil 8, 346 — 9 —3, 142 +1, 459 —3, 224 wo 100 67 | 
470 19 TSS 6, 438 —TH4 —2, 274 +660 —2, 398 91 152 6 | 
480) 123 1,140 4,065 —1,017 —i, 185 —275 —1,350 139 220 142 
490 312 1,616 2, 326 —1,304 — 304 —1,118 —515 208 312 208 | 
} 
500 631 2,309 «1, 360 —1, 768 +389 —2,012 +102 | 323 464 sz | 
510 1, 159 3, 600 791 —2, 441 +1, O87 —3, 122 +691 503 fy SOL } 
520 1,890 4. 8N3 | —2, 993 +1, 751 — 4,090 +1, 206 70 oD ht) 706, 
530 2, 612 §, 676 2 —3, 064 +2, 140 —4, 405 +1, 643 862 1, 107 857 
540 3, 205 5, 996 12 —2, 746 +2, 316 —4, 197 +1, 871 Wt 1,174 Gas 
550 3,815 5,915 44 —2, 100 +2, 317 —3, 552 +1, 976 995 1,163 990, | 
iO 4, 333 5, 503 20 —1,170 +2, 176 —2, 533 +1, O86 was 1, 089 Qo 
570 4,704 4,772 10 —8 +1,04 —1,201 +1, 908 952 953 Y4is | 
580 5, 057 3,801 5 a +1, 255 +1, 536 +293 +1, 740 870 770 S66 
590 5, 132 2,724 6 see +2, 408 +1, 124 +1, 704 +1, 514 757 54 7H 
600 4, 938 1,745 4 +3, 194 +744 +2, 728 +1, 263 63 376 629 
610 4,435 94 2 +3, 441 +449 +3, 160 +1, 003 3 228 w2 
620 3, 663 510 1 +3, 152 +252 +2, 904 +763 Sal 129 380 | 
630 2, 702 233 2, 469 +130 +2, 387 +531 265 6s 265 
640 1, 860 w +1, 766 +4 +1,726 +351 175 “4 175 | 
650 1, 168 33 +1, 135 +30 +1, 116 +214 107 16 107 
Hino) 676 10 pees +666 +14 +657 +122 61 8 61 
670 358 3 +354 +e +350 +64 32 4 32 
fis0 191 J Srey eae +190 +3 +188 +34 17 2 17 
690 ¥3 +93 +1 +92 +16 8.2 1 Ss 
700 46 aa +46 +1 +46 +58 4.1 * 4.1 
710 24 +24 rm +23 +4 21 21 | 
720 12 +11 +i +2 1.0 9 | 
730 6 = +6 +1 5 4 | 
740 3 oe eran Maca eM LY con) +3 a 3 | 








* For protanopic and tritanopic vision the Miller theory does not state 
rigidly that the optic nerve excitation must follow eq 4 and 6, though this is 
the simplest prediction. Either r-g or y-b, but not both, may be zero. 


lower left quadrant the deuteranopic luminosity 
according to eq 3b, the protanopic luminosity 
according to eq 4 together with the luminosity 
contributions of the chromatic sensory processes 


yR and gY. 
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If they are not zero, they have wavelength distributions proportional to the 
chromatic sensory processes. 


VI. Protanopic Luminosity Function 
In the ICI system the standard luminosity 


function is represented by the second function, 
Y; and from eq 3b it may be seen that the Miller 
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Ficure 6. Response functions according to the three stages of the Miiller theory. 


Upper left: Processes in the initial photosensitive stage (same as Young theory, see eq 1 and fig. 5); lower left: components in the luminosity function (#—s) 
both for normal and deuteranopic, W4, and protanopic, W,, vision (see eq 3a); upper right: chromatic retinal sensory processes (see eq 2a); lower right: 
chromatic processes in the optic-nerve fiber stage (same as the Hering theory, see eq 3a). 


theory can be formulated, as he claimed, in such 
a way that the difference between the w-excitation 
and the s-excitation gives the normal luminosity 
function. This function has already been shown 
to be as satisfactory a representation of deuter- 
anopic luminosity as it is for some normal lu- 
minosity functions, because the deuteranopic lu- 
minosity functions fall within normal limits [32]. 
In these two respects this formulation of the 
Miiller theory conforms exactly to that previously 
worked out in accord with the Kénig theory [4, 32]. 
However, it was noted previously that eq 21 
contradicts eq 18; so it remains to be seen whether 
the prediction of protanopic luminosity by this 
formulation of the Miiller theory is as acceptable 
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as that by the Kénig theory. By inserting the 
constants in eq 4 it is found that this formulation 
of the Miller theory requires protanopic lu- 
minosity to be given by: 


(w-s),=0.0075 P,+-0.1921 P.+-0.0003 P3. (4a) 
The previous formulation of the Kénig theory 
yielded the equation: 

W,=0.0075 P,+-0.1902 P.+-0.0023 P;. (14a) 


Figure 7 is a plot of these functions adjusted 
approximately to unit maximum, together with 
upper and lower limits of available data on 


luminosity functions of protanopic and protanom- 
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Protanopic and protanomalous luminosity fune- 
tions. 


Ficure 7. 


The solid curve, Wy, corresponds to a wavelength function derived in ac- 
cord with the later Kénig form of three-components formulation for prota- 
nopic luminosity [32]; the dotted curve, (w—s),, is based on the Miiller theory. 
The arrows indicate maximum and minimum luminosities of 12 protanom- 
alous and six protanopic observers. Note that these data support both 
functions about equally well. 


alous observers [32]. It will be seen that these 
data support both functions about equally well. 


VII. Chromatic Thresholds, Normal and 


Tritanopic 


We are now in position to inquire whether the 
Miller theory offers a basis for explaining the 
chromatic-threshold data of Abney [8] and Priest 
Brickwedde [9] referred to earlier. According to 
the Miiller theory the ability of an observer to de- 
tect a slight variation in chromaticity from a central 
chromaticity, such as that of the light from a car- 
bon are, would depend upon the excitation of the 
chromatic sensory processes. The amounts of 
the excitation of these processes corresponding to 
any color specified in terms (X, Y, Z) of the 1931 
ICI standard observer can be found from eq 2b. 
And, in particular, they have been found for the 
spectrum colors for unit irradiance and are plotted 
in the upper right quadrant of figure 6. Both of 
the above sets of data are given, however, in lumi- 
nous units (luminance of the field just yielding a 
chromatic difference from carbon-are light in 
Abney’s work, or luminance fraction required to 
be mixed with sunlight to produce a color just 
noticeably different from sunlight in the work of 
Priest and Brickwedde). We should expect to 
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compare with them, therefore, the excitations of 
the chromatic sensory processes corresponding to 
the colors of a spectrum of constant luminance; 
that is, we should expect to find the chromatic 
thresholds in luminance terms to correspond to 
the reciprocal of some combination of yR/(w—s) 
and gY/(w—s). The exact form of combination 
would seem to be expressible in terms of the 
probability of a chromaticity difference being 
discriminated as a function of the probabilities of 
each of the two independent chromatic processes 
becoming effective considered separately. For the 
present purpose it is sufficient to take tentatively 
the combination as the square root of the sum of 
the squares; that is, assume, for the moment, 
that the effective chromatic excitation for large 
fields and high luminance is proportional to: 
[(yR)?+- (gY)*|'?/(w—s). For experimental con- 
ditions, such as restricted angular size of field or 
low luminance, that make the normal eye respond 
more or less like a tritanopie eye, the effective 
chromatic excitation may be assumed to be pro- 
portional to [(yR)?+-f? (gY)*]"?/(w—s), where f is 
the relative effectiveness of the gY-rB process 
compared to the yR-bG@ process. In general, we 
would compare to the chromatic thresholds 
dB/dE, expressed in luminous terms, the reciprocals 
of these assumed effective chromatic excitations 
so as to study the validity of the relatior 


dB/dE=k(w—s)/((yR) +f? GY)”, (24) 


where & is the constant required to adjust the 
theoretical function to the units in which the 
chromatic threshold is expressed. 

Abney’s data have been found to agree fairly 
well with eq 24 for k=0.0020 and f=0.04; see 
dotted curve of figure 1. The course of the experi- 
mentally determined function is followed well, 
except for wavelengths greater than 590 my 
where the predicted threshold is considerably 
lower than that found experimentally. As far as 
is known, no explanation of these data has previ- 
ously been suggested. It should be pointed out 
also that more recent determinations of the 
chromatic threshold by Purdy [46], and Otero, 
Plaza, and Casero [47] are quite at variance with 
these data, and indeed with each other. They 
show neither the sharp peak at 570 my nor the 
decline to small values near 450 mp. Needless to 
say, they are quite unexplainable by the Miller 
theory. The data by Abney and Watson, how- 
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ever, may be summarized by saying that they 
conform fairly well to the Miller theory for a 
retinal region in which the normal g¥, rB process 
is 96 percent ineffective. Tritanopia corresponds 
to complete ineffectiveness of this process. 

The data of Priest and Brickwedde have been 
found to agree well with eq 24 for k=0.05, and 
f=0.5; see dotted curve on figure 2. The degree 
of agreement is quite comparable to that obtained 
by a coordinate system adjusted empirically to 
represent such data [18]; see solid curve. In this 
case, the less complete data by Purdy are in sub- 
stantial agreement and are also shown. It should 
be pointed out, however, that these data have been 
corrected to refer to the standard luminosity func- 
tion by multiplying them by the ratio of the 
standard luminosity function to that found by 
Gibson and Tyndall. It is probable that an 
improvement in the theoretical account of other 
psychophysical data by means of the Miiller 
theory would result from revaluation in terms of 
an observer based on the Gibson-Tyndall experi- 
mental mean [18] luminosity function instead of 
on the standard observer. However, we may say 
that the Priest-Brickwedde data correspond well 
to the Miller theory for a retinal region in which 
the normal g¥-r2-process is 50 percent effective. 

It is concluded that the Miller theory affords a 
good explanation of chromatic thresholds in terms 
of a gradual approach to tritanopic vision. A 
thorough study of the implications of the Miller 
theory for chromaticity sensibility of all types 
such as that carried out by Stiles [22] for the three- 
components theory would seem to be worth while. 


VIII. Summary and Conclusion 


By taking into account the metamers known to 
be characteristic of protanopic, deuteranopic, and 
normal vision as well as data on the stimulus for a 
neutral color and the stimulus for a color of unitary 
vellow hue, the spectral variations of the responses 
for each of the three stages of the Miiller theory of 
vision have been evaluated as functions of wave- 
length. 

These response functions are shown to yield an 
account of normal, protanopic, and deuteranopic 
vision that differs in no essential respect from the 
simpler explanation yielded by the Kénig form of 
three-components theory. They differ in their 
explanation of tritanopic vision by requiring the 
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tritanopic luminosity function to be slightly 
higher in the short-wave end of the spectrum than 
normal; the three-components explanation re- 
quires it to be slightly lower in this part of the 
spectrum, 

The chromatic response functions of the second 
stage of the Miiller theory are shown to lead to a 
satisfactory and convenient account of the 
approach to tritanopia exhibited by the normal 
eye in viewing small fields or fields of luminance 
near the chromatic threshold. 

It is concluded that the qualitative ideas of 
Miiller lead to admissible and consistent coordi- 
nate systems. The Miiller theory shows how the 
three-components formulation of Young, Helm- 
holtz, and Kénig (first stage) and the opponent- 
colors formulation of Hering (third stage) may 
both be accepted, and the explaining power of both 
be simultaneously utilized. The intermediate 
stage is also a promising and powerful theoretical 
tool. The quantitative consistency of the Miiller 
ideas and the success demonstrated in accounting 
for tritanopic confusions made by normal observers 
does not, of course, prove the Miiller theory to be 
completely, or even basically, correct. Alternate 
explanations are possible. There are important 
gaps in our knowledge of retinal chemistry and 
conduction and integration of nerve impulses 
that, if filled, might disprove the Miiller theory 
and require adoption of an alternate account. 
Furthermore, several aspects of the Miiller expla- 
nation, though admissible in the present state of 
our knowledge, seem implausible and unlikely to 
be born out by future work. At the very least, 
however, the Miiller theory must be viewed as a 
forward step, and the coordinate system suggested 
by the second stage has practical value regardless 
of any of these future theoretical developments. 
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Rate of Shrinkage of Tendon Collagen: Heat, Entropy, 
and Free Energy of Activation of the Shrinkage of 
Untreated Tendon; Effect of Acid, Salt, Pickle, and 
Tannage on the Activation of Tendon Collagen 

By Charles E. Weir 





From isothermal measurements of the rate of shrinkage of tendon collagen, it is con- 
cluded that shrinkage is a rate process involving a reaction of first order. The effects of 
tension and length of soaking before shrinkage are studied. Average values of heat, entropy, 
and free energy of shrinkage of untreated tendon are obtained by application of the theory 
of absolute reaction rates to the data. These values are found to be 141 keal/mole, 349 cal/ 
mole deg, and 24.7 kcal/mole at 60° C, respectively, with standard deviations of 15 keal/mole, 
43 cal/mole deg, and 0.6 keal/mole at 60° C, respectively. Measurements of shrinkage tem- 
peratures’ are discussed. The effects of pH, salt, tanning, and pretannage treatments on 
the activation process are studied. Heat, entropy, and free energy decrease in acid and 
alkaline media. Salt solutions cause a decrease in heat and entropy but an increase in free 
energy in concentrated solutions. The effect of salt and acid in combination is complex, but 
in concentrated salt solutions an increase in free energy is observed. All tannages investi- 


gated, with the exception of chrome tanning, appear to reduce the entropy more than the 











Sr 





heat, thereby increasing the free energy. Chrome tanning increases heat, entropy, and free 


energy markedly. An interpretation of the results is given. 


. I. Introduction 


The shrinkage of leather and collagen, when 
heated in an aqueous medium, has been the 
subject of extensive investigations. The history 
of such studies is outlined adequately by Me- 
Laughlin and Theis [1]! and Hobbs [2]. The 
measurement of shrinkage temperatures consists 
of heating a water bath containing the specimen 
and observing the temperature at which initial 
shrinkage occurs. The dimensions of the test 
specimen and the rate of heating are prescribed 
(3). Theis and Schaffer [4], lacking a reliable 
analytical method for determining formaldehyde 
combined with collagen, proposed the measure- 
ment of the temperature of initial shrinkage as a 
measure of the degree of “leathering.” Recent 
investigations by Bowes and coworkers [5] have 





} Figures in brackets indicate the ijiterature references at the end of this 
paper. 
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shown that increased chrome content results in 
increased shrinkage temperature. 

Theis and coworkers |6, 7, 8] and Chater [9, 10] 
have investigated many factors that influence 
shrinkage temperature but, as far as can be ascer- 
tained, all investigations, with one exception, were 
concerned with the initial appearance of shrinkage 
under conditions of continually increasing tem- 
perature. The shrinkage temperature of chrome- 
tanned leather is generally found to be above 
100° C, and therefore it cannot be attained with 
an ordinary water bath. Chater [11], Harnly and 
Parker [12], and Merry [13] have published data 
showing that chrome-tanned leather having a 
shrinkage temperature exceeding 100° C will 
shrink progressively if the temperature of the bath 
is maintained at 100° C for a period of time. 
Chater [14] has presented evidence that the 
amount of shrinkage is a function of the initial 
temperature of the bath. Critical examination 
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of the foregoing results tends to indicate that 
shrinkage does not occur at a fixed temperature, 
although shrinkage has been regarded as analogous 
to a phase change occurring at a fixed tem- 
perature [15]. 

Recent measurements at this Bureau [16] of the 
expansivity of leather and collagen showed that 
shrinkage occurred at a temperature many degrees 
lower than that ordinarily measured as_ the 
shrinkage temperatures of the specimens tested. 
These measurements proved that the apparent 
shrinkage is a result of anistropy of the material 
and that the net change results in an increase in 
real volume. One measurement made during 
shrinkage showed that the volume change approxi- 
mately followed the first-order rate equation. 
This investigation was performed by observing 
specimens in thermal equilibrium with water over 
long periods of time. The results indicated that 
there was no fixed temperature that might be 
called the shrinkage temperature, but rather that 
shrinkage was a rate process occurring over a 
range of temperatures. In the previous investi- 
gation involving dilatometrie techniques, it was 
not possible to determine the temperature coeffi- 
cient and the related thermodynamic quantities 
of the rate process, and the present series of 
studies was initiated. 

The present investigation, utilizing tendon col- 
lagen as the experimental material, shows that the 
shrinkage phenomenon is a rate process having an 
extremely high temperature coefficient, measure- 
ments being possible over a temperature interval 
of approximately 10 to 15 deg C. The heat, 
entropy, and free energy of activation of the shrink- 
age process have been evaluated for untreated 
tendons and for tendons subjected to processes 
commonly used in the tannery. 


II. Apparatus and Method of Measure- 
ment 


The equipment required for the measurements 
consisted of a temperature-controlled water bath, 
suitable devices for suspending the specimen, and 
a cathetometer for measuring the length of the 
specimen. 

The water bath was improvised from a cylindri- 
cal glass aquarium, 12 in. in diameter and 18 in. 
in depth. The bath was equipped with a con- 
trolled intermittent heater of 500-watt capacity, 
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a variable continuous heater of 1,000-watt capac- 
ity, a thermoregulator of the toluene-mercury 
type, and an electrically driven stirrer. The 
thermostat was capable of maintaining any desired 
temperature between 30° and 95° C to within 
+0.1°C. At temperatures exceeding 65° C, the 
surface of the water was covered with a layer of 
light mineral oil to reduce excessive evaporation. 

A glass cylinder, closed at the lower end, was 
clamped in the bath in an upright position. The 
cylinder was a 500-ml graduate and was closed 
at the top with a two-hole stopper containing a 
thermometer and a glass rod. The thermometer, 
which was graduated to 0.2° C, was placed to 
record the average temperature of the specimen 
during shrinkage. The specimen was bound to 
the glass rod by fine copper wire and was placed 
under tension by the addition of a weight at the 
lower end. The weight was made up of an 
“alligater” clip to which was fastened a brass 
weight. The weight was essential to prevent 
curling of the specimen during the initial stages of 
the shrinkage. The effect of the weight will be 
described later. 

Measurements of length were obtained with a 
cathetometer graduated to 0.05 mm. Readings 
were made of the heights of the wire at the top of 
the specimen and the jaw of the “alligator” clip 
where it clamped the specimen. 

Test specimens consisted of kangaroo tail ten- 
dons, which were approximtely 1 mm in diameter. 
These tendons represent nearly pure collagen in a 
highly oriented form [17]. The tendons had been 
subjected to no chemical treatment. A confining 
sheath, which is occasionally found on some ten. 
dons, was found to restrict shrinkage and was 
therefore removed prior to test. 

Measurements were made as follows. The 
graduate was filled with distilled water, and the 
bath was adjusted to the desired temperature. 
A length of tendon (approximately 10 em and ex- 
cluding ends of tendons, which give erratic results) 
was fastened to the glass rod, anc the alligator 
clip with its appended weight was clamped to the 
bottom of the specimen. The assembly, consisting 
of rod, tendon, stopper, and thermometer, was 
then inserted into the graduate, and the stopper 
was seated securely. At the same instant a stop- 
watch was started. An initial reading of the po- 
sition of the jaw of the “alligator” clip was made 
as rapidly as possible, that is, within 5 sec. This 
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reading represents the position of the bottom of 
‘he specimen at zero time. The position of the 
‘op of the specimen was then recorded, and further 
readings were made of the position of the alligator 
clip until shrinkage ceased. The rate of shrinkage 
varies with time and temperature, and the inter- 
val between readings was varied. Numerous read- 
ings were made in the region of half-shrinkage to 
permit accurate interpolation of the data. 


Ill. Results and Discussion of Measure- 
ments on Untreated Tendon 


1. Shrinkage Curve 


Routine measurements were made on specimens 
that had been soaked in water at 5° C prior to 
testing. The results obtained on shrinkage at 
60° C of specimens soaked for varying periods of 
time are shown in figure 1. In this figure the 
residual lengths in centimeters, that is, the lengths, 
l, at time, ¢, minus the fully shrunken length, /., 
are plotted as ordinate and the time in minutes as 
the abscissa. Coordinates of this nature were 
chosen in analogy with applications of the first- 
order reaction equation. The curves illustrate the 
effect of varying the length of soaking period 
prior to shrinkage. It is observed that the dry 
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Figure 1. Effect of length of soaking on shrinkage. 


@, Soaked 0 br: D, soaked 1 hr; @, soaked 2 hr: ©, soaked 24 hr. 


specimen elongates before shrinkage occurs, prob- 
ably because of absorption of water. No pre- 
liminary elongation is observed after soaking for 
| hr or longer. Relatively little change is to be 
noted in either shape or position of the curve 
after the 1-hr soaking period. All subsequent 
measurements were made on tendons that were 
soaked in distilled water at 5° C for 24 hr prior 
to test. This procedure was adopted for con- 
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venience, since equilibrium appears to be attained 
in this period, and specimens prepared one day 
may be tested the following day. 

In figure 2, typical shrinkage curves at various 
temperatures are shown. The effect of a few 
degrees variation in temperature is pronounced, 
In general, it has been found that a 2-deg rise in 
temperature results in approximately a twofold 
increase in rate of shrinkage. This behavior 
differs from that of most ordinary chemical reac- 
tions in which a temperature rise of 10 deg C 
approximately doubles the speed of the reaction. 

The curves of figure 2 indicate that shrinkage is 
not a simple process. The shrinkage is seen to 
accelerate initially, to reach an approximately 
constant rate, and then to decelerate. As shown 
by the figure, the periods of acceleration and 
deceleration are shortened as the temperature 
increases. 











° 10 20 30 40 50 60 
TIME, MINUTES 


Figure 2. Effect of temperature on shrinkage. 
@, 58.4°C; ©, 6.1°C; B, 62.0°C 

Following complete shrinkage another process 
occurs, which is not shown in the figure, since it 
is beyond the scope of this investigation. This 
process involves reeclongation of the shrunken 
collagen and is presumably associated with solu- 
tion of the shrunken material, weakening of the 
structure, and resultant stretching. This process 
is of concern only in that it seriously limits the 
range of measurements, rendering impractical any 
measurements requiring more than a few hours. 


2. Characteristics of Shrinkage 


The mechanisnr of shrinkage, as related to the 
structure of collagen, is of great theoretical in- 
terest but is beyond the scope of this investigation. 
Observations made during shrinkage may give a 
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clue to the mechanism. <A tendon, when first 
immersed with the bath at a temperature that 
causes a moderate rate of shrinkage, is completely 
quiescent. The finely corrugated regular surface 
of the original tendon gradually becomes slightly 
distorted at fixed points, the distortion being 
caused by tiny nodules that are formed inside the 
tendon. The formation of nodules marks the 
onset of shrinkage. A nodule grows internally, 
produces a large lump, and then appears to gelat- 
inize and form a gelatinous ring around the ten- 
don. The gelatinous ring is invariably oriented 
at an angle of approximately 30° with the axis of 
the tendon and is of a much larger diameter 
(two to three times) than that of the tendon. 
Tendon adjacent to the gelatinous ring appears to 
be pulled into the ring from both directions. 
Nodules, which appear to be formed at a rate 
proportional to the temperature, grow, coalesce, 
and thus cause shrinkage. Most shrinkage is 
associated with nodule formation and growth, and 
it is most probable that the transformation of the 
adjacent tendon into gelatinous material is of 
minor importance in causing shrinkage. The 
drawings in figure 3 illustrate the course of the 
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Ficure 3. Behavior of tendon during shrinkage. 


shrinkage. These observations are remarkably 
similar to photographs made of shrinkage of fibrils 
of submicroscopic dimensions by Nutting and 
Borasky [18]. 


3. Effect of Tension on Shrinkage 


Measurements were made to determine the 
effect of tension on the experimental results. 
Shrinkage tests were made by using weights of 
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3.6, 14.4, and 40.2 g. onthe tendons. The results 
of these measurements at 65.3° C are shown in 
figure 4, measurements at other temperatures 
yielding comparable results. It is noted that in- 
creasing tension apparently causes decreased rate 
of shrinkage and an increasing induction period. 
This effect will be discussed subsequently. 
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Ficure 4. Effect of tension on shrinkage. 
D. Load 3.6 g, 65.3°C; O, load 14.4 g, 65.2°C; @, load 40.2 g, 65.4°C. 


4. Mathematical Treatment of Results 
(a) Rate Constants and Heat of Activation 


The results shown in figures 1, 2, and 4 indicate 
that, except at low temperatures, where shrinkage 
is quite slow, the shrinkage follows approximately 
an exponential function of the time after a period 
of induction. Such a relationship may be written 

l=(4,—l.)e~** +1.., (1) 
where / is the length of tendon at time ¢, /. is the 
completely shrunken length, J is the initial length, 
and K is the rate constant for the shrinkage. 
Equation 1 is of the form of a first-order equation 
[19]. If such a process is investigated, a graph of 
log (/—/..) as ordinate and ¢ as abscissa will yield 
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Ficure 5. First-order plot of shrinkage data. 
O, Load 3.6 g, 65.3°C; O, load 14.4 g, 65.2°C; @, load 40.2 g, 65.4°C. 


a straight line of slope —K. Figure 5 shows the 
results given in figure 4 plotted in this manner, 
and it is seen that following a period of non- 
conformity, approximate linearity is obtained. 

Objections to the use of rate constants obtained 
in this manner are as follows: (1) Graphs of data 
at lower temperatures show negligible linearity 
and, therefore, do not yield values of K. This 
behavior is probably related to the fact that a 
physical measurement of length is being utilized 
to follow a complex chemical process that probably 
involves consecutive reactions. (2) At higher 
temperatures adequate linearity is attained, but 
shrinkage occurs so rapidly that it is difficult to 
obtain sufficient data to permit accurate evalua- 
tion of K. Therefore a more useful measure of 
the rate of the reaction has been utilized, the time 
of half-shrinkage, designated by t:2 and defined as 
the time, ft, at which the length of the tendon, J, is 
equal to (4)+/.)/2. 

Let the energy, or heat, of the activation pro- 
cess be AH7. Then the temperature dependence of 
KX [19] is given by the equation 


~~ 


where in addition to previously defined terms, R 
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is the gas constant, 7’ is the absolute temperature, 
and C is a constant. A similar relationship holds 
if In(1/t*) is substituted in eq 2 for In A, where t* 
is defined to be any corresponding time in the 
course of the reaction. The corresponding time 
is taken to be a measure of the extent of the 
reaction, and in the present study the time corre- 
sponding to half-shrinkage t, is used. Since the 
evidence that has been presented indicates that 
the shrinkage follows, in part, the law of a first- 
order reaction and, since it has been found that in 
general protein denaturation follows this law [21], 
it is assumed that the shrinkage process is a dena- 
turation reaction of first order and that the times 
of half-shrinkage obtained from experimental data 
represent times of half-chemical reaction. It is 
seen from eq 2 that the latter portion of the 
assumption is unnecessary to obtain correct 
values of A//. It will be required, however, in 
the theory of absolute reaction rates to obtain 
entropy and free-energy values. 

The validity of use of t,, is justified by linear- 
ity of graphs obtained by plotting log (1/ty) as 
ordinate and 1/7 as abscissa. The linearity ob- 
tained implies either that t, is a true measure 
of a corresponding point in the reaction or that 
AH and ty vary in a complicated manner to 
produce linearity. The latter possibility is remote 
and therefore rejected as improbable. Values of 
ti, may be obtained readily by interpolation of 
data and will be used hereafter to characterize 
the reaction velocity. Values of AH will be 
obtained as the slopes of lines obtained from 
graphs of log (1/t,,) and 1/7 as described later. 


(b) Entropy and Free Energy of Activation 


The shrinkage process can be measured only 
over a short temperature range, which implies a 
large heat of activation. On the basis of the 
theory of absolute reaction rates [20], in order for 
such a process to occur with measurable velocity 
at ordinary temperatures, a large entropy increase, 
AS, must be associated with activation. In other 
words, the free energy of activation, AF’, defined 
by the relationship, 


AF=AH— TAS, (3) 


must be of a much smaller magnitude than AZZ. 
To obtain these quantities from the experimental 
data, use was made of the theory of absolute 
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reaction rates. The fundamental equation of 
this theory is 


K ATK’, (4) 


where A is the reaction velocity constant, &, 
Boltzmann’s constant, A, Planck’s constant, 7, 
the absolute temperature, A’, the equilibrium con- 
stant of the reaction involving formation of an 
“activated complex,’ and «, the transmission 
coefficient. This factor will be constant under 
the conditions of this experiment, and in the 
absence of information concerning its magnitude, 
will be assumed to be unity. As this reaction 
involves disorganization of the structure of 
collagen and therefore may be considered to be a 
decomposition, « probably is approximately 1 
[20]. If the term in « is included it will represent 
an uncertainty in AS but will have no effect on 
AIT, as will be seen by a consideration of later 
equations. 

Substituting for AK’ in eq 4 from the relationship 


AF=—RT In kK’ (5) 


we obtain an equation relating A and F, 


J" -SF/RT 
K re ‘ (6) 
Substituting for A F from eq 3, 
7 —(QH/RT-aS/R) 
K s ’ (7) 


h* 


To obtain a relationship involving f,, in place of 
K for reasons previously enumerated, assume the 
equation 

- 0.6932 ' 
K . (8) 
t, 

ry’ ‘¢ >) “) . », 

The constant 0.6932 corresponds to the numer- 
ical relationship between AC and f,, for first-order 
reactions [19]. Substituting into eq 7 and re- 
arranging, the following relationships are obtained: 

0.6932h AF /RT (SH /RT—SS/R) 


kT, ~~ ; (9) 


By taking logarithms of these equations, it is 
seen that a graph of log 0.6932A/k7t, and 1/T 
vields a line that has a slope and intercept pro- 
portional to A// and AS, respectively. 
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Solving eq 9 and 7 for AS the two following 
results are obtained: 


AS=2.303R [ toe Uptlog j+log 0.6932 [+7 


(10) 

AS=2.303R [toe Lptlog K]+7y" (11) 
Equation 11 will vield the absolute value of AS 
for the activation. The essential difference be- 
between eq 10 and 11 lies in the proportionality 
factor 0.6932 of eq 10. Assuming the value of 
0.6932 is in error by a factor of 1,000, an uneer- 
tainty of only 13.8 cal/mole deg in AS results. 
A corresponding uncertainty of only 4.6 kcal/mole 
in values of AF at 60° C is obtained. Derived 
values of A/7 will be unaffected by uncertainties in 
the relationship between A and ty. It seems, 
therefore, that AS and AF obtained through use 
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Figure 6. Application of absolute reaction rate theory to 
data. 
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of ty values may be considered to be good ap- 
proximations to the absolute values of AS and AF. 

It will be noted that the transmission coefficient, 
x, would appear in an added logarithm term in eq 
10 and 11, and unless « differs appreciably from 
unity, would have little effect on AS. 

The results obtained under three different ten- 
sions are shown in figure 6, graphed in accordance 
with eq 9, log 0.6932h/kTt,, being plotted as ordi- 
nate and 1/7 asabscissa. In addition, the reaction 
velocity constants of the measurements with the 
14.4-g weight were obtained graphically and are 
plotted in accordance with eq 7. In the latter 
graph only three points are obtained, since no 
linearity and no velocity constant were obtained 
at the lowest temperature. The results for A// 
and AS shown in table 1 were obtained by appli- 
cation of least squares to the data shown in figure 
6 and are considered to be in fair agreement. In 
practice it was observed that use of the largest 
stress greatly reduces the temperature range ‘in 
which linearity may be found in a first-order 
application. Consequently, it is believed that the 
apparent deviation from the first-order reaction 
is in some manner related to initial stresses existing 
in the tendon from internal and external forces as 
well as to the method of measuring the course of 
the reaction through the change in length. 


TaABLe 1. Effect of tension on activation constants 


f Free energy 


Entropy o of activation 


Data used Heat of acti- 


Load a. ae vation. AH —e at Cc, 
AF w 

g keal’mole —cal/mole deg —— kcal/mole 
$6 tis 136 330 26. 2 
44 fy 136 334 24.8 
44 K 113 27 20.8 
40.2 ths 14 385 25.8 


The results of figure 6 indicate a trend toward 
higher temperatures under larger tension. The 
increase in shrinkage temperature under tension 
has been reported by Woblisch [15], but the experi- 
mental variation is such that the differences noted 
in figure 6 may be due to sampling. 


5. Heat, Entropy, and Free Energy of Activation 
of Untreated Tendons 


The results of measurements of thermodynamic 
activation constants of various batches of tendons 
in distilled water (pH 6 to 7) are shown in table 2. 
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Values of A7/ and AS were obtained by application 
of the least squares to data consisting of f,, and 
1/T values. Values of AF) were obtained by use of 
results for A//7 and AS and application of eq 3 and 
are calculated at 60° C. This temperature was 
chosen since rapid shrinkage of tendons occurs at 
60° C. All measurements were made using the 
14.4-g weight. 


TABLE 2. Activation constants of tendon collagen 


Entropy of Free energy of 
Heat of activation, AH activation, activation at 
As 60° C, AF « 


keal mole cal/mole deg keal/ mole 
136 334 25.0 
136 330 25.2 
145 Sol 25.1 
168 431 24.7 
152 383 24.8 
141 34s 24.9 
1 292 23.2 
128 313 244 
Average i141 3ay 24.7 
Standard deviation — 15 43 0.6 


It is found that, although AH varies between 
120 to 170 keal/mole and AS between 290 to 430 
cal/mole deg., the values of AF remain substan- 
tially constant at approximately 25 keal/mole. 
The size of the values of AH and AS are indica- 
tive of chemical reactions involving large molecules 
and agree in magnitude with similar values previ- 
ously obtained for denaturation of soluble pro- 
teins [21]. The heat of activation is probably 
identified with rupture of bonds during activation. 
The entropy change is probably associated with a 
marked disorientation occurring in the process of 
activation. The variation observed in these 
values may be due to varying stress caused by 
different cross-sectional areas, as indicated by 
figure 6, to varying degrees of crystallinity or 
bonding in different specimens, or to experimental 
errors inherent in determination of quantities of 
such magnitude. 

The process of shrinkage may be assumed to be 
represented by the following reactions: 


Collagen A@Collagen B->Collagen C. 
(inactive) (activated (shrunk 
complex 


The heat, entropy, and free energies obtained re- 
late to the equilibrium process on the left. A 
measure of the position of this equilibrium is given 
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by AF, as shown by eq 5, which relates AF and 
the equilibrium constant of this process. The 
value of AF is then considered to be a measure of 
the temperature interval of shrinkage, as the con- 
centration of collagen B will be determined by the 
value of AF in conformity with eq 5. 


6. Measurement of Shrinkage Temperatures 


The foregoing results show that a true shrinkage 
temperature does not exist, since the process in- 
volved is a rate phenomenon. The high heat of 
activation, however, is probably responsible for 
the fact that shrinkage has heretofore not been 
recognized as a rate process and to a large extent 
minimizes errors in measurements of shrinkage 
temperatures. However, interpretation of shrink- 
age as a process occurring at a fixed temperature 
is subject to all errors inherent in interpreting a 
rate phenomenon as a fixed temperature process, 
plus an additional serious error related to the fact 
that the rate of shrinkage initially accelerates 
from zero. 

By defining the end of the acceleration period as 
that time at which adherence to the first-order 
equation begins, it has been ascertained that it has 
a temperature coefficient approximately equal to 
that of the shrinkage reaction, that is, AH~130 
keal/mole. This means that the initial process is 
actually a portion of the shrinkage, and it will 
invariably cause a shrinkage temperature meas- 
urement to yield too high a temperature, that is, 
a positive error. 

Specimens contained in a water bath heated at 
the rate of 2 deg C/min will necessarily remain 
at a temperature lower than that of the water. 
This consideration also leads to a positive error. 

Dry specimens, when immersed in the bath 
require much longer periods of time before initial 
shrinkage than specimens that are soaked prior 
to test. A shrinkage test specimen placed in a 
bath and tested immediately may give a high 
result. 

All experimental results indicate positive errors 
in measurements of shrinkage temperatures, and 
the results obtained from these measurements 
are doubtless too high [16]. For ordinary test 
purposes, it is believed that the present method 
[3] is adequate, since it is very doubtful if the 
error can exceed 5 deg C. For precise work, 
however, the method leaves much to be desired. 
In such instances it should be possible to define 
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shrinkage temperature in either one of two ways: 
(1) prescribe a permissible change in surface area 
of a specimen of known area subjected to a fixed 
temperature for a specified length of time, (2) de- 
fine the shrinkage temperature as that tempera- 
ture causing half-shrinkage in 1 min. The first 
method of definition would require no apparatus 
other than that normally found in any laboratory. 
The latter, however, would necessitate rate meas- 
urements and is less preferable experimentally, 
although more desirable thermodynamically since 
it, in effect, specifies AF and therefore K’ (eq 5). 


IV. Results of Measurements on Treated 
Tendon 


1. Effect of pH of Equilibrium Solution 


Initial attempts to evaluate thermodynamic 
activation constants of tanned tendons led to 
eyratic and inexplicable results. As tanning in- 
volved occasional use of extreme pH values, it 
was decided to ascertain the effect of pH on the 
activation. In these experiments, the low pH 
values were obtained with sulfuric acid and the 
high pH values with sodium hydroxide. All 
specimens were allowed to equilibrate in the 
solution of desired pH for 24 hr prior to test, and 
the measurements were made by using the equi- 
librium solution as the confining liquid. A large 
excess of liquid was used in every test. Extreme 
acid and alkaline pH values caused enormous 
swelling of the tendons, and early results showed 
the advisability of subjecting such specimens to 
refrigeration at all times prior to testing. Results 
of these tests are shown in table 3 and figure 7. 
Table 3 shows that heat, entropy, and free energy 
of activation decrease in acid and alkaline solu- 
tions. Attention is called to the two latter values 
of the table. These values represent results ob- 
tained on tendons that had been subjected to the 
extreme pH values and then leached with water 
until all acid or alkali was removed. It is ob- 
served that substantial recovery of the initial 
state has been obtained. This result is important 
since it justifies the use of acid or alkali in the 
tanning process, provided thorough washing is 
effected. 

The nature of the shrinkage of highly swollen 
tendons differs from that of normal tendons in 
that the acceleration period is negligible and possi- 
bly entirely missing. There is considerable 


Journal of Research 


4 








yy 











180 T 














120 


OF KCAL/MOL 























60 








oe T T ‘ar et i ae 


OS,CAL/MOL DEG 








PH 


Figure 7. Effect of pH on heat and entropy of activation. 
O, AF; @, AS. 


deviation from linearity under these conditions 
in the graph of log 0.6932h/k7t,, versus 1/T in the 
lower temperature ranges. It appears probable 
that linearity is the limiting condition occurring 
at high temperatures. In figure 8 the linear plots 
are shown from which the data of table 3 were 
obtained. 


TaBLe 3. Effect of pH on activation constants of tendon 
collagen 





Free energy 


Heat of Entropy of of activation 


pil activation, | activation, 
AH 4s 


keal/mole | cal/mole deg kcal/mole 


18 67 150 | 17.4 
3.0 ‘ 79 175 20.6 
3.6 l4l 363 19.8 
4.2. 145 368 22.5 
6.8 : 168 431 24.7 
7.2 152 383 24.8 
10.8 _. . 132 327 23.5 
12.5 77 177 18. 4 
1.8 leached with water to 7.2 142 358 23.3 
12.5 leached with water to 7.3. 125 305 23.0 


Figure 8 shows indirectly that as the pH varies, 
the range of temperatures in which measurements 
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Ficure 8. Effect of pH on rate of shrinkage. 
@, pH 1.8; 0, pH 3.0; ©, pH 3.6; @, pH 4.2; O, pH 6.8; ©, pH 7.2; ©, pH 
10.8; @, PH 12.5 
can be made, that is, the shrinkage “‘temperature”’ 
hereafter denoted as JT, and defined as_ that 
temperature causing half-shrinkage in 1 minute, 
varies in direct proportion to AFyo. The variation 
of T, with p#/ is shown by the shift of the linear 
plots parallel to the 1/7 axis. This finding is in 
agreement with observations made by Theis [6]. 
The dependence of the activation constants on 
pH implies that activation involves reaction with 
or loss of hydrogen ions [20]. The neglect of the 
role of H*, implicit in application of the absolute 
reaction rate theory without considering the re- 
action of H*, is equivalent to choice of variable 
standard states as regards hydrogen-ion concentra- 
tion. Further measurements on the effect of 
pH on activation may be expected to yield valu- 
able information concerning the nature of the 
activation process. 
2. Effect of Sodium Chloride 
Since salt is widely used in the tanning process 
as a swelling deterrent, it is of interest to ascertain 
the effect of salt on the activation constants. 
In these tests the pH was maintained at approxi- 
mately 6.2, adjustments being made with dilute 
sulfuric acid when necessary. Solutions were 
made by weighing the desired quantity of chemi- 
cally pure sodium chloride (NaCl). Tendons 
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were soaked in the solutions for 24 hr prior to test 
and tested with the equilibrium solution as the 
confining liquid. Results of these tests are shown 
in table 4 and figure 9. 
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Ficure Oo Effect of sodium chloride on rate of shrinkage, 
pH 6.3. 
O05 M: @, 0.1 M, @ 10M; @, 2.0 M; @, 4.0 M; @, saturated. 

Tapie 4. Effect ef sodium chloride on activation of tendon 


collagen 


Free energy 


Conentrstimgt | Ackttys | sctva:_ | activation, | ©/gratoaton 

. tion, All AS AFs " 

mole liter keal mole cal/mole deq keal mole 
ono “0 141 wy 24.7 
os “3 114 Was 23.5 
1 fi. 2 WS 217 23.2 
oO 6.3 us 218 22. 4 
2.0 6.2 SS 176 24. 1 
40 f.1 ST 185 25.8 
Saturated 60 75 My 25.3 


that A/7 and AS 
decrease with increased concentration of sodium 
chloride with one exception, that is, 4 M/ sedium 
chloride. This well 
within the limits of variation noted for untreated 
tendon. The value of AF remains substantially 
constant but appears to proceed through a mini- 
mum in the region of 1 \/ sodium chloride. The 
variation of AFyy is beyond the experimental 
variation and is probably real. In figure 9 it is 


It is observed in table 4 


exception, however, — lies 
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noted that 7, likewise reaches a minimum in this 


region of concentration of sodium chloride. This 


result corroborates findings by Theis [6]. 


3. Effect of Pickling 


The effect of a combination of acid and sodium 
chloride was determined in a similar manner. One 
concentration of acid was used, and the salt 
content was varied. Sulfuric acid of pH 1.8 was 
used because of the low activation constants 
observed at this pH and the fact that maximum 
swelling of collagen occurs in this region of pH, 
as shown by Highberger [22]. Tendons were 
soaked in the solutions for 24 hr and tested using 
the equilibrium solution as the confining liquid. 
Desired concentrations of sodium chloride were 
obtained by weighing chemically pure NaCl. 
The results of these tests are given in table 5 and 
figure 10. 
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Figure 10. Effect of sodium chloride on rate of shrinkage, 
pH 7.8. 


dD, 0.1 M; @, 0.5 M; OO, 1.0 M; @, 2.0 M; ©, 4.0 M; @, saturated. 


Tasie 5. Effect of acid and salt on activation constants of 


tendon collagen 


. eat of entropy of ‘ree ener 
pe nny ae hag — B.. = tal ~ tives ion, on tiv ition 
All AS wEerc, AFe 

Mole/liter keal, mole cal/mole deg keal mole 
00 1.8 67.0 10 17.4 
1 Ls $1.2 Mi 2 12.4 
5 1s my “oY 16.3 
10 1S 12h 10 23.3 
20 t% 161 405 26.2 
0 1S 18 oat 2.9 
Saturated 1S we 212 20.1 


It is of interest to compare these values with the 
results previously obtained with sulfuric acid and 
sodium chloride individually. Bearing in mind the 
experimental variation previously obtained with 
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untreated tendons, the general trends with increas- 
ing concentration of sodium chloride appear to be 
as follows: (1) A// increases, passes through a 
maximum, and then decreases, (2) AS follows a 
similar behavior, (3) AF%o appears to increase 
monotonically. Comparing these results with 
those obtained previously, it is seen that the 
behavior of the combination of sodium chloride 
plus sulfuric acid toward collagen is not a simple 
additive property of the individual effects. The 
results indicate that pickling is complex in nature 
and possibly involves two competitive reactions 
involving collagen, sulfuric acid, and sodium 
chloride. 

It is known that acid and salt individually cause 
swelling of collagen [22, 23]. It was observed that 
the swelling caused by 0.1 MM sodium chloride in 
sulfuric acid of pH 1.8 was far greater than that 
resulting from sulfuric acid of pH 1.8 alone. This 
is probably the cause of the low constants found 
for the 0.1 MM sodium chloride solution of pH 1.8. 


4. Effect of Tanning 


Tanning of the tendons has presented difficul- 
ties, particularly in view of the fact that extraneous 
masking salts are undesirable because they may 
introduce unknown effects. Tanning was effected 
in most instances by immersing the specimens in a 
solution of the tanning agent and adjusting the pH 
to the optimum value. Following tannage, speci- 
mens were thoroughly washed and tested by using 
distilled water as the confining liquid. Some speci- 
mens were reserved for analysis, consisting in most 
instances of an ash determination. Tannage with 
iron and aluminum was found to necessitate mask- 
ing salts and, to date, no satisfactory vegetable- 
tanned tendons have been prepared. This does not 
imply that vegetable tannage does not occur, but 
for purposes of test a uniform tannage must be ob- 
tained, and the uniformity presents difficulties. 

Tannages were secured by the following 
methods: formaldehyde: 15-percent solution of 
formaldehyde, pH 6 to 7 [24]. Vanadium: satu- 
rated solution of sodium vanadate, pH 3 [25]. 
Iron: solution of ferric sulfate and citric acid, pH 5 
{26}. Aluminum: solution of aluminum sulfate and 
sodium citrate, pH 5 [27]. Uranium: solution of 
sodium uranate (complex carbonate formed), pH 5 
|28]. Zirconium: solution of zirconyl sulfate, pH 2 
[29]. Chromium: two-bath process consisting of 
immersion in a solution of potassium dichromate 
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followed by reduction with a solution of sodium 
bisulfite [1]. Bating: trypsin treatment in satu- 
rated solution of calcium carbonate at 35° C for 
24 hr. 

The results of tests of tendons subjected to these 
tannages and treatments are given in table 6. 


Tape 6. Effect of tanning and pretannage treatments on 
activation constants of tendon collagen 


Heat of | Entropy of Free energy 


Tannage or treatment Sy activation, activation, of activation 
aT as as at 60°C. AF ae 
Percent keal/mole cal/mole deg | kcal mole 

Bated 161 40s 24.8 
Thiosulfate 131 31s 24.6 
Bisulfite 14 421 23.7 
Dichromate 129 2st 23.7 
Formaldehyde 113 252 29.0 
Zirconium 7.0 137 wT 35.1 
Iron 5.3 lly 282 25.4 
Vanadium 12.3 118 pa) | 25, 2 
Aluminum 1.5 82.1 72 24.8 
Uranium 3.0 102 231 25.3 
Chromium 0.06 (CrO,) 1H 404 26.5 
Do 13 (CrgQs) 208 522 5 

Do .57 (CrxOs) 323 S62 36.1 

Do 1.02 (CreOy3) 390 1,041 4.1 

Do 1.98 (CryO,) 325 S41 45.2 


Three important conclusions are to be drawn 
from the results of table 6 as follows: (1) Chrome 
tanning alone drastically increases A/7, AS, and 
AF. Further, the amount of chrome required to 
produce this effect is small, that is, about 0.6 
percent of chromium oxide (Cr,Q;).. (Chromium 
contents were determined by perchloric acid and 
by colorimetric methods [3, 30]). (2) Tanning 
with zirconium and treating with sodium di- 
chromate, sodium bisulfite, sodium thiosulfate, 
and trypsin have no appreciable effect on AZ7, and, 
of these, zirconium alone increases AF, an increase 
resulting by decrease of AS. (3) Tanning with 
formaldehyde, vanadium, iron, aluminum, and 
uranium appears to decrease A/7 but, except for 
formaldehyde, does not affect AF yo appreciably. 
In some instances the decrease in A/7 is doubtful 
since it is small. 

Tanning with most metallic tanning agents in- 
creased the ash content only slightly, and in most 
instances probably mediocre tannage was obtained. 
However, in view of the extraordinary changes 
caused by 0.6 percent of Cr,Oy, it is seen that other 
metallic tans differ from chromium. 

The graphical results from which A/7 and AS 
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Figure 11. Effect of tannage on rate of shrinkage. 
@. Alum tanned; @, iron tanned; ©, zirconium tanned; O, formaldehyde 
tanned; @, chrome tanned, 0.13 percent CROs; @, chrome tanned, 0.57 
percent C R,0;; @, chrome tanned, 1.03 percent C R20s. 


were determined are shown for some tannages in 


figure 11. 
V. Discussion 
1. Hydrothermal Stability of Leather and Collagen 


Shrinkage is a rate process and, as such, presum- 
ably proceeds at all temperatures at a rate gov- 
erned by the temperature, previous treatment, and 
environment as shown by the foregoing results. 
To obtain an idea of the behavior at room tem- 
perature a graph was made in which the logarithm 
of the time of half shrinkage in minutes, log ¢,,, is 
plotted as ordinate and the reciprocal of the 
absolute temperature, 1/7, as abscissa. The 
measurements are limited to a short temperature 
interval by virtue of the large value of A// but, 
assuming the linearity is valid beyond this range, 
it is possible to extrapolate to room temperature 
and obtain the time of half shrinkage at such 
temperature. The graphs are shown in figure 12. 
Because of the nature of the slopes and intercepts, 
it is not possible to show all extrapolated curves 
on the same graph, but the extrapolations of all 
results shown in figure 12 yield the approximate 

25° C shown in table 7. 


values of t,, at 25 
The results of table 7 must be considered 


approximations if only because of the length of the 
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Ficure 12. Effect of tannage on half-shrinkage time. 


@, Zirconium tanned; ©, formaldehyde tanned; ©, chrome tanned; 
©, alum tanned; 0, untreated, pH 7.2; @, sulfuric acid, pH 1.8; , sodium 
hydroxide, pH 12.4. 





necessary extrapolation, but, even as approxima- 
tions, interesting conclusions result. It is noted 
that extreme pH values cause fairly rapid shrink- 
age, that is, small time of half-shrinkage. It is 
for this reason that specimens in alkaline and 
acid solutions were cooled as previously described. 
Tanning, with the exception of alum-tanning, 
drastically increases the time of half-shrinkage at 
room temperature. This increase is due to two 
factors: (1) A shift of the line to lower values of 
1/7, that is, increased 7,, and (2) an increase in the 
slope of the line, that is, increase in AH. This 
latter behavior is noted only in the case of chrome- 
tanned tendon. A comparison of the relative 
efficiency of these two processes in increasing ty, 
is afforded by comparing the values of ty, for 
chrome and zirconium tannages. The zirconium 
tannage caused the largest shift in 7., whereas 
the chrome tannage increased AH. It is noted 
from table 7 that the latter effect is many fold 
more effective in increasing ty. To date, however, 
no treatment other than chrome tanning has been 
observed to produce an increase in A//. 

TABLE 7. Time of half-shrinkage at 25°C 





Treatment Time of half- 





shrinkage, ts, 
min 
Chrome-tanned. 2x10 
Formaldehyde-tanned 3x10" 
Zirconium-tanned 5x 10% 


Alum-tanned 8x 10° 


Untanned, pH 7.2 5x10" 
Untanned, pH 1.8 15x10 
Untanned, pH 12.4 8x<10° 


Recent results obtained in this laboratory by 
Kanagy [31] indicate that shrinkage at slightly 
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elevated temperatures and high relative humidity 
may be an important factor in the deterioration 
of leather. From this point of view, it would 
appear that the most stable material would result 
from chrome-tanned leather having the highest 
possible shrinkage temperature. 

The results of table 7 show that the time of 
half-shrinkage is drastically reduced when tendons 
are subjected to extreme values of pH. Assuming 
that hide behaves in a similar manner, it is con- 
cluded that exposure of untanned skins to extreme 
values of pH may cause irreparable damage to 
the finished leather, unless a swelling deterrent 
such as salt is present. 

2. Interpretation of Heat, Entropy, and Free 
Energy of Activation 
(a) General Considerations 

Theories of the structure of collagen and the 
tanning process have been published by many 
authors. An excellent bibliography is given by 
McLaughlin and Theis [1]. Essentially, the pres- 
ent picture of collagen assumes three types of 
forces: (1) valence forces along the peptide chains, 
(2) salt bonds (ionic linkages) between carboxyl 
and amino groups occurring in side chains of 
acidic and basic amino acids, and (3) hydrogen 
bonds between carbonyl groups and neighboring 
amino or imino groups. The particular theory of 
tanning of interest here is the “bridging’’ theory, 
the origin of which is in controversy [32]. This 
theory depicts tanning as a reaction in which a 
tannin molecule reacts with portions of two ad- 
jacent peptide chains to form a “bridge.” The 
resultant increase in shrinkage temperature is 
conceived as due to the increased thermostability 
conferred by the cross-linking tannin molecules. 
Shrinkage is then postulated as due to rupture of 
cross-linkages by thermal agitation and a subse- 
quent collapse of the polypeptide gridwork [1]. 
The evidence for the rupture of linkages during 
shrinkage has been given by Gustavson [33] 
and Theis [34]. 

Gelatinization involves rupture of few peptide 
linkages [35], whereas the shrunken collagen 
possesses a larger number of active groups than 
the original material [33, 34]. Therefore it is 
probable that the activation process consists of 
rupture of key cross-linkages in the collagen 
molecule. On the assumption that the same 
cross-linkages must be broken irrespective of the 
initial configuration, that is, the structure of the 
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activated complex is essentially the same, regardless 
of tannage, environment, etc., it is possible to obtain 
interesting information from the foregoing results. 

Superficially, the results show that the activa- 
tion reaction (eq 12) involves absorption of a 
tremendous quantity of heat. The total heat ab- 
sorbed in the transformation of collagen A into 
collagen C has been measured calorimetrically 
and is of the order of 12.5 cal/g [36]. On the 
assumption that no heat is involved in the irre- 
versible process (collagen Be>collagen C), this 
figure combined with the value of heat of activa- 
tion 141 kcal/mole gives a molecular weight of 
11,300 for collagen, a figure of the order of mag- 
nitude expected for a protein. However, the 
assumption involved in the calculation is most 
probably not valid. It appears that no further 
progress along these lines is possible without ex- 
tensive calorimetric investigations. 

The size of the heat of activation is not dis- 
turbing, since large molecules are involved in this 
reaction. As shown above, an assumption of a 
molecular weight of 11,300 yields the correct 
value for the heat of shrinkage. Inasmuch as the 
molecule may be much larger, the values of heat 
and entropy involved are not excessive. 


(b) Heat of Activation 


The linkages ruptured in activation are probably 
either salt or hydrogen bonds. Consider the 
following facts about AH: (1) A/T decreases with 
increasing concentration of sulfuric acid (to pH 
1.8); (2) AH7 decreases with increasing concentra- 
tion of sodium chloride, even in concentrated solu- 
tions; (3) AH decreases with increasing concentra- 
tion of sodium hydroxide (to pH 12.5). 

Since electrolytes may be expected to act pri- 
marily at ionic linkages and AH appears to de- 
crease in acid, base, and salt solutions, the effect 
may be attributed to a salt-type interaction. 
From this consideration plus the fact that a large 
number of hydrogen bonds must be postulated to 
account for the magnitude of the AH value ob- 
served, it is believed that A/7 measures for the 
most part the strength of the salt linkages in the 
collagen molecule. 

The results obtained may all be explained by 
interaction of solutes at salt bonds combined with 
swelling of the protein structure. However, the 
effects of electrolytes on hydrogen bonding might 
conceivably cause similar effects. In the light of 
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the results obtained with sodium chloride, it is 
believed preferable to conclude that salt bonds 
rather than hydrogen bonds are primarily respon- 
sible for the variations observed in AH, and 
further that salt bonds contribute at least 70 kcal 
to the total 140 keal observed for AH. 

If the bridging theory of tanning is correct, it 
is to be expected that A// will increase on tanning. 
Chromium is the only material investigated that 
causes such increase in A/7. Other tanning mate- 
rials, such as iron, uranium, vanadium, and 
formaldehyde, appear to decrease AH, although 
the variation is so small as to be doubtful in view 
of the experimental spread found previously. 
Aluminum decreases A// and, therefore, behaves 
as a salt interacting with the salt bond. 

Chromium drastically increases AH/, and the 
increase may be pictured as due to either strength- 
ening existing cross-linkages or formation of new 
linkages. 

Therefore it is believed, on the basis of the pres- 
ent experimental evidence, that A// measures 
the salt-bond cohesion energy and that chromium 
alone, of the materials studied, can form any cross- 
links between polypeptide chains. Highberger 
has presented evidence of a mechanical nature 
that formaldehyde does not form cross-linkages 
[37]. 

Inspection of results obtained with chrome- 
tanned tendons shows that about 1.0 percent of 
chromium oxide is sufficient to produce the appar- 
ent maximum cross-linking and that further in- 
crease in chrome content produces small effect. 
This amount of chromium is far less than the 
amount required to react with the excess acidic 
and basic groups present in collagen [35]. There- 
fore, it is necessary to conclude that only a fraction 
of these groups present possess the required spatial 
orientation to either form salt-bonds, which are 
strengthened by chromium, or react with chro- 
mium itself to form new cross-linkages. The 
amount of chromium combined in excess of 1 per- 
cent probably combines with other groups of the 
collagen in the same manner as the other tanning 
materials studied, that is, affects AS. 


(c) Entropy of Activation 


The entropy of activation is identified with the 
degree of disorganization or disorientation occur- 
ring in the activation process. This term will, in 
effect, measure the crystallinity of the collagen, 


30 





prior to activation, assuming similar final states 


as regards collagen B. Thus the addition of 1 
percent of chromium produces more orientation, 
a result in agreement with the bridging theory. 
Acids, bases, and salts, as well as all tanning 
materials, aside from chromium, in general appear 
to decrease the AS factor. This means that less 
orientation results under these conditions and 
implies reaction of collagen with these materials 
in a more or less random manner that may involve 
rupture of salt bonds, that is, decreased AH. The 
results obtained with aluminum are interpreted 
to mean a rupture of salt bonds and a correspond- 
ing disorientation sufficient in magnitude to offset 
the energy lost by rupture. The free energy in 
this case remains unchanged. Zirconium, on the 
other hand, tends to react in such a manner that 
few salt bonds are broken (normal A//). How- 
ever, a loss in orientation is found to occur since 
AS is smaller than for untreated tendon, as shown 
by the increased value of AFyo. In chrome tan- 
ning, a marked increase in A// occurs, but the 
increase in AS is correspondingly less, since AF x 
increases. A portion of the chrome, that is, 
possibly the whole amount > 1 percent of Cr,O; 
and a portion of the initial 1 percent, can be con- 
sidered to add nothing to the orientation. This 
means some loss of organization due to chrome 
tanning as well as to other types of tanning. 
These conclusions are borne out by observations 
that tanning reduces crystallinity, as shown by 
X-ray diffraction studies [38]. 


(d) Free Energy of Activation 


The free energy of activation is a measure of 
the position of equilibrium in the activation re- 
action. Since shrinkage will ensue rapidly, pro- 
vided sufficient numbers of activated molecules 
are formed, AF, will then be a measure of the 
shrinkage “‘temperature,” T7,. AF, will vary 
with temperatures by virtue of its dependence on 
the factor TAS. It is found that rapid shrinkage 
occurs at temperatures at which A¥’, ~ 25 keal/ 
mole, that is, rapid shrinkage may take place only 
if the activation equilibrium has been displaced 
sufficiently to provide required numbers of activ- 
ated molecules. 

The foregoing conclusion implies that treat- 
ments that reduce AS or increase AH will increase 
AF» and consequently the shrinkage “tempera- 
ture.’”’ Chromium tans principally by the latter 
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method, that is, increases AH up to 1 percent of 
Cr,O; and by the former in amounts exceeding 
this percentage. All other tannages tested appear 
to tan by reduction of AS. This reduction in 
AS may be pictured as resulting in an increased 
T, by steric hindrance phenomena. The amount 
of steric hindrance may be assumed to be propor- 
tional to the tannin content-——-percentage of 
Cr,0, in chrome tanning, pH in acid solutions, etc. 

Since 7, is determined by AF;, graphs were 
made using AF) as the ordinate and logarithm 
of the percentage chrome as the abscissa. It was 
found that a linear relationship resulted. If the 
logarithm of the concentration of solute, [NaCl], 
or the pH is used as the abscissa, straight lines are 
also obtained. Figure 13 shows typical results 
obtained for such graphs. 
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Ficure 13. Effect of chromium oxide content and sodium 
chloride on free energy. 


@,. Chrome tanned; ©, pickled. 


AF is related to the equilibrium constant of the 
activation process, A’, by the equation, 


AF;=—RT |\n kK’. (13) 


The results shown in figure 13 imply that the 
logarithm of the equilibrium constant is linearly 
related to the logarithm of the chromium content, 
concentration of solute, ete. 

Proceeding by speculation, assume K’=(C 
[(Cr.O;]", where C and n are constants. Then 
substituting in eq 3, the following equation is 
obtained: 

AF,=—RT In C [Cr,0,}". (14) 


In the case of rapid shrinkage it has been ascer- 
tained that AF> is nominally constant and ap- 
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proximately equal to 25 keal/mole at that tem- 
perature considered to be the shrinkage ‘“tem- 
perature.” Substituting this value for AF; and 
rearranging, 


25 < 10° og C ‘ 
log [CriOs]=— 9 393 Ratt ec. (15) 
The concentration of Cr,O; is not known, but it 
will be assumed proportional to the percentage 
of Cr,O;. On this basis if the logarithm of the 
percentage chrome is plotted against the recip- 
rocal of the absolute temperature of shrinkage, 
a straight line is expected. The data presented 
by Bowes, et al. [5] have been tested in this 
manner, and a straight-line results. Bowes 
showed that a straight line was obtained by 
using the temperature instead of the reciprocal 
of the absolute temperature but found it neces- 
sary to mention possible negative shrinkage 
temperatures, which might arise. The graph, 
according to eq 5, eliminates the difficulties of 
negative shrinkage temperatures and at the same 
times vields a satisfactory linear plot. Inasmuch 
as ¢ is proportional to 1/7 over the limited tem- 
perature range involved, the linearity using 
either ¢ or 1/7 is not unexpected. No con- 
clusions are possible concerning the slope of the 
line obtained, since the unknown factor n is 
involved. By using the experimental value 
given by Bowes, it is found that the value of n is 
approximately 1. 

No conclusions concerning the mechanism of 
shrinkage may be made from the present results 
at this time. It is believed, as a result of obser- 
vations during shrinkage, that the shrinkage 
reaction involves water either as a reactant or a 
catalyst. Experiments are in progress to ascer- 
tain the role of water in the shrinkage. 

The experimental results obtained in this 
investigation are in substantial qualitative and 
quantitative agreement with results of studies 
made on soluble proteins. Excellent summaries 
of these researches are given by Eyring and Stearn 
39] and Neurath, et al. [21]. 





The author is indebted to A. A. Maryott for 
valuable suggestions concerning the treatment of 
the data and the manuscript. 
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Kinetic Study of the Reaction of Carbon Adsorbents 


with Oxygen’ 


By William V. Loebenstein,? Leland F. Gleysteen,* and Victor R. Deitz 


A study was made of the kineties of oxidation of bone char and other carbon adsorbents. 


The measurements were limited to temperatures between 250° and 390° C. The results 
were expressed in terms of the fraction of carbon burned per hour at known values of oxygen 
concentration, temperature, and total rate of flow for activated carbons and bone chars. 
The rate of oxidation could be expressed either in terms of simple first-order or retarded 
first-order relationships. A nomograph was used to obtain a rapid estimate of the total 
fraction of carbon burned at any time for a given rate constant and degree of retardation. 
The relationship between the reaction-rate constant and the oxygen concentration was linear 
for a given adsorbent under the same conditions of temperature and flow. Effect of varia- 
tion in temperature on the rate of oxidation of a new bone char disclosed deviations from the 
Arrhenius equation even over a 100-deg interval. The influence of time, of oxygen concen- 
tration, and of temperature on decarbonization is discussed in relation to the conditions 
necessary in industrial operations for attaining the optimum uniformity in the carbon 


content of bone char. 








I. Introduction 


A study of the reaction of oxygen with carbon 
discloses many fundamental problems. ‘Two com- 
bustion products, CO, and CO, are always found 
simultaneously but in varying relative amounts, 
depending upon the experimental conditions 
employed. Heats of adsorption of oxygen on 
carbon are known to be very high initially but 
decrease rapidly as the amount of oxygen adsorbed 
increases. This phenomenon was observed as 
early as 1924 by Beebe and Taylor [1]* and by 
Garner and Blench [2]. The large range of acti- 
vation energies (20 to 70 keal/mole), which have 
been reported in the literature, may be ascribed 
to the complex reaction mechanism rather than 
to particular details in the experimental inves- 
tigations. 





! This investigation resulted from a joint research project undertaken by 
the United States Cane Sugar Refiners & Bone Char Manufacturers, a greater 
part of the refining industry of the British Commonwealth, and the National 
Bureau of Standards. 

? Research Associate at the National Bureau of Standards, representing 
the cooperating manufacturers. 

! Present address, Owens Illinois Glass Co., Toledo, Ohio. 

‘ Fieures in brackets indicate the literature references at the end of this 
Paper, 


Reaction of Carbon Adsorbents with Oxygen 


Aside from the theoretical interests in the 
reaction of oxygen with carbon adsorbents, there 
is an important practical application. The regen- 
eration of carbon adsorbents, which have been 
used for the purification of solutions, is accom- 
plished by heating in the presence of a limited 
amount of oxygen or oxygen-containing gases at 
elevated temperatures. The heating is intended 
to remove the adsorbed organic matter by oxida- 
tion to volatile compounds and restore as com- 
pletely as possible the adsorptive properties of the 
surface. 

After the exhaustion of the decolorizing power 
of carbon adsorbents in the refining of sugars, it 
has been demonstrated that it is technically pos- 
sible to regenerate almost all kinds of carbon 
adsorbents. However, there are large differences 
in the loss of weight realized during regeneration 
of adsorbents of different origin. This has an 
important bearing on the economic feasibility of 
the regeneration. The results given in this paper 
indicate a striking difference in the magnitude of 
the rates of oxidation of bone chars and of the 
activated carbons. 
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The raw materials from which carbon adsorbents 
are prepared vary considerably. The final prod- 
ucts may be viewed as carbonaceous residues of 
varying stability resulting from a destructive distil- 
lation of the raw material followed by a limited 
oxidation under variously controlled conditions. 
The physical and chemical properties of these 
residues vary and, consequently, the exact speci- 
fication of such residues is extremely difficult. 
The reproducibility of a particular product 
derived from a given raw material is not too easily 
maintained even when the requirements of a com- 
mercial application are considered. It is neces- 
sary, therefore, in order to meet the exact require- 
ments of fundamental research, that the desired 
experiments be made with aliquot parts of a 
well-mixed sample. The materials used in obtain- 
ing the results reported in this paper have been 
so prepared. Data on activated carbons and new 
bone char are reported in this paper. The results 
for spectroscopically pure graphite and spent bone 
char are also given for comparison purposes. 

In the experiments described in this paper the 
reaction is studied under conditions that maintain 
a large excess of oxygen in the gaseous phase com- 
pared to that consumed in the oxidation. By 
studying the reaction with respect to the carbon 
consumed it is possible to focus attention on the 
regeneration process. 


II. Experimental Procedure 


The sample for reaction with oxygen was sup- 
ported in a vertical tube with inlet for the gas 
mixtures at the top and the outlet at the bottom. 
This tube was 15 mm in diameter and 40 mm high 
and held samples varying in amounts from 0.5 
to 5.0 g. It was constructed of fused silica for 
the graphite experiments and of Pyrex glass for 
the carbon adsorbents. In order to insure uniform 
conduction of heat the reaction tube was fitted 
snugly into a cylindrical brass block 150 mm tall 
and75mm indiameter. This whole assembly was 
contained in a vertically mounted cylindrical 
furnace with a control thermocouple placed close 
to the windings of the furnace. A second thermo- 
couple was fitted into the center of the brass block 
through a hole situated so that the junction lay 
as close as possible to the outside wall of the 
reaction tube in the vicinity of the sample. 
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The reacting gas mixtures contained in gas 
cylinders consisted of known percentages of oxy- 
gen in nitrogen. The compositions had been 
determined by the Gas Chemistry Section of this 
Bureau. Purified nitrogen was passed through 
the reaction tube during the preliminary heating 
interval. The flow of nitrogen was discontinued 
by the employment of a three-way stopcock 
ahead of the reaction tube at the same instant 
that the oxygen mixture was introduced. Inter- 
vals of time were measured from this moment. 
The gas mixture was metered into the system at 
91 ml/min by means of a flowmeter and main- 
tained at this rate throughout the entire run. 

Preliminary tests showed that the water vapor 
contained in the exit gases, in comparison with 
the carbon dioxide and carbon monoxide, was so 
small that it could be neglected entirely. U-tubes 
of Ascarite-Anhydrone were used in the absorp- 
tion train to determine carbon dioxide. At the 
end of each predetermined time interval, one tube 
was removed from the train and a second tube of 
the same dimensions immediately inserted. The 
side arms of the U-tube were bent downward so 
as to fit into mercury-seal couplings and thus 
maintain airtight connections. The replacement 
of the tared U-tubes required but a few seconds. 
By alternating the twin U-tubes in this manner, 
it was possible to follow the course of the combus- 
tion by gravimetric determination without inter- 
rupting the conditions of flow and yet have a 
quantitative cumulative measurement. The re- 
maining gases, which contained residual oxygen 
and carbon monoxide along with nitrogen, were 
next passed through a hopealite catalyst, which 
converted the carbon monoxide quantitatively to 
carbon dioxide. This, in turn, was determined 
gravimetrically by means of a pair of U-tubes of 
Ascarite-Anhydrone in exactly the same manner 
as for the carbon dioxide. 

Each sample, prior to its controlled combustion, 
was heated to constant weight at 110° C in an 
atmosphere consisting either of helium or purified 
nitrogen [3]. 

The constancy of temperature during a com- 
bustion was maintained manually for periods of 
8 hours or less and by automatic means for longer 
periods of time. In the latter case, an indicating 
potentiometer controller was employed. Manu- 
ally, it was possible to hold a desired temperature 
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to within less than +1 deg C. A temperature 
control of about +3 deg C was readily attainable 
by automatic means. 

The total weight of carbon burned was com- 
puted from the yields of carbon dioxide and 
carbon monoxide. The rate of reaction was 
expressed as that fraction, f, of the initial carbon 
present in the sample which burned per hour. 
No attempt is made in this present Jreatment to 
interpret individually the amounts of carbon 
dioxide and carbon monoxide produced. 

The rate of flow during the experiments was 
chosen so that the oxygen concentration was not 
significantly decreased by its consumption. Furth- 
ermore, the residual oxygen concentration was 
sufficient for the quantitative oxidation of the 
carbon monoxide to carbon dioxide. This is a 
requirement in the analysis of the 
gaseous reaction products. Moreover, it was 
thought that by keeping the oxygen in excess, the 
interpretation of the reaction mechanism might 
be simplified. The oxygen concentration in the 
exit gases was verified later with a specially 
calibrated Pauling oxygen meter. 


necessary 


III. Kinetics of the Reactions 


The rate of combustion is expressed throughout 
this paper with respect to carbon. The control of 
the carbon content of bone char during the re- 
vivification process is of vital importance in char 
house operations of sugar refineries. Hence, the 
study of the combustion with respect to carbon 
lends itself most readily to practical interpretations. 

If the over-all combustion is of the first order 
with respect to carbon, the rate at which carbon 
burned should be proportional to the fraction of the 
carbon remaining, thus 


df 
5 =k(l—f) l 
dt I), ) 
where f is the fraction of initial carbon burned in t 
hours and & is the constant of proportionality 
(reaction rate constant). Upon integration, eq 1 
becomes 


2 303 logo = f kt, (2) 


with the constant of integration vanishing upon 
substitution of the initial conditions that f is zero 
when ¢ is zero. It will be seen from figures 1 and 


Reaction of Carbon Adsorbents with Oxygen 








007 
| 154% 


tt Ape 
iit Ms. aa 


aes 


Vn ae a 


YY 56 See 

















LOG if-1) 


> 
7 
‘ 








MO 


| 
002 - -— —_+—___ 4 
| 




















0.01 EE ee = 
° 1 | 
° 10 2.0 3.0 40 5.0 60 
TIME HOURS 
Ficgure 1, Reaction of activated carbon 6 with varying 


oxygen concentrations. 


The ordinate is logy 1/(1—f), where f is the fraction of carbon oxidized. 


The abscissa is the time in hours. The rate of flow is 91 ml/min, and the 


temperature is 390° C, 


2 that a linear plot of log 1/(1—/) against ¢ applies 
to graphite and the activated carbons, but not to 
any of the bone chars investigated. 

The combustion of bone chars will be shown to 
obey a retarded first-order relationship identical 
in form with that proposed by C. N. Hinshelwood 
[4] for the kinetics of the catalytic decomposition 
of nitrous oxide, in which it was assumed that the 
reaction was retarded by the oxygen produced. 
This does not necessarily imply that the retarda- 
tion in the combustion of bone chars can be attrib- 
uted to a similar mechanism. 

In the present treatment, the differential 
equation for combustion of bone chars may be 
written as follows: 


df k(i—f) ‘ 
dt 1+0f’ (3) 


where the constant 6 is a measure of the degree of 
retardation during a given combustion. This is 
also similar to an equation suggested by J. D. 
Lambert [5] for the oxidation of charcoal where the 
reactiop was investigated with respect to oxygen 
consumed. It can be seen from eq 3 that if no 
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Ficure 2. Oxidation of carbon adsorbents and graphite in 
6.1 percent of oxygen and a rate of flow cf 91 ml/min. 
The left ordinate scale is logy 1 (1—f), where f is the fraction of carbon oxi- 
dized, and the right ordinate scale is percentage oxidized (100/). A, Carbon 5 
(390° C); B, bone char 2 (280° C); C, earbon 6 (390° C); D, graphite (711° C); 
FE, carbon 4 (390°C); PF, graphite (600° C); G, earbon 4 (320° C 


retardation exists, 6 becomes zero and the expres- 
sion reverts to the special case of eq 1. It is also 
apparent, even with a finite value of 6 other than 
zero, that the limiting value of eq 3 as f approaches 
zero is the same as eq 1. Equation 3, upon in- 


tegration, vields 


1 l bh f k 
7 My +b ¢t1+0 


(4) 


which may be used by plotting 1/f In [1/(1—f)] 
against f/f to determine both k and 6. Equation 4 
may be rearranged, however, into the form 


l l l l » 
, In If k—b ; [ (in 1p) ‘| (5) 


which lends itself more readily to a direct evalua- 
tion of 6 and &, since these two quantities emerge 
immediately from the slope and intercept, re- 
spectively. 
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IV. Application of Simple First-Order Ki- 
netics to Carbon Adsorbents Other Than 
Bone Char 


Table 1 lists the values of f computed at the 
periods of time indicated during the course of five 
experiments with different samples of the same 
activated carbon 6. Each of the five experiments 
was conducted with a different percentage of 
oxygen in nitrogen, as listed in the first column. 
All five combustions were performed at the same 
temperature, 390° C, and at the same rate of flow 
of reactant gases, 91 ml/min. The data from 
table 1 are plotted in figure 1. The corresponding 
values of & computed from the slopes of these lines 
are recorded in the last column of table 1. Addi- 
tional curves are presented in figure 2 for activated 
carbons 4, 5, and 6. Some results are also in- 
cluded for samples of spectroscopically pure 


. ; . ° - ° 
graphite at 600° and 711° C. 
TaBLe 1. Combustion of an activated carbon 
[Reaction of activated carbon 6 with different oxygen-nitrogen mixtures at 
390° C and 91 ml/min flow showing the variations in the fraction of carbon 

burned, f, with time] 


Fraction burned at (hours) Rate 
= con- 
Oxy- ~— ° 
gen (cal- 
0.5 1.0 1.3 0 2.5 3.0 3.5 $0 45 5.0 - 
cu- 
lated) 
Per- 
cent 
2.2. 0.0129 0. 0246 9. 0367 0. 0496 0. 0630 0. OS63 0.118 0. 026 
34 .0208 .0868 . 0511 .0656 . 0706 wW7 0.120 . 134 Ol 
48 0197 .0347 20468 (0596 20749 0.0878 . 101 2118 0. 139 029 
6.1. .0231 .0418 .0583) . 0745 107 138 036 
15.1 0888 .0665 .0850 .120 . 146 . 060 


V. Application of Retarded First-Order 
Kinetics to Bone Char 


When values of log 1/(1—/f) were plotted against 
t for bone chars, the resulting graphs were never 
linear, but invariably showed a marked curvature, 
which was concave to the time axis. Such a curve 
for a typical bone char (char 2 at 280° C) is shown 
in figure 2 (filled-in circles) as compared to other 
activated carbons and graphite (open circles). 
Even in some of the latter cases a slight deviation 
from linearity is apparent, although the effect is 
not nearly as marked. The percentage of the 
original carbon oxidized at any time, ¢, for each of 
the experiments plotted on figure 2 can be esti- 
mated by reference to the vertical axis at the 
right. 
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Before attempting to apply the concepts of 
retarded first-order kinetics (eq 4 and 5) to the 
experimental data, it is necessary that zero time 
be known with more precision. One reason for 
this is that the time variable occurs in both co- 
ordinates of the proposed graphs for testing eq 4 
and 5. Although this complication did not exist 
in the simple first-order test as given in eq 2, a 
close examination of figures 1 and 2 will disclose 
that even in this case most of the lines do not pass 
through the origin. The displacement of these 
lines to the left of the origin in figure 1 may be 
due to the presence of chemisorbed oxygen in the 
samples. This remains despite the preliminary 
heating in pure nitrogen prior to the introduction 
of gaseous oxygen. There is also a tendency to 
displace the lines to the right of the origin. This 
is due to the dilution of the reaction mixture with 
the pure nitrogen and the time lag involved in 
sweeping out the dead space. The actual dis- 
placement observed suggests the dominating effect 
of chemisorbed oxygen. 

It is possible to avoid these errors by comput- 
ing a steady-state zero time from the original 
data for each experiment with bone char. The 
beginning of the steady-state condition was made 
to agree with a time at which the Ascarite-Anhy- 
drone tubes were changed early in the experiment. 

Figure 3 is a typical plot of data obtained for 
the controlled combustion of a sample of new 
bone char 2 according to eq 4 after corrections * 


' The weight of carbon, Aw, burned up to this time was subtracted from 
the original weight of the sample and also subtracted from the original weight 
of carbon, W, present in the sample. In the same manner, Sw was sub- 
tracted from the accumulative weights of carbon Jurned, w, subsequent to 
this time in order to get the corresponding accumulative steady-state values 
The steady-state fraction f, can then be expressed by the relationship 

w—Aw)/(n Aw). 
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Fiagure 3. Oxidation of new bene char 2 in 6.1 percent of 
oxygen and at 91 ml/min. 


The abscissa is f/t and the ordinate 1/t log [1/(1—/)], where / is the fraction 
of carbon oxidized, and ¢ is the time in hours 
had been made for steady-state conditions. A 
similar illustration for a spent bone char (char 
15) is given in figure 4. 
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Ficure 4. Oxidation of spent bone char 15 in 6.1 percent 
of oxygen and at 91 ml/min, 


The plot is similar to that in figure 3. 


TasBie 2. Combustion of a bone char 


[Reaction of new bone char 2 with different oxygen-nitrogen mixtures at 280° C and 91 ml/min flow showing the variations in the fraction of carbon burned 
with time. The data are referred to steady-state conditions for each combustion.} 





Fraction burned at (hours)— Retar- Rate 


Oxygen meen pearomenee nae - aces meeoneen t 


0.5 1.0 1.5 2.0 2.5 3.0 3 
Percent 
2.2 ‘aia “ .. 0.0672 0.0829 0.0985 0.124 
34 .. 0.0280 0.0527 .0720——. 0906 0. 123 
48 .| .03815 .0559 | .0768 Oo78 £116 . 130 145 
5.6 ; ..| .0865 0661 =. 0905 oe 151 . 167 
7.0 ..| .0449 OS17—. 110 w--| «180 . 181 
9.1 aie .-| .0423 | .0756 .1083 .127 weal . 169 
15.1 0551 . 0957 . 160 eee - 212 


dation con- 
factor, stant, 


4.0 4.5 50 5.5 6.0 6.5 7.0 b 
it 0. 149 0.172 0. 191 . R 7.3, 0.057 
0.150 aan 0.173 0. 197 — 0. 218 , 008 
| .161 . . 187 ‘ . 231 6.9 . 069 
. 184 . 2038 - 332 |... . oe lietaucs 6.7 O83 
. 223 . . . 260 . 295 “ . 326 3.2 OST 
6.5 . 097 
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Table 2 lists the steady-state values of f for 
new bone char 2 at 280° C and 91 ml/min of flow 
at various oxygen concentrations of the reactant 
gas mixtures. In the last two columns of the 
table are the calculated values of 6 and k com- 
puted graphically when the data were plotted 
according to eq 5. 


VI. Nomograph 


If it is desired to estimate the fraction of carbon 


f that is burned at any time, ¢, assuming both 6 


and k& to be known, the solution of eq 4 or 5 
becomes rather laborious. By the use of a nomo- 
graph, however, this may be accomplished in a 
matter of seconds. Such a nomograph has been 
constructed, figure 5, which may be utilized for 
any first-order combustion described in this paper 
regardless of the temperature, oxygen concentra- 
tion, rate of flow, or degree of retardation. A 
straight edge is laid across the sheet connecting 
the particular value of 6 located on the vertical 
axis at the right with a value of kt located on the 
seale at the left. Where this straight edge crosses 
the curved line the corresponding value of f is 
obtained. For example, consider first the com- 
bustion of activated carbon 6 in 2.2 percent of 
oxygen at 390° C and 91 ml/min flow, as shown 
in the top row of table 1. In this case k=0.026 
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Ficure 5. Nomograph fer the rapid solution ef equation 
4 or 5. 


The left vertical scale is expressed in units of k t, where & is the reaction rate 
constant (initial fraction of carbon burned per hour), and ¢ is the time in 
hours. The right vertical scale is given in terms of 6, the degree of retardation. 
The intermediate scale is in units of f, the fraction of carbon burned. 
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and 6=0. Since it is desired to estimate the frac- 
tion of carbon burned after 4.5 hr, kt=4.5 < 0.026= 
0.117. The nomograph gives directly an f-value 
of 0.11. This compares quite well with the 4.5-hr 
value of f (table 1), which is 0.113. 

As a second example, consider the combustion 
of new bone char 2 in 3.4 percent of oxygen at 
280° C and 91 ml/min flow. These data are found 
in the second row of table 2. In this case 6 and k 
were calculated to be 7.1 and 0.063, respectively. 
It may be of interest to estimate the fraction of 
carbon which would be expected to burn in 1 hr 
under this particular set of conditions. For ¢t=1, 
kt=0.063, and the nomograph gives an f-value of 
0.052. The corresponding value of f in table 2 is 
0.0527. It is thus seen that the nomograph has an 
exceedingly wide range of applicability. It is 
possible to extend either or both of the vertical 
axes in order to accommodate still larger values 
of b and k, since each of these two axes is marked 
off in a linear scale. 


VII. Influence of Oxygen Concentration 
and Temperature on the Oxidation of 
Some Carbon Adsorbents 


When the specific reaction rate constant k was 
plotted against the corresponding oxygen percent- 
age in the reaction mixture under the same con- 
ditions of temperature and rate of flow, the 
resulting points were again best correlated by a 
linear relationship. Curve A in figure 6 shows the 
dependency of k upon oxygen percentage in the 
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Ficure 6. Rate constant k as a function of oxygen concen- 
tration, 


The line A is for the oxidation of xctivated carbon 6 at 91 ml/min flow and 
390° C (see table 1). Line B is for the oxidation of new bone char 2 at 91 
ml/min flow and 280° C (see table 2). The ordinate is the rate constant k, 
and the abscissa is the oxygen percentage. 
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TABLE 3. 


Combustion of a new bone char in 6.1 percent of oxygen in nitrogen at a total velocity of 91 ml/min. The data are 


referred to steady-state conditions for each combustion, 


Fraction burned at (hours)— 


Temperature, ° C a Pe we ————— ae 
0.5 1.0 15 2.0 2.5 3.0 
250_. 0.0106 0.0199 0.0274 0. 0416 
280 ° O47 0580 0. 0727 UST4 0. 100 
3a om 0843, 44 v0 - 223 . 260 
325 OR60 me Bs | 200 277 310 


reaction mixture for the combustion of activated 
carbon 6. The data used for this graph were 
obtained from the first and last columns of table 1. 
Curve B was constructed in the same manner by 
using the data for new bone char 2 as given in 
table 2. The fact that neither of these lines 
appears to intersect the origin of coordinates is 
quite striking. The magnitude of this deviation 
is greater than can apparently be accounted for on 
the basis of experimental error. 

The reaction rate constant for the combustion 
of a new bone char was determined at four tem- 
peratures within the range of 250° to 325°C. In 
each case the reactant gas consisted of 6.1 percent 
of oxygen in nitrogen flowing at a total velocity of 
91 mlmin. It was not practical under the 
conditions employed to exceed this range of 
temperatures. When this was attempted, the 
combustion proceeded either too rapidly to obtain 
reliable measurements or too slowly to insure 
constancy of conditions over the period of time 
required to obtain measurable differences in 
weight. The steady-state values of the fraction 
of carbon burned at each interval of time are 
recorded in table 3. Corresponding values of k 
and 6 computed from eq 5 are also tabulated for 
each of these temperatures. 

The dependence of k on temperature is usually 
expressed by the well-known Arrhenius equation 


In k-=—E/RT +1n ko, (6) 


where F is the activation energy (or free energy of 
activation), 2? the gas constant, and fp is a con- 
stant. Figure 7 shows that the logarithm of & 
when plotted against the reciprocal of the absolute 
temperature, gave a smooth curve, which was 
slightly convex to the 1/7 axis. An over-all 
activation energy (or free energy) computed from 
the slope of this curve would seem to vary between 
15 and 35 keal under the conditions of these 





Reaction of Carbon Adsorbents with Oxygen 


81519749 4 





















































SE ee ee Rate con- 
| 5 factor, b stant, & 
3.5 4.0 4.5 5.0 5.5 6.0 6.5 ° 
0. 0564 0. 0664 0. 0772 0. OS61 13.5 0. 0224 
0.123 0.144 0. 164 7.9 . 0495 
6.0 . 221 
. 364 5.3 231 
-04 
\ 
-06 \ 
No 
-08 
\ 
\ 
-1o \ 
“ \ 
= \ 
-12 x arene 
\ 
° 
‘ 
\. 
-14 Vv 
\ 
\ 
-\6 q 
N\ 
‘ 
-1.8 
1.5 1.6 7 18 9 20 
1000/T°K 


Ficgure 7. A plot of In k against 1/T for new bone char 2, 


where T is the absolute temperature. 


experiments. However, the deviation from linear- 
ity definitely points to the need for further inves- 
tigation before any reliable interpretations are to 
be made concerning this apparent anomaly. Also, 
the physical significance of the parameter 6 along 
with its dependency upon the other variables 
present is beyond the scope of this paper. How- 
ever, there is no doubt that the reaction rate does 
depend more or less exponentially upon 1/7. 
Other combustion experiments in which only the 
rate of flow of reactant gas mixture was varied by a 
twofold factor or more seemed to have very little 
effect upon the reaction rate constant. Further- 
more, the weight of the sample was varied under 
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otherwise identical conditions from one combustion 
experiment to another, yet the total fraction of 
carbon burned at each of the similar time intervals 
remained constant. These results, together with 
the linear effect of oxygen concentration upon the 
combustion rate, form a fairly complete picture as 
to the various effects of each factor. 


VIII. Decarbonization of Bone Char as 
Influenced by Time, Oxygen Concen- 
tration, and Temperature 


Uniform decarbonization of service bone char is 
an important factor in industrial installations in 
sugar refining [6]. In practice, methods to control 
decarbonization have entailed deliberate varia- 
tions in time, oxygen concentration, or temperature 
that were involved in the burning process. The 
relative effects of changes in each of these variables 
are illustrated in figure 8. It may be seen in figure 
8, a, that the rate of oxidation diminishes. Figure 
8, b, on the other hand, illustrates that the changes 
in the fraction of carbon burned are always directly 
proportional to increases in the oxygen concen- 
tration. In contrast, the effect of temperature as 
shown in figure 8, ¢, is considerable. The fraction 
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Ficure 8. Influence of time, oaygen concentration, and 
tem pe rature on decarbonization of a bone char. 


The ordinate scale common to the three curves is given in terms of f/. The 
abscissa of curve A is in hours (6.1% of O», 280° C, 91 ml/min); that of curve B 
is in percentage of O » (1.5 hr, 280° C, 91 ml/min); that of curve C is the tem- 
perature °C (1.0 hr, 6.19% of Og, 91 mi/min 


of carbon burned becomes progressively more and 
more sensitive to changes in temperature as the 
latter increases. 

It is thus’seen that temperature is by far the 


most important single factor in influencing 
the extent of decarbonization. The temperature 
range investigated in the laboratory is appreciably 
lower than the peak temperatures (400° to 500° C) 
realized in industrial operations. At these high 
temperatures it was impossible to measure the 
exceedingly high rates of combustion with the ap- 
paratus described. An extrapolation from the 
data in the temperature range 250° to 325° to 
400° C or higher may be seen not to be feasible 
from either figure 7 or 8, c. 

Some quantitative comparisons regarding the 
conditions affecting the variation in the fraction 
of carbon burned are evident in the contents of 
table 4. There are two ranges of f-values, which 
have immediate application to refinery practice. 
In a single cycle of regeneration with the conven- 
tional kiln, less carbon is to be removed than in 
a single passage through a decarbonizer. In the 
first case the order of magnitude of the change in 


f is 0.02 to 0.03; in the second case it may be as 


high as 0.13 to 0.14. In order to realize the 
desired carbon removal by adjustment in temper- 
ature alone (1.0 hr, 6.1% of oxygen, and 91 
ml/min), it is necessary to maintain a temperature 
of 260° and 320° C, respectively. While a fluetu- 
ation of + 10.5 deg can be permitted at the low 
temperature for the desired carbon removal, a no 
greater variation than +1.0 deg is allowed at 
820° C. Correspondingly, to obtain the desired 
decarbonization by adjustment of oxygen concen- 
tration alone (1.5 hr, 280° C, and 91 ml/min), it is 
necessary to hold the concentration within the 
range of +0.85 percent in either case. Finally, to 
obtain the decarbonization by adjustment of the 
time factor alone (6.1% of oxygen, 280° C, and 
91 ml/min), it is necessary to heat the char for 
0.6 and 4.5 hr, respectively. A fluctuation of only 
+7 min is permitted at the lower time, while as 
much as £15 min may be tolerated at the higher 
time, 

These results definitely point to the following 
criteria for uniform decarbonizatior: (1) long re- 
action times, (2) low temperatures, and (3) what- 
ever concentrations of oxygen are required to give 
the desired decarbonization. A rough visual in- 
dication of such uniformity is evidenced by the 
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absence of the familiar black and white appear- 
ance of decarbonized char which is observed in 
most refinery operations. A more precise esti- 
mate of uniformity requires numerous laboratory 
analyses. There is an additional advantage in 
low-temperature operation based on considera- 
tions of the liberation of the heat of combustion. 
Temperature control is dependent upon an ade- 
quate dissipation of this heat which, in turn, de- 


TABLE 4 Conditions affecting fluctuations in } 


Perm issable deviation of any one 
variable 
(Assumed 
nge of f Oxygen 
T inme concen em perature 
tration 
D ce f I 
0.02 to 0.08 +0). 125 +105 418.0 
13 t 14 +. 2%) +f. 85 +1.0 +18 
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pends upon the heat transfer away from the char, 
Because of the poor heat conductivity of char, 
the necessary heat transfer is obtained only at the 
low temperatures, where the heat of combustion 
per unit time can be kept at a minimum. 
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results after 1 year of exposure. 


City, N. J. 


specimens Was unimpaired after 7 years. 


I. Introduction 


During the past 2 decades, the increased use of 
porcelain enamel as an architectural material, 
including its recent application in the low-cost 
housing field [1, 2],' has aecentuated the need for 
adequate data on the weather resistance of various 
types of enamel and also for laboratory tests to 
indicate weather resistance. 

An investigation designed to obtain such data 
was begun by the Enameled Metals Section of the 
National Bureau of Standards in 1939 and was 
planned with the assistance of an advisory com- 
mittee from the industry. Through the coopera- 
tion of 16 manufacturers, 864 1-ft-sq. panels and 
an equal number of 4- by 6-in. laboratory speci- 
mens were prepared. Most of the enamels 
furnished were regular commercial products, but 
they were not, in all cases, enamels that had been 
proved suitable for architectural purposes. On 
the contrary some were not expected to have 
good resistance to weathering. The exposure 
sites selected were Washington, D. C.; St. Louis, 
Mo.; Lakeland, Fla.; and Atlantie City, N. J. 

The present paper, which describes the test 
panels after 7 yr of exposure, is the second prog- 
ress report of this investigation. The first 
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at the National Bureau of Standards in 1939. 
and a like number of 4- by 6-inch laboratory specimens. 


ance of enamels and their resistance to weathering. 


Weather Resistance of Porcelain Enamels Exposed for 
Seven Years 


By William N. Harrison and Dwight G. Moore 


A study of the weather resistance of porcelain-enameled architectural panels was begun 
The study involves 864 1-foot-square panels 


A previous report (1942) gave the 


The present paper describes the condition of the panels after 


7 years of weathering at Washington, D. C., St. Louis, Mo., Lakeland, Fla., and Atlantic 


The results of the 7-year inspection showed a good correlation between the acid resist- 


Where initial coverage was complete 


and no mechanical damage had occurred, protection of the metal by the enamel for all 


Fading of colored enamels occurred only on those panels with very poor acid resistance. 


report, describing the results obtained for 1 year 
of exposure, was published in 1942 [3]. The lapse 
of 6 years between the first and second progress 
reports was due to the impracticability of carrying 
out a regular inspection program during the war 
years, 


II. Conditions Prevailing at Locations of 
ure 


Table 1 lists the exposure locations, and table 2 
gives pertinent data on weather conditions during 
the 7-vr period of exposure at each site. At all 
four locations the racks face south, the panels 
being exposed at 45° from the horizontal. 


Figures in brackets indicate the literature references at the end of this 


paper. 


TABLE 1. Exposure test locations 


Exposure conditions 


City “xposure s 
it Exposure site represented 


Washington, D.C Roof, Industrial Building, 
National Bureau of 
Standards. 

Roof, Union Electric Co. | Temperate, industrial. 
warehouse. 
Ground, Municipal Air- | Semitropical, residential. 
| port. | 

Ground, U. 8. Coast Temperate, “salt air.” 

Guard Station. | 


Temperate, residential, 


St. Louis, Mo_. 
Lakeland, Fla 


Atlantic City, N. J 














TasBLe 2 General weather data for the first years of 
exposure from S. Weather Bureau records 
Annual Annual Aver- 
( Exposure pe i rain sun a 
. fall shine tem per 
- ture 
h I 
Was! Ly. « Db 1939 t ugh N Ls 2 507 57.0 
146 
St. Lou M Ap rt wh Ma ‘ 270 
I 
Lakeland, Fla July 1940 tl wh March i 2.4 72 
17 
Atlantic ¢ NJ Au HM l April 10.0 2 TAl 7 
“4 
Avera mputed from data for actual period of exposure 
*Taken from Tampa, Fla., record rotal sunshine for Lakeland not 


III. General Description of Panels 


The 14 types of enamel included in the study are 
indicated in table 3. For a more complete deserip- 
tion of the specimens and their fabrication, refer- 


ence is made to the first report [3]. 


The supporting racks were constructed of angle 





the general condition of both panels and racks. 
The following remarks summarize their condition 
at each location. 


1. Panels at Washington, D. C. 


All panels and racks were in comparatively good 
condition. Periodic painting of the steel racks 
had been possible during the war years, and no 
serious corrosion had taken place. The panels 
were substantially free from dirt and grime and 
the natural washing by rain was apparently all 
that was required at this location to keep surface 
deposits from accumulating. 

On one panel (V-71) the surface had been ac- 
cidentally fractured in April of 1941, resulting in 
a local exposure of the metal. It was estimated 
that during the interventing 6 yr, corrosion had 
progressed into the metal to a depth of 0.008 in., 
but there was apparently no appreciable penetra- 
tion of the corrosion under the adjacent enamel. 
The size of the fracture had not increased from 
its original 's-in. diameter. 





Ficure 1. 


Partial view of exposure-test installation on roof of Union Electric Co. warehouse in 


St. Louis, Mo., near the 


railway terminal. 


Photograph taken 


iron and, after priming, were painted with alumi- 
num paint. Seven racks were required for each 
location. Figure 1 shows a part of the installation 
at St. Louis. 

Figure 2 is a view of the reverse side of one of 
the panels. The two clips at one side (C) and 
the flange extension at the opposite side (fF) both 
fitted into galvanized channels on the exposure 
rack. The spaces between the panels were not 
caulked but were left open to facilitate removal of 
the panels for inspection. 


At the 7-yr inspection, notes were taken as to 
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n 1941 before heavy deposit of combustion products had formed on panel surfaces (see page 48) 





Figure 2. View of the 
attachment of clips (¢ 


reverse side of a panel showing 
and lower flange extension (F) used 


for fastening panels to racks. 
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resistance lests-on same compositions 


| Average 
Fabri- | 45° initial 























Specimen identification | 
cator specular Pasian, ; 
gloss | Wasting | st. Louis?) Atantic | Lakeland | 
WHITE, GLOSSY, ACID-RESISTANT ENAMEL 
Percent 
A-ltos a 5.70 S14 9.6 4.5 77.5 
A-ll to IS b 5.74 83.2 63.6 87.2 
A-21 to 2s a 6.04 79.8 76.5 S57 
A-31 to 38 b 6.01 74.8 wo.4 05.4 
A-41 to 48 a 5. 85 4.1 SL5 77.8 
4-51 to 5S b 5.81 77.5 66.6 6 
A-fil to 68 a 6.20 74.6 74.8 83.7 
A-71 to 78 b 6.01 70.5 92.9 61.2 WO. 7 
Avg 78. 2 (8.7 81.2 
WHITE, GLOSSY, NONACID-RESISTANT ENAMEL 
B-1 to8 b 5. 26 2.9 “4 37.7 
B-21 to 2s b 5.41 34.5 36.4 2.4 
B-41 to 48 b 5.16 26.5 48.2 37.7 2.4 
B-451 to 68 b 5. 32 34.8 4.4 37.7 
Ave 31.2 40.5 32.0 
- — - : - 
WHITE, SEMIMAT, ACID-RESISTANT ENAMEL 
C-11 to 18 c 4.06 58 > 100.0 2 45.0 
C-31 to 38 c 5.24 61.5 42.0 (4.58 
C-51 to 58 c 5.16 76.8 581 80.5 
C-71 to 78 c 5.32 73.1 2 70.3 
Avg 65.8 52.6 65.1 
WHITE, SEMIMAT, NONACID-RESISTANT ENAMEL 

D-l tos c (5) 
D-21 to 28 € 5. 69 43.2 46.0 41.5 
D-41 to 45 c 5.51 38.0 31.9 36.6 

D-5l to 68 c (4) 
Avg 40.6 38.9 38.2 

BUFF, GLOSSY, ACID-RESISTANT ENAMEL 
d 5.09 2 12.0 57.8 
d 5.43 4.2 51.2 7¥.6 
d 5.44 74.0 70.4 SOLS 75.6 
d 5.35 77.7 71.0 fl.l 
71.0 61.2 66.0 
BUFF, GLOSSY, NONACID-RESISTANT ENAMEL 

F-1 to8 d 5.15 43.0 33.5 21.9 
F-11 to 18 e 4.87 61.2 4.0 24.7 
F-21 to 28 : d 5.02 52.9 76.0 1 35.7 
F-31 to 38 e 5.56 45.5 4 31.8 
F+4lto 4 - d 4.31 37.3 406 40.1 
F-51 to 58 . e 5. 66 33.7 3u.7 4.6 
F-61 to 68 d 44 22. 5 42.1 36.8 
F-71 to 78 e 5. 26 2.0 46.7 25.3 
es 40.9 43.6 32.6 






See footnotes at end of table. _ 
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TABLE 3. Percentage of initial specular gloss retained by 14 enamels for ? years at four exposure locations and results of acid- 
resistance lests on same compesitions— Continued 


Average percentage of initial specular gloss retained at 


A verage ‘ Acid re- 

., ® , j sistance 

Specimen identification Fabri 45° initial } istance 
cator specular Washine- Atlantic Average PEI 

gloss ‘ton 8° | St. Louis? | ¢ City Lakeland for three test ¢ 


locations 


BUFF, SEMIMAT, ACID-RESISTANT ENAMEL 





| 
| 
j 
( 
| Percent 
| H-i tos e 3.06 0.7 62.5 55. 2 50. 5 4 
! H-11 to 18 f 4.81 57.6 70.6 57.0 61.7 A 
| H-21 to 28 e 5. 65 70.7 63.1 “6.9 4 
1 H-31 to 38 f 5.45 71.2 40.1 70.5 60.6 \ 
' H-41 to 48 : e 4.74 75.7 > 100.0 80.6 81.9 79.4 AA 
| H-S1 to 58 f 4.85 72.5 60.6 74.5 42.2 AA 
| H-41 to 6s e 5. 51 74.4 100.0 uM. 6 SO 73.0 AA 
H-71 to 78 f 5. 54 76.4 m2 S14 72.7 AA 
l] 
Ave 9 62.9 71.5 i. 3 
BUFF, SEMIMAT, NONACID-RESISTANT ENAMEL 
K-l to8 f (i D 
K-11 to 18 g (3) D 
K-21 to 2 { 5.41 49.0 75.3 47.3 IAS 45.0 D 
K-31 to 38 e g 5.37 36.4 48.2 27.6 37.4 D 
K-41 to 48 f 5. 20 63.5 73.8 49.5 46.7 53.2 D 
K-51 to 58__.. , ‘ ‘ ru 5.35 36.5 2.4 29.4 36.1 D 
K-61 to 68 : f (*) D 
K-71 to 78 iy (5) D 
Avg ‘ 46.4 46.9 35.6 42.9 
RED, GLOSSY, ACID-RESISTANT ENAMEL 
L-1 to8 £ 5. 5s 76.4 97.1 77.7 $9. 2 “4.4 AA 
L-11 to 18 h 5. 61 74.7 SS. 6 96. 5 &6. ¢ AA 
L-21 to 28 £ 5. 42 78.1 57.4 86.7 74.1 B 
L-31 to 38 h 5.64 79.7 49.6 92.8 74.0 \ 
L-41 to 48 g£ 5. 58 75.7 78. 7 94.1 82.8 AA 
L-Si to [8 h 5. 23 67.6 79.7 S66 78.0 A 
L-4il to && ¥ £.30 73.2 76.3 91.0 0). 2 A 
L-71 to 78 h 4.18 74.0 82.7 78.1 78.3 AA 
Avg 74.9 73.8 00.6 7.8 
RED, GLOSSY, NONACID-RESISTANT ENAMEL 
N-1 to8 h 5.17 | 70.0 | } 62.3 40.2 | 7.5 
N-I1 to 18 » 5. 08 0.9 | 55.0 | 39.9 51.9 Cc 
N-21 to 28 h 5.12 74.3 48.4 | 46.1 56.3 D 
N-31 to 38 k 4.55 70.4 59.2 53. 6 61.1 Db 
N-41 to 48 h £.17 65.0 62.4 35.0 54.1 Cc 
N-51 to 58 k 4. 68 61.4 52.8 52.8 55.7 Cc i 
N-6l to 68 h 4.71 43.0 7.4 32.9 44.4 D 
N-71 to 78 k 5.44 36.4 38. 4 52.9 44.5 M44 D 
Avg.. 60.2 %.3 43.2 54.4 


See footnotes at end of table. 
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id- TaBLE 3. Percentage of initial specular gloss re tained by 14 enamels for ? years at four exposure lecations and results of acid- 
resistance tests on same compositions Continued 


Average percentage of initial specular gloss retained at 
— Acid re- 


Average 


= 8g : ; 
ce —— atten Fabri- | 45° initial ! sistance 
Specimen identification cator specular | washing- astenste Average PEI 
‘ gloss - ® | st. Louis?) “G: Lakeland for three test ¢ 
ton City } 
jocations 
RED, SEMIMAT, ACID-RESISTANT ENAMEL 
Percent 
P-l tos k 5.54 70.3 77.1 S1.0 7H. 1 AA 
P-11 to 18 1 (5) . ° . ° AA 
P-21 to 28 k 1.41 62.6 “0.9 49.1 65.4 9.0 c 
P-31 to 3s ] 3.99 4.4 58.7 42.9 55.3 Cc 
P-41 to 48 . 4.95 i. 2 6.7 75.4 7u.4 A 
P-5S1 to 58 l 3.34 73.8 73.4 4.1 77.1 A 
P-4l to 68 k 4.97 76.2 6.4 92.9 78.5 AA 
P-71 to 78 ] 5. 54 S1.0 81.3 97.2 86.5 AA 
Avg 71.1 67.5 77.0 71.1 
RED, SEMIMAT, NONACID-RESISTANT ENAMEL 
S-l to8 ] (5) D 
S-ll to 1s m ©) D 
S-21 to 2 | 4.51 65.8 51.1 4.1 55.0 B 
8-31 to 38 m 4.51 57.6 47.2 45.1 0 ( 
S-41 to 48 7 l 3.50 58.0 49.8 32.0 46.6 4 
S-51 to 5S m 3.95 4.1 75.8 44.4 is. 1 \ 
S-H1 to 68 l 4.50) 35.2 47.9 30.4 37.8 ( 
S-71 to 78 m 3. 57 70.3 53.3 40.2 4.6 ( 
Avg 61.8 54.2 40.0 52.0 
BLACK, GLOSSY, ACID-RESISTANT ENAMEL 
lr-l to8 m 7.30 56.1 47.8 56.3 53.4 AA 
r-1l to 18 n 6.600 (4.2 90. 7 62. 2 79.6 i. 7 AA 
T-21 to 28 m 5. 38 80.6 75.7 96.3 84.2 AA 
I-31 to 38 n §. 55 70.0 7u.4 03.9 81.1 AA 
r-41 to 458 m 6.09 70.7 65.6 aA. 4 74.9 AA 
T-51 to 58 n 5. 90 ‘9 6.9 &2.5 72.8 A 
r-61 to 68 m 6. 55 59.8 53.9 61.8 58.5 AA 
1-71 to 78 n 6. 56 54.9 77.6 63.1 6.3 70.5 AA 
Avg 65.6 64.3 78.4 w4 
BLACK, GLOSSY, NON ACID-RESISTANT ENAMEL 
V-ltos n 5. 52 63.6 50.8 36.8 3.4 c 
V-ll to 1s a 5. 62 57.0 SiS 43.9 .6 Cc 
V-21 to 28 n 5.76 SS 61.8 41.2 53.9 Cc 
V-31 to 38 a 5.0 40.1 40.5 43.3 41.3 
V-41 to 48 n 5.30 0.7 77.4 67.2 43.8 me Cc 
V-51 to 58 a 5.49 58.0 65.5 0.4 M6 Cc 
V-4l to 6s : n 4. 67 75.9 a9 67.0 67.3 Cc 
V-71 to 78 a 5.45 73.9 70.1 65.8 6.9 Cc 
Avg 0.9 1.4 47.8 56.0 
1 Groups of 8 panels exposed, 2 at each location. Ninth panel kept in 2 Values are average of 2 panels with 2 readings on each panel. 
storage. «Test made on storage panel according to Porcelain Enamel Institute 
2 Cleaned by nonstandard procedure of prolonged scrubbing with commer- standard acid-resistant spot test for flatware issued April, 1940. 
cial scouring powder. $s Full-mat enamels, initial gloss too low for measurement. 


ch Weather Resistance of Porcelain Enamels 47 














2. Panels at St. Louis, Mo. 


The 192 panels at the St. Louis location (on a 
roof near a railway terminal) were found to have 
a dark gray deposit over the entire surface, which 
was sufficiently heavy to make all panels have 
approximately the same appearance, regardless of 
their original color or gloss. An examination 
showed this deposit to consist of fly ash, soot, and 
a tar-like substance. A small sample seraped 
from the surface of one panel of high acid resist- 
ance showed an ignition loss of 31 percent, indi- 
cating that the greater part of the deposit (69 
pereent) consisted of the noncombustible fly ash. 
The deposit originated from soft-coal combustion 
products that were present in high concentration 
in the exposure area. 

The standardized cleaning procedure of wash- 
ing with I-perecent trisodium-phosphate solution 
did not remove this deposit, and vigorous scrub- 
bing with a scouring powder was the only method 
found to give satisfactory cleaning. Unfortu- 
nately, this scrubbing had a polishing action that 
affected the gloss measurements and vitiated their 
reliability as a criterion of the degree of weathering. 

The aluminum paint on the steel supporting 
racks had deteriorated badly at the St. Louis 
location, but corrosion of the steel was not as vet 
pronounced. The previously mentioned deposit 
was also present on the steel racks and may have 
had a protective influence in decreasing the 
corrosion rate. 

3. Panels at Lakeland, Fla. 

Tiny patches of an organic growth were attached 
to all panels at the Lakeland exposure site. The 
patches were largest and most numerous on the 
white enamels, but the black panels also showed 
some evidence of the growth. The patches 
adhered tightly to the enamel surface but could 
be removed by vigorous washing with a wet rag. 
When dry, the patches were gray but when thor- 
oughly wetted they became green or, in some 
cases, brown. 

An examination of these patches was made at 
the Florida Agricultural Experiment Station at 
Lake Wales. Under the microscope the patches 
appeared as clusters of oval-shaped fungus spores 
with occasional algae growth present, especially 
on the white panels. 

On many of the nonacid resistant enamels, 
there were areas having much higher gloss than the 
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remainder of the surface. Examination showed 
these areas to be substantially free of the growth, 
indicating that these organisms produced an acid 
reaction that accelerated the deterioration of the 
enamel surface. None of the panels with acid 
resistance of class B or better showed this localized 
etching 

Jones [4], in explaining the behavior of optical 
glass elements used in the tropics, states that fungi 
produce organic acids that etch glasses subject to 
acid attack. 

Corrosion of the steel racks was only moderate 
at the Lakeland location. The aluminum paint 
had deteriorated but rusting was not serious. 


4. Panels at Atlantic City, N. J. 


At Atlantic City, damage to the installation 
was considerable and was caused by a combination 
of factors including the following: (1) The 1944 
hurricane that struck the site with full intensity, 
(2) moving of the racks to make room for war- 
time expansion at the Coast Guard Station, (3) the 
location in 1943 of a Navy athletic field immedi- 
ately adjacent to the exposure site, and (4) the 
corrosive action of salt air on the steel. Nine of the 
panels were missing, and seven others had become 
detached from the mounting and were on the 
ground. Many of the panels had been damaged 
by impact, probably caused by stones thrown 
from the adjacent athletic field. Nevertheless, 
surface measurements were possible on practically 
all of the panels. 

Observations at this site indicated that both the 
paint and galvanizing failed rapidly under the 
salt-air conditions. Corrosion had seriously dam- 
aged the painted steel racks and the galvanized 
channels were, in many cases, so rusted that 
replacement was necessary. 

The top attachment clips (C in fig. 2) had be- 
come detached on some of the panels. These clips 
had been spot-welded to the panels prior to 
enameling, and, in some cases, this joint was 
protected only by a thin ground coat application. 
On many of the clips this thin coat was insufficient 
to give complete coverage where the clip joined the 
panel, and the severely corrosive conditions had 
caused failure. 

This severe corrosion of the metal at Atlantic 
City was also responsible for a somewhat unex- 
pected type of damage to the enamel surface 
resulting from poor enamel coverage of the backs 
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if panels. On those panels that had been covered 
m the reverse side with only a ground coat and 
hen fired in a box furnace while resting on alloy 
points, poor coverage of the steel occurred at the 
contact points. Corrosion had started at these 
areas and had, in some cases, progressed com- 
pletely through the 16-gage steel to the under sur- 
face of the enamel on the face side. Fracture of the 
enamel on the face immediately above these areas 
occurred as the corrosion approached the enamel- 
metal interface. The resulting fractures resembled 
very large fish scales and were believed to have 
been caused by hydrogen generated during the 
corrosion process. That hydrogen diffusing 
through steel may literally explode the enamel 
from the opposite surface had been demonstrated 
previously by Zapffe and Sims [5]. The fact that 
a number of “‘fish-scale’”’ type fractures were noted 
on the faces of panels over localized corrosion 
areas while some metal still remained beneath 
indicates that the fractures could not have been 
caused by the foreing off of the enamel by corrosion 
products, nor could they be caused by buckling of 
the enamel due to the release of compressive 
strains resulting from a small area of the under- 
lying metal being removed by corrosion. 

Figure 3 shows the face of a panel that was 
damaged by corrosion originating from firing 
marks on the reverse side. Several holes and also 
a number of the previously described fractures are 
apparent on this specimen. It should be pointed 
out, however, that these defects did not occur on 
the panels that had been given a thin second coat 
of enamel on the reverse side, nor had corrosion 
progressed to this extent at any location except 
Atlantie City. 

Another type of surface defect, noted only at 
Atlantic City, was a type of staining that invari- 
ably appeared on the enamel surface adjacent to 
any damaged area that exposed the metal. These 
stains appeared as iridescent films and were 
believed to be analogous to the iron stain on glass 
containers described previously by Marboe and 
Weyl [6]. These investigators found that staining 
of wet glass in contact with iron results from the 
interreaction of positively charged ferrous ions 
with the glass surface, thus leading through oxida- 
tion to the formation of an insoluble, colored ferric 
hydrosilicate. In this study, gloss measurements 
on areas showing these stains gave exceptionally 
high readings. 
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Fiaure 3. Panel A-58 after 7 yr of exposure in the salt 
air at Atlantic City, N. J. 


Surface defects, including hole at right center, were caused by rusting 
through from areas of incomplete coverage on the back of panel. Specimens 
with a thin second coat of enamel on the reverse side were not affected in 
this way. 


Numerous tiny localized areas of iron staining 
were also noted on many of the Atlantic City 
panels. Close examination of these patches showed 
tiny pin holes at the centers, which were allowing 
corrosion of the metal. 


IV. Surface Changes Resulting From 7 
Years of Weathering 


Reference has been made in the preceeding sec- 
tion to corrosion of the steel where it was not ade- 
quately covered by enamel. This part of the 
report pertains to changes that occurred as a result 
of weathering on the enamel surfaces where cover- 
age was complete. 


1. Microstructure 


As shown in figures 4 and 5, photomicrographs 
taken with a metallographic microscope indicated 
practically no change from the original surface 
microstructure of the enamels with high acid 
resistance after | and7 vr of weathering. Enamels 
of poor acid resistance, on the other hand, showed 
considerable progressive deterioration. The“ pits” 
on most of these enamels (illustrated by several 
photomicrographs in the first report [3]) were no 
longer a characteristic feature of the microstruc- 
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Ficure 4. Photomicrograph (600) of typical black, 
glossy, acid-resisting enamel showing no fundamental 
change in microstructure with progressive weathering. 


S is storage panel, T-69. E-1 and E-7 show areas on a duplicate panel, 
T-61, after exposure for 1 and 7 yr, respectively, in Washington, D. C. 
Specks are probably mill additions, 
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Figure 5. Photomicrograph (600) of a red, semimat 
enamel of good acid resistance, showing no appreciable 
change in surface microstructure from weathering during 
a 7-yr period. 

8, storage panel, P49. E-1 and E-7 show areas on a duplicate panels 
Particles 


P-Al, exposed for 1 and 7 
exposed at surface are undissolved material added to give a mat finish. The 


yr, respectively, in Washington, D. C 


difference in size of the exposed particles is not considered significant. 
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Frat RE 6. Photomicrograph x 2 of a black, nonacid 
resisting enamel showing large surface change after? y 


of weathe ring. 


S is storage panel, 1-19. E is duplicate panel, 1-11, exposed 7 yr in Wash 
ton, D. ¢ Large light areas in E are alteration pro‘ucts resulting from 
1 chemical change in the surface layer 


ture, and deterioration of the surface had pro- 
gressed considerably bevond the pitting stage. 
Two examples of these enamels are illustrated in 
figures 6 and 7. 

The photomicrographs in figure 6 provide a 
comparison between the storage panel (V—-19) 
ind the duplicate panel (V—11) after 7 yrs of 
veathering. The surface in the photomicrograph 
f V—11 is that of a relatively heavy gel layer, 
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which formed during the 7-yr exposure. This 
gel is apparently a hydrated layer, rich in silica, 
resulting from the leaching of soluble constituents 
from the enamel surface. The layer was much 
softer than the underlying enamel and could be 
scraped from the surface with a sharp blade. 





Figure 7 Photomicrograph <20) of a red, nonacid- 
resistant, semimat enamel pre pared with a two-frit combi- 
nation showing a grain boundary type of disintegration 


afte r weathe rind. 


S is storage panel P-29. E is duplicate panel P-22 after exposure for 7 yr 
it Washington, DD. ¢ Dark gray material surrounding lighter colored 
grains in S is made up of the more fusible of the two ename!l frits rh wore 


fusible frit shows greater weathering effects and appears us a light material 
bounding darker grains in E rhis light material appears under the micro 


scope to be shattered 
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Thickness measurements, as made with a magnetic 
thickness gage, before and after seraping, indicated 
a gel layer of 0.0008 in. on the Washington panel, 
V_11, and 0.0012 in. on the St. Louis panel, V—15. 
The surface film on all V-11 to V—18 panels had a 
vitreous texture and still retained an average of 
50.6 percent of the initial gloss. Fading was pro- 
nounced on these panels at all locations, the color 
changing from black when first installed to a light 
gray after 7 yrs. 

A portion of the gel-like layer was scraped from 
the St. Louis panel, V—-15, after cleaning. The 
ignition loss of the scrapings upon heating to 
1,650° F was 23.4 percent. The material after 
heating was a brown, partially sintered powder 
having no resemblance to a porcelain enamel. The 
loss in weight of the heated material was believed 
to be due almost entirely to loss of water. 

The photomicrograph of panel P-22 (fig. 7) 
shows a peculiar type of surface break-down. 
The original enamel was a glossy red of class C 
acid resistance. According to the frit manu- 
facturer, the enamel represented by panels P-21 
to P-28 was prepared from a_ two-frit combi- 
nation. From the results obtained it is apparent 
that one of the frits was much more resistant to 
weathering than the other, and the two frits did 
not blend during firing into a homogeneous glass. 
The frit of poor weather resistance is shown by the 
light areas in figure 7, E, which appear in figure 
7, 5, as dark boundaries to the relatively large 
grains. Another enamel (panels P-31 to P-38) 
prepared from the same frits, but applied by a 
different fabricator, did not show this effect. In 
this case the enamel may have been more finely 
ground and fired for a longer time. 


2, Fading 


Observations of color difference made on all 
panels during the 7-yr inspection indicated that: 


? Test for Acid Resistance of Porcelain Enamels; Part 1 Flatware. 
Issued by the Porcelain Enamel Institute, 1010 Vermont Ave., NW., Wash- 
ington, D.C. In the commercial test, which separates ename! according to 
classes, a small pool of 1-percent citric acid is placed on the specimen for 
15min at 8% F, The degree of attack is then evaluated by visual methods 
by using such characteristics as visuel strain, blurring of image, and ease of 
removal of a pencil mark. Class AA shows no visible effect from the treat 
ment and is the most resistant, with class A, class B, class C, and class ID 
following in that order. Enamels falling in the latter two classes are not 
considered as acid resistant 

A research test for acid resistance is included in the same pamphlet. The 
research test is quantitative in that the loss of 45° specular closs is measured 
for each specimen after a 15-min immersion in the 10-pereent citire acid at 
aoe PF, 


52 


1. No colored enamel of class AA or A by the 
Porcelain Enamel Institute acid-resistance test 
(see footnote 2) showed any noticeable change in 
color. 

2. No colored enamel of class B acid resistance 
showed any objectionable fading. 

3. Practically all deeply colored enamels of 
classes C and D acid resistance showed color 
change, and in most cases this fading was suf- 
ficiently pronounced to be considered objection- 
able. 

4. The colored full-mat enamels that showed 
pronounced fading at the end of the first year did 
not change much in appearance during the next 
6 yrs. 

5. Fading of the colored enamels of poor acid 
resistance was almost equally pronounced at all 
four exposure locations, even beneath the surface 
deposit that formed at St. Louis. 

Reflectance measurements for color difference 
were not attempted at the 7-yr inspection. De- 
terminations of this type are currently under 
study, and it was believed that color difference 
data taken at the time of the inspection would be 
of questionable value. A standard method of 
measuring color difference is expected to be estab- 
lished in time for the 10-yr inspection in 1950, at 
which time quantitative comparisons of stored 
panels and exposed panels may be made. 


3. Specular Gloss 


Specular-gloss measurements were made on each 
panel at the exposure site, by using the Hunter 
Multipurpose Reflectometer [7] adjusted for a 45° 
angle of incidence. Measurements were made at 
two fixed locations near the center of the panel 
immediately after the cleaning operation, which 
consisted of (a) washing with a warm 1-percent 
solution of trisodium phosphate, (b) thoroughly 
rinsing with tap water, and (¢) drying in air. The 
initial measurements were standardized 
against a liquid film [8]. The 7-yr data were 
obtained with the same multipurpose reflectom- 
eter and with the same liquid-film standards. 

The changes in 45° specular gloss that had 
occurred at the end of 7 vr of weathering are 
summarized in table 3. The values as given for 
Washington, Lakeland, and Atlantie City are for 
the most part complete, but only a few represent- 
ative values are reported for St. Louis. The St. 
Louis measurements were all made on surfaces 
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that had been cleaned by the nonstandard pro- 
cedure of a vigorous and prolonged scouring, which 
in effect constituted a polishing action that affected 
the gloss and thus invalidated the gloss measure- 
ments on these specimens. 

In general, it was noted that the 45° surface 
gloss changed at a faster rate in the earlier stages 
of exposure than later. Figure 8 illustrates this 
effect for three Washington panels that were 
chosen as being representative. In all three 
panels there was considerably more change during 
the first 3 years than during the next 4, and in the 
case of panel F-11, there was more change during 
the first vear than during the next 6. 

The relation between acid resistance and weather 
resistance, as measured by the percentage of 
initial specular gloss retained, is given in table 4. 
It will be noted that there is a direct correlation 
between the class of acid resistance by the PEI 
commercial test (see footnote 2) and the average 
percentage of gloss retained. 
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Ficure 8. Curves showing decrease in percentage of gloss 
retained with increasing exposure time for three typical 
panels exposed at Washington, D. C. 
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TABLE 4. Summary of gloss data showing the average per- 
centage of specular gloss retained after ? years of weathering 


for various classes of acid resistance 


Average percent of initial specular gloss 
retained § at 


Number of Acid re- hetate mE sh Lag Ee . Ww 
enameis sistance | | 
averaged ! class 3 Wash- Lake- Atlantic ~y 4 } 
ington land | City heentions 
29 AA 73.3 81.6 69.6 74.8 
15... A 69.7 70.0 62.3 67.3 
3 B 70.7 4.2 0.2 61.7 
24 C 9 43.1 §2.8 0.3 
16 D 44.4 36.0 46.3 42.2 


' Each enamel represented by 2 panels at each location, 

? From spot tests made on 12- by 12-in. storage panels, using the porcelain 
Enamel Institute standard acid-resistance spot test for flatware 

§ Percentage of gloss retained for St. Louis panels not included because 
tightly adhering surface deposits made gloss measurements unreliable. 


The results at Lakeland showed the same 
excellent correlation between the class of acid 
resistance and gloss retention that was shown by 
the average. At Washington and Atlantic City 
there was a reversal of two values, involving class 
B enamels in both instances. Only three class B 
enamels were included in the investigation. If 
more enamels of this class had been included 
a better average value would have been obtained 
and these reversals might not have occurred. 

Table 4 shows only moderate differences in the 
perentage gloss retained at the three locations after 
7 yr of exposure. The class AA and class A 
enamels were somewhat more affected by the 
conditions existing at Atlantic City than at 
Washington and Lakeland. On the other hand, 
the semitropical conditions at Lakeland and the 
accompanying fungus growth appear to have 
produced more surface deterioration on class C 
and class D enamels. 


V. Effectiveness of Scouring in Restoring 


Gloss 


The effectiveness of scouring with a commer- 


in restoring gloss to two weathered 


cial cleanser ‘ 
panels was investigated during the Lakeland 
inspection. A red, nonacid resistant panel (N- 


54) was found to increase from 52.8-percent gloss 


3 Examination showed the cleansing powder to consist of finely ground 
feldspar with a soap admixture 
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retained after trisodium phosphate cleaning to 
60.4 percent after 5 min of scouring with the 
powder. <A red, acid-resistant composition (panel 
L-13) on the other hand, decreased from 96.5 to 
88.3 percent. Indications were that a more pro- 
longed rubbing with the commercial scouring 
powder did not materially change these values. 
The continued, severe scouring required to clean 
the surface deposit from the St. Louis panels, on 
the other hand, showed a polishing action that 
raised the gloss readings of some of the enamels 
of better acid resistance to even higher than their 
initial values (see table 3). 
VI. Accelerated Weathering Tests 

In the earlier report [3] reference was made to 
a carbon dioxide test that produced a type of sur- 
fuee breakdown of the finish coat closely related 
to the deterioration that occurred with the poorer 
enamels during the first year of weathering. 
Gloss measurements made before and after the 
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Figure 9. Curves showing relationship between the average 

percentage of gloss retained after weathering and the per- 

centage of gloss retained after treatment of duplicate speci- 

mens by both the CO®* test and the PEI research test for 
acid resistance. 


Letters on curves show class of acid resistance by PEI spot test. The gloss 
values for class I) enamels were too low for measurement after the PEI 





research test treatment, 
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carbon dioxide treatment gave gloss-retention 
values that correlated well with the percentage of 
gloss retained by the exposed panels. 

It was hoped that this correlation might con- 
tinue with a more prolonged exposure of the 
panels, but the results of the 7-yr inspection 
indicated that the carbon dioxide test is not a 
sensitive criterion for predicting the durability 
of an enamel exposed to weather. As presently 
constituted, the test does not sufficiently differ- 
entiate between the best and the poorest enamels. 
This is brought out in figure 9, which shows the 
relationship between the average percentage of 
gloss retained by the weathered panels after 7 vr 
and the average percentage of gloss retained 
after treating the duplicate laboratory specimens 
in accordance with the carbon dioxide test and 
with the Porcelain Enamel Institute research 
test for acid resistance. It will be noted from 
these curves that the 7-yr resistance to weather- 
ing has a considerably better correlation with 
the PEI research test for acid resistance (see 
footnote 2) than with the earbon dioxide test. 

As pointed out earlier in section IV, 3, the PET 
commercial test for acid resistance (see footnote 
2) also correlates well with the average percentage 
of gloss retained after 7 yr of weathering and can 
undoubtedly continue to function satisfactorily as 
an acceptance test until such time as more com- 
plete data are available. 


VII. Discussion 


In the first report of this investigation [3] there 
was considerable discussion regarding the mechan- 
ism by which porcelain enamel surfaces are affected 
by weathering. The inspection after 7 vr has not 
appreciably altered these concepts except perhaps 
with regard to pitting. The pitting that was pre- 
viously illustrated by a number of photomicro- 
graphs is now recognized as being only an early 
manifestation of surface deterioration. On longer 
exposure, these same surfaces show an advanced 
stage of attack, an extreme case of which is illus- 
trated by panels V-11 to V-18 (see fig. 6). The 
entire surface of the enamel on these panels has 
been altered by chemical change. 

It should be strongly emphasized that the type 
of surface breakdown referred to in the preceding 
paragraph occurs only on those enamels of ex- 
tremely poor weather resistance, which can be 
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eliminated from use by the criterion of acid re- 
sistance. All enamels of class B acid resistance 
or better were still in excellent condition at the 
latest inspection. 

At the end of the first year of exposure, it was 
noted that weathering effects were more pro- 
nounced at St. Louis and Washington than at 
Atlantic City and Lakeland, and this was as- 
cribed to differences in combustion gases in the 
respective atmospheres. At 7 yr, these differences 
in severity were no longer in the same order. 
After the longer exposure time, conditions at 
Atlantic City seemed the most severe on the acid- 
resistant enamels, whereas Lakeland conditions 
were most severe on the nonacid resistant com- 
positions. The deposit that formed on the St. 
Louis panels made a reliable comparison with the 
other exposure sites impossible. Any future tests 
in such areas should include arrangements for at 
least semiannual cleaning of the specimens. 

The cause of the gradual decrease in gloss of the 
enamel surfaces with weathering, as illustrated in 
figure 8, is believed to be a slow leaching of slightly 
soluble constituents from the enamel surface, 
leaving a gel-like layer rich in silica. The rate at 
which this alteration occurs is a function of the 
composition of the enamel and the conditions to 
which the surface is exposed. After 7 yr, the 
thickness of the altered layer may be compara- 
tively heavy (up to 0.0012 in. thick), as is the case 
with enamel V—11 to V—18, or it may be very thin 
on the enamels of high acid resistance. 

This same type of surface deterioration also 
occurs with glass. According to Jones [4], when 
glasses containing less than 60 percent of silica 
are subjected to weathering, the divalent ions 
from the glass surface go into solution and are 
replaced by hydrogen ions from weakly acid water, 
thus resulting in a hydrated-silica surface layer. 
That glasses of higher silica content also may show 
this same gel formation on longer exposure is 
shown by Laubengayer [9] in his study of a soda- 
lime-silica glass (66.5 percent of silicon dioxide) 
that was entombed in a burial vault in Cyprus for 
approximately 1,800 yr. Laubengayer found 
that this glass was incrusted with a white flaky 
material that consisted mainly of hydrated silica. 
Also, other investigators [10] have found indica- 
tions of the presence of a gel layer on old window 
glass and old glass tubing. 

Fading of the colored enamels with poor acid 
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resistance is also probably caused by the leaching 
and subsequent hydration of the surface layer. 
The resulting hydrated film on colored enamels 
usually assumes a lighter color than the original 
surface and gives a faded appearance to the panels. 
The degree of color change depends on the thick- 
ness and composition of the film and on the original 
color of the enamel. Unlike the panels studied 
by Sweo [11], in this investigation enamels having 
the deeper colors, rather than the pastel shades, 
showed the maximum color change. The etching 
at the bottom edges of the panels, as reported by 
Sweo, was noted on a few panels at Lakeland. 
The method of mounting the panels was such as 
to prevent water from collecting in pockets, but 
the heavy rainfall at Lakeland (see table 2), 
together with the resulting high humidity, prob- 
ably allowed water to remain along the bottom 
edge for prolonged periods during humid weather. 

Protection of the surface of nonacid-resistant 
colored enamels by a very thin application of a 
clear overglaze enamel has proved only moder- 
ately successful. In the earlier report [3] this 
method of protecting the surface appeared to show 
marked promise, but after 7 yr, surface pits had 
begun to appear in several of these overglazed 
surfaces. Heavier applications of the overglaze 
would undoubtedly overcome this weakness. 

One of the more important observations made 
during the 7-yr inspection was the great impor- 
tance of good enamel coverage when panels are 
exposed to salt air. The salt air conditions at the 
exposure site at Atlantic City are probably more 
severe than at most seacoast installations, and in 
commercial practice the reverse sides of the panels, 
where corrosion began, would not be exposed as 
they were in this investigation. It is the authors’ 
recommendations that on all seacoast installations, 
special care be taken to insure complete enamel 
coverage of the backs of panels and attachment 
lugs. The same precaution should probably also 
be followed on tropical or subtropical units, and as 
an engineering safety factor the practice could well 
be universal. 

A discussion of the mechanism of weathering 
naturally results in bringing to the fore the enamels 
of relatively poor weather resistance. The im- 
portant fact is, however, that enamels of good acid 
resistance, applied to obtain good coverage, were 
still in very good condition in all locations after 
7 yr. Hudson and Banfield [12] working in Eng- 
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land, also found acid-resistant enamels to be in 
excellent condition after 5 yr in industrial or non- 
industrial atmospheres, or after 2 vr submerged in 
sea water. 

On the basis of both the l-vr and the 7-yr data, 
the authors feel that the recommendations as made 
in the earlier report continue to be valid. These 
recommendations were : 

1. Where appearance is an important factor, 
full-mat enamels of the type included in this 
investigation should not be used for outside 
installations, as they tend to accumulate and 
retain a dingy film and to fade. 

2. Enamels of acid resistance less than class B 
(PEI test) should not be used in any architectural 
installation where general appearance and absence 
of fading are important. An acid resistance of 
class A or class AA is to be preferred. 


VIII. Summary 


A second inspection of the 784 1-ft.-sq. porce- 
lain enameled panels of varying types exposed at 
Washington, D. C., St. Louis, Mo., Lakeland, Fla., 
and Atlantic City, N. J. was completed during 
1947. The observations made during this inspee- 
tion, which represents 7 yr of exposure, may be 
summarized as follows: 

1. Good correlation existed between acid resist- 
ance and the percentage of initial gloss retained, 
the enamels of best acid resistance retaining the 
highest percentage of their original gloss. 

2. No noticeable fading of enamels of class AA 
or class A acid resistance occurred, nor was there 
objectionable fading of class B enamels. Prac- 
tically all class C and class D colored enamels, 
however, showed very noticeable color change. 

5. In contrast to the results of the first-year 
inspection, which showed greatest weathering 
effects at St. Louis and Washington, the 7-vr 
data indicate that the conditions at Atlantic City 
were slightly more severe than elsewhere on the 
acid-resistant enamels, whereas Lakeland condi- 
tions were most severe on the nonacid-resistant 
compositions. 

4. The salt-air Atlantic City 
caused considerable corrosion of those parts of the 
panels that were incompletely covered by enamel. 
This corrosion caused failure of attachment lugs 
and in some cases failure of enamel on the face 
by rusting through to near the enamel-metal 
interface from areas of poor coverage on the back. 
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conditions at 


Specimens with a thin second coat of enamel on 
the back were not affected in this way. 

5. Where the initial coverage was complete on 
all parts of the panel and where no mechanical 
damage had occurred during exposure, protection 
of the metal against corrosion was unimpaired on 


all specimens after 7 vr of weathering. 
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On the Precision of a Certain Procedure of Numerical 
Integration 


By Harry D. Huskey 
With an appendix by Douglas R. Hartree 


An example of numerical integration is given that shows very systematic effects in the 
less significant digits. This lack of randomness gives rounding-off errors that exceed the 
predicted standard deviation by a factor of three. 

The example considered in this paper shows that systematic rounding-off errors can 
occur in numerical integration, irrespective of the number of digits kept in the contribu- 
tions to the integral. In the appendix this phenomenon is examined, and criteria are set 


up to detect the eases in which it may arise to a serious extent. 


I. Introduction 
























The use of numerical methods has led to the 
study of the accumulation of errors in computa- 
tions by various people.’ In this paper we apply 
formulas developed by Rademacher to the errors 
involved in the integration of simultaneous linear 
differential equations. The system chosen for 
this application is 


r’(t)—y(b), 


y’ (t)=—ax(t). 

The results of integrating these equations were 
easily checked by comparison with the sine and 
cosine tables published by the National Bureau of 
Standards.* 

The errors involved in the numerical integration 
of these equations arise from two sources. One, 
called the truncation error, arises from replacing 
the differential equations by difference equations; 
the other, a round-off error, comes from the 
rounding-off procedure used in the computation. 
Formulas developed by Rademacher account for 





' F. Schlesinger, Astron. J. 3@, 183 (1917); D. Brouwer, Astron. J. 46, 149 
(1987); H. Rademacher, On the accumulation of errors in processes of inte- 
gration on high-speed calculating machines, Proceedings of a Symposium 
on Large-Seale Digital Calculating Machinery (Harvard University Press, 
Cambridge, Mass., 1948). 

? Tables of sines and cosines for radian arguments (National Bureau of 
Standards, 1940) MT4; Tables of Circular and Hyperbolic sines and cosines 
for radian arguments (National Bureau of Standards, 1939) MT3. 


Errors in Numerical Integration 


the truncation error. The rounding-off error can 
be estimated in a statistical manner, provided the 
dropped digits are randomly distributed. Rade- 
macher suggests that this random property ts 
satisfied provided the increments involved in the 
integration are not too small. We shall exhibit 
an integration where this assumption is satisfied, 
but the dropped digits vary from zero to four and 
back to zero over a range involving nearly three 
hundred steps in the integration. This causes the 
error to increase by a factor of twenty and to 
become almost three times the standard deviation 
as given by Rademacher’s formulas. 

In certain other runs the error exceeds the pre- 
dicted standard deviation by a small factor. In 
two of these cases results were tabulated every 
five or ten steps in the integration and a frequency 
count of the digits taken. Standard statistical 
tests indicate that these numbers did not consist 
of randomly distributed digits. 

These results show that one must be very careful 
in applying error estimates based on an assump- 
tion of randomness. To be safe it is best to use 
the estimates for the maximum rounding-off error. 


1. Rademacher Theory 
(a) Heun Method 


We now indicate the method of solution studied 
by Rademacher and give the formulas developed 
by him. He starts with the system. 
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x’ (t)=f(zr,y) 
(1) 
y' (th =g(a,y) 


and the solution is to be found for an interval 
toSt=T by application of the Heun method. 
That is, having found 2,_, and y,_; as approxima- 
tions to the solutions at t;-;=to+(j—1)(At), the 
following formulas give x, and ¥;. 


TF yy + At f(2j-1,.Y5-1) 
(2) 
Y= Yn-rt Mteg(ayn.yy)), 
At a 
T=HotGVlawatS(y)) 
(3) 


wy +4 a late, iW-d+9(2y¥))) 


(b) Definitions 


Let us make the following definitions: 

(a) Let a(t), y(t) be solutions of eq 1 satisfying 
the condition that r(f,)—s, and y(t,)=y,. 

(b) Let x,, y;, j= 1,2, . . ., n, be the numbers 
obtained by successive application of eq 2 and 3. 

(c) Let A(t), w(t) be generic notation for solu- 
tions of the system 


d\/dt = — (Of /dr)A— (Og oa 
(4) 

dujdt= — (Of /Oy)A— (Og/Oy) u , 
(d) Let u(t)=2(t))—r, and e(t))=—y(t)—y 
j=l, ..., nn. The numbers «(t;) and e(t;) are a 


measure of the truncation error in each step of the 
integration. 


(c) Truncation Error 


Rademacher derived the following formulas for 
the truncation error: 


M(T)u(T) + p(T)r( 7) ~ r | 
( 
| 
) 


san ,(At)*[ A(t)” 


1 )+-p(ty” On 


(5) 
- (At)? if IN’ (f(t) +p’ (d)y” (t) dt 


The truneation errors u(7) and ¢(T) can be sep- 
arately obtained from eq 5 by applying the proper 
terminal conditions to the solutions A(t) and u(t) 
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of eq 4. For example, 
letting \( 7’) 


u(T) can be found by 
land u(7T)=0 


(d) Rounding-Off Error 


Thus far, it has been assumed that all computa- 
tions are done exactly. In actual computing, this 
is not the case. The accumulators or registers of 
the computing machine accommodate only a lim- 
ited number of digits. 
written as 


Thus eq 2 and 3 should be 


r=27)-,+Atf(r 


> > 


m 


+ D> tint io 


m 


1Yj-1) 4 


Y= Yjp-rt At-eg(ay,yy-a) 


> > 


m 


Lj TL; +a Uf (a; YW, ) + f(x, Yi )] t 
(7) 


+-g(tj-1,Yj-1)] 


‘ 4 
+ De ym! m 


m 


OO a 
WH=Watolg@wy) 


The €,, satisfy €,,<0.5. The coefficients r2 
depend not only upon the equations to be solved 
but upon the explicit procedure or order of oper- 
ations in the process of solution. In the following 
discussion quantities with bars above them repre- 
sent the actual numbers stored in the registers or 
accumulators of the computing machine. Although 
the analysis can be carried through using eq 6 
and 7, Rademacher makes the simplifying assump- 
tion that 

At f(7),9)) = At f(F,,7,). 
This means that he assumes that f(2;,y;) can be 
computed sufficiently accurately so that when 
multiplied by At any inaecuracies it may have are 
lost in the digits that are dropped. Thus, 
can be replaced by 


eq 7 


rs=7; +“ SLA: DAA Fy iv] 4 €, 10% 
(8) 
Ve yy o alga + g(t y- 1) d)+€210~ 


Note that if the parentheses are removed in eq 8 so 
that there are four multiplications, then there are 
four rounding-off terms, say €}m (8= 1,2; m= 1,2). 
As in the case of the truneation error let us 
make the following definition: @—2z,—7, and 
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TABLE 1. 


Sine-cosine rounding errors 
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Letting 7,=%(T) and 7,=7(7), the expression 


for the rounding-off error is 


A(T) ul T) + a | T)r(T) 


n 
— 107 Sed; +€ jou). (9) 
1 


}= 


From the inequalities | €,,, =0.5, m 
imum possible value of the rounding-off error is 


1,2, the max- 


n 
MT )ulT)+ wl Tye(T) S107 SSCA, + u,)) 


j=1 


u(t) jdt. 


107% (7? 
fees | A(t) 
Mt J, 
However, if the ¢,,.. m=1,2, are random variables 
then the standard deviation of the rounding-off 


(10) 


error is 
of (TT) u( T 


za 
At) 
3 


| u( T)r(T)|~ 
27 
| | [r7(?) H uti 


II. Example 


(11) 
\ 


1. Sine-cosine Integrations 


To check the theory developed by Rademacher 
the system 
rit y(t) 
(12) 


y(t x(t) 


was integrated on the Electronic Numerical In- 
tegrator and Computer’ The range 0.1=t=0.9 
radians was chosen as the integration interval, 
since neither function was zero in that interval. 
(While the function is near zero the increment 
At f(7,%) is small and might lead to a systematic 
effect in the rounding-off.) All computations were 
done to 10 decimal digits. About 10 values of Af 
were used ranging from 2107° to 2*«10~° A 
run “A” was made with the parentheses appearing 
in eq 8 removed; this gives four round-offs per 

A run “B” was made with the 
this gives two round-offs per in- 
tegration step. The runs are 
tabulated in table 1. 

The first entry in each rectangle in table 1 is 
the run .1 standard deviation for the respective 
angle and increment as given by eq 11. For run 
B this standard deviation should be divided by 


integration step. 
parentheses in; 
results of these 


Underneath are the residual errors (after 


two. 


'The “ENIAC ” was built by the Moore School of Electrical Engi- 
neering of the University of Pennsylvania and is now located at the 
Ballistic Research Laboratories of Aberdeen Proving Ground 
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the truncation error is removed) for the various 
A typical entry (such as the 
residual error of 15 for the cosine in run A, angle 
equal to 0.7 radians, and increment of 2 107°) is 
found as follows: 


runs and functions. 


Integration result 0.76484 19311 


Truneation error 
as given by eq 5= +-0.00000 02577 
0.76484 21888 


True value 0.76484 21873 


Residual error 0.00000 00015 

The most interesting feature in the table occurs 
in run «4 for the sine with an increment of 2 10 
For the angle changing from 0.5 to 0.6 radian the 
residual error jumps from +-0.00000 00008 — to 
—(.00000 000190. This integration was rerun, 
and results were printed more frequently. It was 
found that most of the disturbance occurred be- 
tween 0.5211 and 0.5264 radian. Table 3 gives 
the results over this range with printings at every 
five integration steps, and table 2 exhibits a 
tvpical five steps between the values of table 3. 


TABLE 2. Sample step in the sine-cosine integration 
Mf = 0.00002 0. 52250<1 <0. 52200 
¥ 
Myi-1 I > (el +-yi-1) 
‘ f 
y=» st u=t : [ro+ ] 
t : y . 
sin? cos f 


(0 19004 81273 0. 49004. 81273 0.86657 42703 0. S657 42703 

0.52250 SHO5T 4 173314 49004 8 GOSRl0 

| S658 4 sO") S4587 — 49906 5 SHH56 42805 

| 10006 4586 . 40906 4586 SH656 42801 SHH56 42801 

252 SHH 4 1 73312 49906 5 vOsl4 

| SHO55 4 4008 2TSUS {0008 2 S55 43077 

| {0008 2TSO7 {0008 2707 SHOH55 43076 SHH55 43076 

52204 S55 4 1 73310 -- 490908 2 GOS16 

| sO 4 40910 O1207 —10 0 soo 45200 

| 19910 01206 10010 01206 S64 43258 86654 43258 

22S awa 4 i 73308 499100 Y9s820 
' 

| sHOSS 4 sell 74514 4117 SSS 45458 

[ #0911 74518 #0911 74513 SHHAS 45436 SHS 45436 

“2258 SHOSS 4 1 73306 —4911 7 g9s24 

| StW52 4 40013 47819 91S 4 . 86H52 43612 

52260 . 49913 47818 49015 47818 SHO52 43611 . 86652 43611 
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TABLE 3. 





Sine-cosine integration 


St =0.00002 


Cosine 


0. 48817 72474 0. 87274 45090 
48404 O747S 
48002 17591 
49079 32802 

. 49166 43102 
402535 484909 
443540 48068 
$0427 44485 
49514 35083 


87176 (4098 
87127 60521 
S707 48234 
S7029 27237 
SOUTY YT538 
SH930 5Y143 
SONS] 12054 
$0601 2000S SOSS1 56275 
JU68S O1S0S SHOTS] GISIZ 
#40808 OVLUS 

$9705 36938 

$9704 04033 

49722 72278 

49731 30873 

49740 O7418 

49748 74913 . 86747 11525 
$9757 42358 . 86742 13995 
49766 09753 S6737 16378 
49774 T7098 S6732 18673 
49783 44345 86727 2USS81 
SOT22 Z004 
S6717 25039 
SO712 26087 


49792 11638 
40800 TSS35 
49809 4507S 
. SOT0T 28850 
. 86702 30625 


SHH9T 32313 


49818 13073 
{UK26 SOLIS 
$9835 47113 
49844 1405s 


JUSS SOUSS 


. SHH92Z 33916 
SOORT S41 
40861 4770S SHOR? BOSSO 
4YST0 14593 


4487S S1S38 SHIT? BO45T 


SHOHTT SS202 
4USS7 48055 SHOOT 40025 
4US06 14678 
40004 S1275 d SHOH57 42708 
SHOH52 43611 
StH 44454 
SHO42 45169 


SHH62 41708 


49013 47818 
49922 14513 
44930 SOTSS 
409039 47153 SHO37 45817 
s0048 13408 SOOO2 40580 
40956 TYTYUS SHH2T 40855 
1906.5 46038 SHH22 472435 
49974 12233 
49US2 78378 
$9001 44475 
SOOO) LOSTS 


SOOLT 47546 
sOO12 47761 
SHH0T 47880 
SHH02 47482 
SWOS TH51S SOSOT 47887 
MOLT 42458 


MOW265 OS353 


SO502 47757 


BERE 


SOUSST 47540 


Thus. most of the change in error occurs in an 
interval of about 0.0053 radian. This represents 
about 260 integration steps and over a thousand 
round-offs. In one half of the multiplications, the 
digit being dropped in the rounding-off process 
(see the sixth digit in the cosine values of table 3) 
changes gradually from zero up to four and back 
to zero again. 

The errors in the sine values (see column headed 
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Cosine 


0. 50034 74108 SHSS2 47255 
SOS SUNS 4845 

. K52 05738 40308 
OOH) T1433 45805 
SOO89 ST0TS 45153 
SOOTS O2673 4416 
SOUS OS21S f 43593 
5OOVS 33713 SO547 42683 
SOLOS GYLSs SO542 41086 


NOLL2 O455: ; SHAST 4060035 


50121 2U808 g SH532 S433 
5129 O5148 NY) wiz 

S138 GO0438 

MOAT 25633 57 SO517 35 
NOLSS GOTTS SH512 33 
SO164 SASTS vy SHOT 32% 
S173 20918 ry sH502 30509 
MOIS S5SY13 7 SMUT JSS24 
SOLO SOSSS N42 20061 


SOIYY LSTS3 rv SOHST 25015 


O07 SOSUS NH482 2207S 
M216 45393 S477 BUSH 
M225 10138 S472 ISO 
M233 T4833 r SMOT 16354 
W242 SY4ATS SH40i2 13074 
W251 04073 ; SO457 11507 
NO259 GSHDS S452 OND 
A268 33113 SH47 OOS 14 
MV276 Y7T5AS si442 OS5SS87 


AOISS 61953 S437 OOTT4 


M204 26208 S431 GTS74 
50302 BO593 S426 (4880 
SOSLL S4ASSS S421 GISIO 
5320 19033 SH416 SSH58 
MIS2S ASLTS SO411 S413 
W337 47273 SHH406) S208] 
50346 11318 SOL 78604 
M34 T5313 86306 T5150 
MBS 3h SHOUT TG 


0372 UOS153 rr SHOSH OTSOR 


SOSS0 66998 ry SHSS1 4129 
50389 B0T7U8 7 SUSTO HOTU 
SOB07T W453 SOS7T1 344 
A406 SS2338 iy SH3H6 52321 
MMS 21873 5 S661 48213 
5423 S45 SO356 44018 
MMS2 480903 SASS 3YT36 
M441 1247S 72 SHS46 35569 
M449 T5913 SAS4 1 SOU14 


MMAS SY2US SH856 BAST4 


E in table 3) may be wrong by up to plus or minus 
five units, since they were obtained by subtracting 
the integration results (listed in table 3) from nine 
digit values of the true results taken from Tables 
of Circular and Hyperbolic Sines and Cosines.* 

It will be noted that the last digits of the sine 
values in table 3, generally speaking, are alter- 


nately 3’s and 8’s. This can be traced to the 


4 See footnote 2. 
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alternate 2’s and 7's or 1’s and 6’s in the fifth digit 
of the cosine values. This is another warning 
against unconsidered assumption of randomness 
in the less significant digits of numbers involved in 
computations. 

Rademacher asserts there will be statistical 
independence of the dropped digits in the round- 
ing-off process provided (At/2) f(%)-1,%)-1) and 
(At/2) f(F,,7,) differ in the place 10. The example 
given here shows that Rademacher’s condition is 
not sufficient. In fact, inspection of table 2 
shows that the increments At/2 f(z,y) may differ 
in the place 10°* and yet be alike for a large 
number of integration steps in the place 10°*"!. 

There are other examples listed in table 1 
leading to large residual errors. For example, 
consider in run A the sine and cosine values for 
At=2x10™ and the sine values for At=110™. 
More frequent tabulation of results shows a steady 
increase of residual error with no such jumps as 
described above. 


The author thanks various members of the staff 
of the Ballistics Research Laboratories and of the 
Moore School who assisted in running the prob- 
lem, and extends his appreciation to L. 8. 
Dederick of the Ballistics Research Laboratories, 
whose permission and cooperation made this paper 


possible. 
Ill. Appendix 


Note on Systematic Rounding-off Errors in 
Numerical Integration 


By D. R. Hartree * 


In this paper, which summarizes the results of a numeri- 
eal study of truncation and rounding-off errors in the 
numerical solution of a differential equation by a step-by- 
step process, Huskey has exhibited a case in which 
rounding-off errors in a sequence of successive contribu- 
tions to the solution are systematically of one sign and 
approximately equal in magnitude, although the leading 
digit rounded off is the sixth significant figure in each con- 
tribution. The result is that the rounding-off errors build 
up to a total substantially greater than would be estimated 
in the basis of a random distribution of rounding-off errors 
in the individual contributions. The purpose of this note 
is to examine this situation further and to establish 
criteria for identifying the conditions in which it is likely 
to occur, so that steps can be taken to deal with it, as for 
example by carrying an extra significant figure temporarily 
in the course of the solution. 

Consider the numerical evaluation of f ydt, k decimals 
being kept in the caleulation. 


5On the stat? of the Institute for Numerical Analysis. 
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Systematic rounding-off errors occur when the leading 
digit rounded off remains the same in a number of suc- 
cessive contributions to the integral; that is, when for 
successive contributions, last integer digit of 10**'yé¢t is 
the same. When this occurs, [last integer digit of 
6(10**'y5t)]=0, or Nast integer digit of 10**+'y(dt)?]—0; that 
is 


10n—0.5< 10**!y(5t)?< 10n + 0.5, (13) 


for soine integer n. 
or other of ¢ if 


This will usually oceur for some value 


10* max | y| (dt)? >1; (14) 
it will also occur if 
10* max yg) (6t)?< 1,10. (15) 


The range Ay of y over which the inequalities (eq. 13) 
are satisfied is 
Ay= 1/[10**'(8t)"], 


and the number N of intervals required to cover this range 
is given approximately by 


N(8t) | gj) =Ay, 
80 that 
N= 1/[(104*'(8t)* |g ). (16) 


The accumulation of systematic errors is only serious if NV 
is greater than 3 or 4, that is, if 


4-104*!(8t)3 gy <1; (17) 
for this not to occur 
(6t)? > 1/[4- 10#*1 | g|. (18) 


The inequalities (eq. 13 and 17) together provide a 
criterion for identifying the situations in which accumu- 
lation of systematic rounding-off errors may be dangerous. 
Such a situation may arise in any numerical integration, 
not only in the golution of a differential equation, the con- 
text in which it was first found by Huskey. The inequality 
(eq. 17) shows how much more likely it is to arise with 
small values of the integration interval (é6t) than with 
large values. 

In the case considered particularly by Huskey, y=cos ¢, 


k= 10, 6t=2-10-5, y=sin ¢ so that eq 13 becomes 


10n—0.5< 40 sin t< 10n +0.5; 


this is satisfied for a range of t in the neighborhood of sin 
t=, which is just the region in which the phenomenon 
does occur; and it happens to be particularly marked in 
this case, since the digit which is rounded off systemati- 
eally happens to be a4 over a considerable range. Also 


10#+1(5¢)3 | | 7. 10-4, 
so that, from eq. 16, N is about 1,400, and the inequality 
(eq. 17) is very far from being satisfied; hence it is not 
surprising that the phenomenon arose in a marked form. 
The condition (eq. 18) suggests that with the value of 


k= 10, the interval length should certainly be greater than 
s= 10, 


Wasuineton,, July 17, 1948. 
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t t of Leath th Synthetic R 
a reatment of Leatner wl etic Mesins 
uw O- By Rene Oehler and Timothy J. Kilduff 
| 
Laboratory methods have been developed for treating leather with snythetie resins. 
15) The first method of treatment may be accomplished by immersing the leather specimens in 
13) monomer and subsequently polymerizing im situ; the second method consists of immersing 
the specimens in polymer solutions; and the third method involves immersion in a solution 
of partially polymerized liquid resin, which is further polymerized in the leather itself. The 
ge last method was developed in cooperation with the Thiokol Corporation, using Thiokol LP- 2. 
Abrasion resistance of vegetable-tanned crust leather treated with n-butyl metha- 
erylate, polymerized in situ, is increased by approximately 75 percent, as measured on the 
' National Bureau of Standards abrasion machine. Treatment with Thiokol LP-2 increased 
16) abrasion resistance by approximately 30 percent. However, treatment of vegetable-tanned 
oN crust leather with polybutyl methacrylate solution-polymer resulted in no improvement in 
: abrasion resistance. One-half-hour water absorption data show a reduction of one-tenth 
= : to one-third that of untreated leather, depending upon the kind of treatment and leather. 
- By proper choice of monomers and polymerization conditions, stiffness may be con- 
trolled to meet the needs of different types of leather. 
1) Water vapor permeability of shoe upper leather is reduced 60 percent, but is still eon- 
sidered ample for foot comfort. 
a 
.u- 
As I. Introduction concerning their use as impregnants for leather. 
mn, 7 a The chief reason for this lack of information is 
ie Preliminary research, conducted at the National that ordinary, commercial resins do not penetrate 
“4 nace of Standards wip boss, apueiens bee —F leather to any considerable extent. This lack of 
" acrylate monomers might be successfully poly- penetration is borne out in a paper by Niedercorn 
merized in leather within leather. However, it and Thayer [1],? who studied 25 different resins of 
t, was not possible to resume work on this project several types. They found that the maximum 
until 1946. Because of the excellent penetration resin content of treated sole leather was only 6.9 
of monomers into leather and the promise of suc- percent, Which resulted in only a slight improve- 
cessful polymerization ” situ, to give resins of ment in abrasion resistance. In fact, they found 
- Suitable molecular weights, 1 was hoped that the “that an oil treatment resulted in a greater and 
in ; : more consistent improvement than did the resin 
i resistance of sole upper leather could be improved. ecatenenin. 
50 This paper presents the first results of a ayste- Chadwick [2] reviewed the whole field of plastics 
; matic study of the general problem of impregnating as applied to leather. He proposed, among other 
F leather with resins and plastics. things, that since polvisobutylene is compatible 
BE Barer work as reste in many wellAnown igh tan wasce nay ight vata 
% wa pe “wal et nag’ Og ag rea ously replace part of the waxes and oils normally 
of ypes ’ . used in sole leather; he also suggested that certain 
n 4 This report is made as a part of the Leather Research Program sponsored thermosetting resins that were only partially 
1y the Research and Development Branch, Military Planning Division, of a ceaene- cena 
he Office of the Quartermaster General, United States Army. This program ? Figures in brackets indicate the literature references at the end of this 
4 * under the Advisory Direction of the National Research Council. paper, 
h Treatment of Leather with Synthetic Resins 63 
815197 49 5 























polymerized might be further cured by heat after 


impregnation; he suggested, too, that properly 
plasticized thermoplastic resins, such as polyvinyl 
chloride, polyvinyl chloride-acetate copolymer, and 
polyvinylidene chloride should be especially useful 
for sole leather because of their outstanding 
abrasion resistance. 

Virtue [3] has reported 50-percent improvement 
in wear resistance by impregnating finished taps 
with & to 12 percent of maleic anhydride modified 
alkyd, or modified phenolic resin. 

The wide differences in type of resins used and 
suggested by the above authors present the prob- 
lem of selecting suitable resins with which to work. 
The salient points considered in making selections 
are discussed in the next section. 

II. Selection and Properties of Resins 

The resins selected to begin this study were the 
acrylates and methaerylates, because of ready 
availability, rather wide range of properties within 
the series of esters and relative ease of poly- 
merization. The should have 
properties which will improve the wear of leather 
without seriously detracting from the desirable 
characteristics of leather. For example, a resin 
to be used for sole leather should have a greater 
abrasion resistance and may be more viscous than 
a resin to be used for shoe upper leather. Specific 
characteristics to be considered in selecting a 
resin are: proper flexibility of the treated leather, 
toughness, high abrasion resistance, high exten- 
water 


resin selected 


sibility, low water absorption, and low 


penetration. 


Tas.e lI, Prope rlies of methac rylate ester polyme rs 
Nor Ner- 
Properties =. Ethyl — A. —, 
’ pyl mityl : 
Density at 25°C ; 1.99) Lt 106 1.05 1.02 
Hardness (Pfund) at 25°C 2 Ml 100 l 20 
Thermal yield point, °C 125 65 3s w 70 
Tensil strength, Ib in. 9,000 «5,000 4,000 | 1,000 | 3,400 
Impact strength (Dynstat), ke-cm em? 10.5 wn 6.5 15 1.6 
Refractive index 1.490 | 1.485) 1.484 | 1.483 | 1.477 
Extensibility, 5-mil films (percentage | | 
at break 4 7 5 230 | 2 
Toughness, area under load-elongation 
curve, expressed in arbitrary units ys 174 76 | 1,000 | 23 





The principal properties of the common poly- 
methacrylates have been reported, as shown in 
table 1, by Strain, et al. [4]. This table shows 
that polybutyl methacrylate should be the most 
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satisfactory of the group for treating sole leather. 
Of course, no single resin or any other material 
can have all the desirable properties that may be 
enumerated. Consequently, compromises must 
be made to obtain the most desirable properties. 
These compromises can best be achieved by the 
addition of plasticizers or by polymerizing two or 
more monomers to form copolymers or inter- 
polymers with properties intermediate between 
those of the separately polymerized monomers. 


III. Methods of Treating Leathers 


The three methods used for impregnating leather 
with resins are immersion of leather in monomer 
and subsequent polymerization jn sifu; immersion 
in polymer solutions followed by evaporation of 
the solvent; immersion in partially polymerized 
resin followed by varying degrees of cure, depend- 
These methods 
have been studied and will be discussed one at 
a time. 


ing upon the desired results. 


1. Immersion in Monomer Followed by Polymerizc- 
tion in situ 

This method of impregnating leather with 
resins essentially consists in immersing the weighed 
leather specimens in a monomer containing benzoyl! 
or lauroyl peroxide catalyst for 15 minutes, wiping 
off excess solution, reweighing, transferring to a 
suitable bomb or closed system, heating overnight 
at a given temperature, evaporating the excess 
monomer, and weighing again to determine the 
amount of resin. Several variations were tried in 
the type of leather treated, the kind and amount 
of catalyst used, the temperature and time of 
heating, the different monomers or combinations 
thereof, the use of oils or other solvents with 
monomers, the use of liquid monomer (without 
catalyst) in the polymerization chamber, the 
relative size of polymerization chamber, and 
different materials for polymerization chambers 

Preliminary experiments disclosed no consistent 
relationship between the amount of resin in the 
treated sole leather and different catalyst con- 
centrations, different vegetable tannages, chrome- 
tanned specimens, and different oils used as dilu- 
ents or plasticizers, or both. However, vegetable- 
tanned crust leather does have the advantage of 
large pore size and a large void space or capacity 
to absorb monomer. A temperature range of 70° 
to 80° C was regarded as about the maximum 
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polymerization temperature that could be used 
without damaging vegetable-tanned crust leather. 
Accordingly, most of the data presented are con- 
cerned with vegetable-tanned crust leather in 
which polymerization was carried out at 70° C. 

Table 2 shows that the mean resin content 
obtained in the leather specimens decreased as 
the volume of the polymerization chamber in- 
creased from 19 to 54 ml. For chambers larger 
than 36 ml, the resulting resin content falls off 
rapidly, and all chambers larger than 54 ml are 
sufficiently large to give more or less uniformly 
low results. These low results may be due to the 
evaporation of monomer from the leather, or to 
the inhibiting action of oxygen, which has been 
reported by Barnes [5], Heuck [6], and Kolthoff 
\7|. These factors will, of course, be reduced if 
the ratio of specimen volume to total volume is 
nearly 1. Another way of reducing evaporation 
from the leather is to add liquid monomer to the 
polymerization chamber. Both ways of improve- 
ment will be presented. 

Figure 1 clearly shows the effects of adding 
liquid monomer to the polymerization chamber. 


In both groups of data, the points scatter quite 
uniformly about their respective averages. It 
should be noted that both the average and the 
run-to-run variation are considerably higher in 
the case of added monomer. However, the vari- 
ability within runs is quite uniform and the same 
for both sets of data; only two points fall outside 
reasonable range limits. 


TABLE 2. Results of in-situ polymerization of ethyl acrylate 
in vegetable-tanned crust leather heated at 70° © for 16 


hou rs 


(Composition of impregnant: ethyl-acrylate-neatsfoot oil 50:50, with 5-per- 
cent benzoyl! peroxide catalyst) 


Mean resin 


content of 
specimens 


Range of resin 


Polymerization chamber content 


Percent Percent 
Test tube, 19 ml * 13.4 to 23.9 19.7 
Test tube, 36 ml * 13.6 to 17.5 16.0 
Test tube, 54 ml * 7.4 to 15.3 11.4 
Jar, 75 ml» 3.8to 7.7 6.0 
Jar, 125 ml» 3.2to 64 4.5 
Pint jar * 25to 4.0 3.2 


* Single '4- by 3-in. specimens in each of three tubes, 




















> Three '9- by 3-in. specimens in each jar. 
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Ficure 1. Effect of liquid monomer in polymerization chamber on in situ polymerization. 


Each point represents one run consisting of three specimens, 1 by 3 in., of vegetable-tanned crust leather impregnated with ethyl acrylate-neatsfoot oil, 
50:50, and 5-percent benzoy! peroxide, heated in pint jars for 16 hours at 70° C, 
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Several attempts were made to reduce evapora- 
tion of monomer from the leather by using sand 
as a filler for the excess volume in the polymeriza- 
tion chambers. The amount of resin obtained in 
the leather specimens varied widely from run to 
run and from specimen to specimen within runs. 

Table 3 shows the inconsistent results obtained 
when a large steel bomb is used to accommodate 
full-sized soles. These variations are typical and 
were repeated several more times than shown in 
the table. Sometimes nearly all the liquid mono- 
mer had polymerized, whereas hardly any poly- 
merization occurred within the leather itself. At 
other times practically all the liquid monomer was 
recovered with almost no change in weight of the 
leather specimens during the heating period. The 
fact that liquid monomer remained in the poly- 
merization chamber apparently precludes any 
possibility of a pressure phenomenon influencing 
the polymerization. 


TABLE 3. Results of in-situ polymerization of ethyl acrylate 
in vegetable-tanned crust leather heated at 70° C (except as 
noted) for 16 hr 


[Composition of impregnant, 50.50) ethy! acrylate and neatsfoot oil] 


Specimens Mean 
Cata- treated resin 
lyst ore con- 
Polymerization chamber con- Range of tent 
. ‘ resin content 
centra- Num- of 
tion her Size speci- 
mens 
Per- Per- 
cent in Percent cent 
{ 3 )2by6 
5 jsar7e lo to 16.4 12.0 
Steel bomb (8 in. in diameter | 4 | taps. j 
’ : ! 6 12 ‘ 
by 14 in. long). | nd }2 by 6 l. 5 to 24.2 18.2 
| 7 | taps. } 
Steel bomb+30 ml monomer * 5 & 2by4 28to 194 15.4 
Steel bomb (preheated) +30 5 6 2bhy4 B7to 19s 13.9 
ml monomer.* 
| 2 4 2by4 12to2.5 1.7 
Steel bomb+30 ml monomer ' ‘ 4 2by4 3Sto lls 7.0 
| 6 4 2by4 12810253 17.4 
| 3 4 2by4 38050 4.5 
Steel bomb+30 ml monomer © 6 4 2by4 O8TOSS 4.6 
| 12 4 2by4 3.6to 15.2 91 


* Polymerized overnight at 100° C 
© All 12 specimens were in chamber simultaneously, 


Table 4 illustrates the excellent results obtained 
when there is practically no excess volume in the 
polymerization chamber over that required for the 
leather. Not only is the mean resin content higher, 
but the run-to-run variation is very much reduced. 
A larger chamber to accommodate full soles with 
a minimum of excess volume is being built to 


Table 


investigate this avenue of approach further. 
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4 also indicates that a temperature of 50° C. is 
sufficient to secure consistent results of about half 
the resin content obtained at 70° C. These lower 
results are probably sufficient for sole leather, as 
will be pointed out below. 


Polymerization in situ, using 2-in. disks of 
5-percent lauroyl peroride 


TABLE 4. 
vegetable-tanned crust leather, 
catalyst, and 2-in. pipe fittings as the polymerization 

chamber 


Polymerized at Polymerized at 
>. ar Cc 


70° € 
Mean Mean 
Impregnant Range of resin Range of resin 
resin con- — con- resin con- con- 


tent of four ‘tent of tent of four tent of 
specimens four specimens four 


per run speci- per run speci- 
mens mens 
Per- Per- 
Percent cent Percent cent 
~ 6toS.3 520 |, 
n-Butyl methacrylate 49. 9to 4.3 52.5 (25.3 to 29.7 27.6 
|s0 8to 50.2 M0 | 
45. 8to 60.7 52.2 
n- Butyl methacrylate - ethy | 
<- bee ylate - ethyl Jas gtosa.9 51.3 }20.5t0282) 227 
acrylate (¢ , 
, 146.3 t0 51.4 48.5 | 
42.9 to 44.6 43.2 
n~- Butyl methacrylate - ethyl | | 
<P geet , 46. 2to 49.2 47.7 722.2t ‘ 27.3 
acrylate (75:25). = 
. 47. 5to 611 $2.0 
18. 0to M1 aa 
n - Butyl methacrylate - ethyl | “ | 
ao 43.110 55.9 8 728. 3to31.2 20.9 
acrylate (50:50), a 
43.6 to 59.2 7 


Additional data show that oils and diluents 
alone or together have varying effeets on in situ 
polymerizations. Neatsfoot and castor oil do not 
appear to have an appreciable effect on the extent 
of polymerization, but they do act as plasticizers 
and diluents. Polyethyl acrylate by itself shrinks 
and distorts the leather specimens. However, 
when it is copolymerized with n-butyl methacryl- 
ate, or when used with about 10 percent of oil, no 
shrinkage or distortion Methyl ethyl 
ketone is satisfactory as a diluent, but Stoddard 
solvent, perchloroethylene, Marcol JX (a water- 
white mineral oil), carbon tetrachloride, dioxane, 
ethanol, and ethyl acetate inhibit practically all 
polymerization. 

Test data on treated leather 
tion IV, which compares other treatments, as well 
as this one, with untreated leather. 


results. 


will follow in see- 


2. Immersion in Polymer Solutions 


Leather may be impregnated with resins by 
means of certain polymer solutions. In order to 
use polymer solutions, the molecular weight of the 
resin and the viscosity of the solution must be 
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sufficiently low to obtain penetration into the 
lea ther. Impregnating solutions for this purpose 
may be obtained by making solution-polymers in 
a solvent, especially for use in leather, or by using 
solutions of resins, which are known to have low 
molecular weight. 


(a) Specially Prepared Solution-Polymers 


The preparation of special solution-polymers is 
quite easy, and by using different dilutions of the 
prepared resin, any desired resin content up to 50 
percent has been obtained. Briefly, the method 
of preparing these solution-polymers consists in 
placing a solution of one or more monomers with 
lauroyl peroxide catalyst in a flask, and refluxing 
until the desired resin solids are obtained. The 
extent of penetration depends upon the nature of 
the monomer used, the molecular weight of the 
polymer, the viscosity, and the solids content of 
the resin solutions. For sole leather a resin is 
desired that will add materially to abrasion 
resistance (requires moderately high molecular 
weight) and yet penetration of the resin must be 
secured (requires low molecular weight or vis- 
cosity). Hence, a compromise between high and 
low molecular weight resins must be found. 
Variations were made in a number of factors that 
influence molecular weight, such as solvents, 
temperature, catalyst concentration, monomer 
concentration in solution, ete. The prepared 
resin solutions were used by immersing the leather 
specimens for ': hour, wiping off the excess solu- 
tion, and evaporating the solvent at 80° C for 
4 to 6 hours. 

Typical changes in polymerization conditions 
and resultant effects on the amount of resin ob- 
tained in treated vegetable-tanned crust leather 
are shown in table 5. The different reflux pe- 
riods for run 1 show that resin solids or 
polymerization yield increases with prolonged 
heating. As the solids increase, the viscosity also 
increases, Consequently, run 1 shows that the 
Inean resin contents of treated specimens cor- 
responding to the different reflux periods are all 
approximately the same. The last column of 
table 5 lists the ratio of resin found in the treated 
leather to the resin calculated. The calculated 
resin is the product of the solution pick-up and 
resin solids in solution. Ratio values of 95 percent 
snd higher indicate that the resins are sufficiently 
‘ow in molecular weight and viscosity for rather 


Treatment of Leather with Synthetic Resins 


complete absorption. Ratios in excess of 100 
percent are probably due to the difficulty of re- 
moving the solvent from the resin in the interior of 
the treated leather. If one considers that the 
absorption is complete in all four parts of run 1, 
it would be expected that the resin content would 
increase with increasing solids, but the fact that 
the resin contents do not increase with increasing 
solids indicates that it is the increased viscosity 
that offsets the increased solids content. 


TaBLe 5. Preparation of solution-polymers and impreqna- 
, y Pp 
tion results 


Mean 
’ ratio of 
; Resin —— resin 
Monomer used with equa solids of ,, found 
Run volume of toluene as Re- Amount solu- COPMLEDT EG resin 
num- ~ hens otien flux ~ aly : of three aie 
ber solvent unless other- time “#talyst tion- treated ,C#leu- 
wise noted poly- eneck lated for 
mer = a three 
, speci- 
mens 
Hours Percent Percent Percent Percent 
1 | | 37.5 18. 6 loo 
3] | 39.5 19.3 1s 
l n-Butyl methacrylate 20 ‘ > : 
| 2 | 44.9 20.5 wi 
| 4] | 48.8 19.6 108 
. g 0: 
s do [oP os 30 W3 asd 
| 2 | | 43.8 4.2 55.1 
29. 22 9 
3 do ji l 5 : . 4 
| ww) | 46.1 11.0 86.1 
4 do* 1 1.0 35.5 17.1 wl 
5 Ethylaerylate 1 0.5 46.0 15.1 83.2 
6 do 1 1.0 46.3 22.8 95.1 
7 do l 2.0 w.0 23.1 “3.4 
s Ethy] acrylate: 
Unpressed specimens 1 0.2 41.2 16.2 74.4 
Pressed specimens 1 un 41.2 7.4 o4 
9 Ethyl acrylate 
Unpressed specimens 1 5 45.0 23.6 “4.3 
Pressed specimens 1 .5 45.0 16.7 &3.3 
10 Ethyl acrylate 
Unpressed specimens 1 1.0 0 30.7 92.0 
Pressed specimens 1 1.0 0 21.3 a6. 2 
ll Ethyl acrylate 
Unpressed specimens 1 20 0.0 a0 95.3 
Pressed specimens 1 2.0 0 194 au. 3 
12 Ethyl acrylate, n-butyl 
methacrylate (50:50) 1 1.0 0.1 16.3 ws 
37.1 
13 Ethyl acrylate, »-butyl 2.4 s.0 
: “ 1 1.0 {diluted 
methacrylate (20:80) | ge | 12.6 wi. 6 
50.50 
“4 Ethyl acrylate, n-butyl 
methacrylate (25:75) 1 0.5 37.6 16.6 4.0 
15 Ethyl acrylate, n-butyl 
methacrylate (25:75) * 4 5 31.3 15.4 wi. 


* Benzene solvent used. 
>» Ethyl acetate solvent used. Extra solvent necessary to keep resin 
dissolved. 


Run 2 of table 5 shows that when using 0.2- 
percent catalyst instead of 2.0 percent as in run 
1, the molecular weight and viscosity are consider- 
ably increased so that only about half as much 
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resin is absorbed, and the ratio of resin found to 
resin calculated is reduced from 100 to 88 percent. 

Run 2 also illustrates the difference in results 
obtained by treating leather with a solution- 
polymer of different solids content. Both the 
resin content of leather and the ratio of resin found 
to resin calculated are lower when treated with a 
solution-polymer of higher solids content. The 
reduction of both values points to increased vis- 
cosity directly or indirectly as an explanation for 
the lower resin content. Since Strain found that 
prolonged heating of polymethyl methacrylate 
solution-polymers decreases rather than increases 
molecular weight, then the lower resin content 
of leather corresponding to the higher solids con- 
tent cannot be attributed to an increase in molec- 
ular weight during the 20-hour reflux period. 
The ratio of resin found to resin caleulated corre- 
sponding to the higher solids content is only 55 
percent, or about 30 percent less than for the 
solution-polymer of lower solids content. This 
wide difference in ratios suggests some preferential 
absorption of the solvent for the solution-polymer 
of higher solids content. = 

Runs 8 through 11 appear to show the effect of 
decreasing molecular weights with increasing cata- 
lyst concentration in a different way. Compari- 
sons were made between leather specimens pressed 
at 3,000 Ibin. and unpressed specimens. The 
pressed specimens had a smaller volume, smaller 
pore size, and consequently less room for absorb- 
ing resin solution. It would be expected, there- 
fore, that the resin content of the pressed speci- 
mens would be less than that for the unpressed. 
More important, however, is the apparent relation 
between catalyst) concentration or molecular 
weight and the ratios of resin found to resin cal- 
culated. In the case of 0.2-percent catalyst, this 
ratio for the unpressed specimens is about 20 
percent less than that for the other three resin 
solutions of this series (runs 9 through 11). In 
other words, a larger proportion of the resin solu- 
tion pick-up consisted of solvent for the resin pre- 
pared with 0.2-pereent catalyst, which tends to 
show that the molecular size for viscosity of this 
resin solution was too large for absorption. The 
ratios for the remaining three resins of this series 
are all about the same, which also tends to show 
that 0.5-percent catalyst: produces a molecular 
size that is sufficiently small to allow very good 
absorption by the unpressed leather. This one- 
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step difference in ratios for unpressed leather further 
tends to show that a molecular size smaller than 
that obtained with 0.5-percent catalyst does not 
materially alter the ratios obtained. 

The results obtained from specimens that were 
pressed at 3,000 Ib/in’ prior to treatment with 
the same solution-polymers suggest additional 
relations between the molecular size and_ vis- 
cosity, resin contents, and ratios of resin found to 
resin calculated. In the case of pressed specimens 
treated with resin prepared with 0.2-percent 
catalyst, only slightly more than half of the 
calculated resin was actually absorbed by the 
leather. As the catalyst concentration increases 
and the molecular weight of the resin decreases, 
the ratio differences between pressed and un- 
pressed specimens exhibit two downward steps. 
These steps appear to indicate that as the molee- 
ular size and/or viscosity decreases, the large 
pore size of unpressed crust leather becomes less 
and less necessary for good absorption. I+ 
should also be pointed out that no substantial im- 
provement in ratio differences was obtained for 
resins made with 2.0-percent catalyst as com- 
pared to 1.0 percent, as shown by the results of 
runs 10 and 11. 

Run 13 illustrates the effect of diluting a some- 
what viscous resin in order to obtain treated 
leather of different resin content. It will be 
noted that slightly more than half as much resin 
is obtained in leather when it is treated with a 
solution-polymer that has been diluted with an 
equal volume of solvent. This is explained on 
the basis of reduced viscosity. The solution- 
polymer was diluted with an equal volume of 
toluene. The resin content obtained, however, 
was slightly more than half that obtained for the 
undiluted solution-polymer. A suggested ex- 
planation for this deviation from proportionality 
is based on reduced viscosity, and a higher pro- 
portion of solution solids absorbed by leather 
treated with the diluted solutions, as shown by the 
values for the ratios of resin found to resin cal- 
culated. It is apparent, therefore, that the resin 
content of treated leather may be varied at will 
to any value lower than that obtained with un- 
diluted solution-polymer, although the actual 
value of resin content obtained with the diluted 
solution may not be predicted precisely from the 
extent of dilution. 

Similar relations and effeets may be observed in 
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the remaining data of table 5, which have not been 
discussed in detail. It should be remembered, of 
course, that when prepared under corresponding 
conditions, polyethyl acrylate is softer and its 
solutions less viscous than polybutyl methacry- 
late. As would be expected, copolymers will have 
properties intermediate between those of the 
separately prepared polymers. 


(b) Solutions of Resins of Known Molecular Weights 


In order to have some definite reference points 
from which to judge the effects of molecular weight 
variations on resin impregnating solutions, poly- 
methyl methacrylates of known molecular weights 
were used for treatment, as indicated in table 6. 
The impregnations were carried out in the same 
way as those described above for specially pre- 
pared solution-polymers. Viscosities of the solu- 
tions increased with molecular weight, as would be 
expected. A resin of corresponding molecular 
weight but softer and less viscous would probably 
be absorbed better than polymethyl methacrylate. 
However, these polymethyl methacrylate solu- 
tions do demonstrate the same relations between 
molecular weight and resin content of treated 
leather, as discussed under specially prepared 
solution-polymers. Table 6 shows that when 
polymethyl methacrylate is used, as described 
above, and has a molecular weight of 11,000 or 
16,000, about 14 percent of resin is absorbed. 
However, for a polymer with a molecular weight 
of 43,000, only about one-tenth as much resin is 
absorbed. The use of methanol indicates that it 
is only slightly helpful in reducing viscosity and 
increasing resin content. 


TaBLe 6. Impregnation with polymethyl methacrylate resins 
of known molecular weight 


Mean Ratio of 


| = resin 
Mol | I omnes found to 
Molecular weight . , tesin rag resin ca!- 
of resin Solvent solids — culated 
epeci- for three 
moe speci- 
—_— mens 
Percent Percent Percent 
11,000 Toluene » 15.3 110 
16,000 do ‘ ‘ a» 13.6 107 
45,000. do » 1.4 oo 
43,000 Toluene-methanol (3:1) on 20 24 47 
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3. Immersion in Thiokol LP 2 Followed by Special 
Cures 


Another type of resin solution that has promise 
for treating sole and shoe upper leathers is a 
polysulfide rubber dissolved in’ toluene. The 
Thiokol Corporation has made exploratory treat- 
ments of leather specimens with Thiokol LP. 2. 
In considering the test data that follow, it should 
be recalled that these experiments are preliminary 
ones and do not represent any development work 
whatever. The processing appears to be simple 
and after the elimination of certain difficulties, 
should lend itself to easy commercial application. 

Essentially, the method depends upon impreg- 
nating the leather specimens with a low molecular 
weight resin solution, by immersing, wiping off 
the excess solution, and curing or further poly- 
merizing jn situ. For shoe uppers and other soft 
leathers, a “‘chain-terminator” is used to maintain 
a relatively low molecular weight and soft resin, 
which does not materially reduce flexibility. 

Treated vegetable-tanned crust leather con- 
tained 34-percent Thiokol LP-2; factory finished 
taps, 18 percent; and degreased Army retanned 
upper leather, 50 percent. Contrary to previous 
experience, treatment with Thiokol LP 2, as de- 
scribed, contributes practically no odor to the 
leather specimens. 


IV. Properties of Treated Leathers 


Comparative tests were made on different types 
of treated and untreated leathers. The tests in- 
cluded abrasion, water absorption, water-vapor 
permeability, and relative stiffness. 

All tests were carried out in accordance with 
Federal Specification KK-L-311, Leather and 
Leather Products; General Specifications (Meth- 
ods of Sampling, Inspection, and Tests). 


1. Abrasion Resistance 


Results of abrasion tests are shown in figures 2, 
3, and 4. Figure 2 shows comparative results of 
vegetable-tanned crust leather that was treated 
with solutions containing equal amounts by 
volume of the monomers indicated and neatsfoot 
oil followed by polymerization im situ. The 
improvement shown by extrapolating the curve 
for polybutyl methacrylate to the line correspond- 
ing to 60-percent abrasion loss is about 75 percent, 
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REVOLUTIONS, ABRASION MACHINE 
Figure 2. Effect of polybutyl methacrylate and polyethyl 
acrylate (polymerized in situ) on abrasion resistance of 
vegetable-tanned crust leather. 


Mean values of four specimens for each treatment plotted. ©, n-Buty! 
. ethy! acrylate (8.4%); X, untreated. 


methacrylate (20.10); 
and about 10-percent improvement for polyethyl 
acrylate. This wide difference in abrasion im- 
provement is probably due to the fact that poly- 
ethyl acrylate is much softer than polybutyl 
methaerylate, and also because there was only 8 
percent of the former resin in the treated leather 
compared to 20 percent for the latter resin. Figure 
3 shows that the improvement found was not due 
Figure 4 shows the improve- 
ment in abrasion Thiokol 
treatment of factory-finished sole leather. The 
improvement of specimens treated with Thiokol 
LP 2 is about 30 percent. Attention is drawn to 
the fact that for all the curves the points are quite 
Indirectly this linearity 


to the neatsfoot oil. 
resistance due to 


clos’ to a straight line. 














0200 
wo 
w 0.150 60%, ABRASION LOSS —> ona a 
2 
S 0100 |} —-_____- 
+ A | 
3 
o” 
2 
§ a0so |— 

P J 
° 4000 8000 12000 16000 


REVOLUTIONS, ABRASION MACHINE 


Figure 3. Effect of pelybutyl methacrylate (polymerized 
tn situ) and neatsfoot oil on abrasion resistance of vege- 
table-tanned crust leather. 


Mean values of four specimens for each treatment plotted. @, Buty! 


methacrylate (20.7°7); , neatsfoot oil: X, untreated, 
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REVOWTIONS , ABRASION MACHINE 
Ficure 4. Effect cf treatment with Thiokol LP-2 on 
abrasion resistance of factory finished sole leather 


Mean values of three specimens for each treatment plotted.@, Factory 
finished; ©, leached, factory finished; §§, Thiokol treated, factory finished; 
©, leached, Thiokol treated, factory-finished. 


implies that penetration of the resins is more or 
less uniform throughout. Uniform penetration is 
to be expected for leather treated with monomer 
and subsequently polymerized, as the monomer 


has about the same viscosity as kerosine. The 
linear abrasion curves for soles treated with 
Thiokol LP 2. indicate uniform penetration, 


thereby showing that the molecular weight is 
suitable for adequate penetration. 

Abrasion tests on leather treated with solution- 
polymers of specially prepared polyacrylates, 
polybutyl methaerylates, or copolymers thereof 
do not show any improvement in abrasion resist- 
ance over that of untreated leather. This differ- 
ence from the improved abrasion resistance shown 
by the same polymerized monomers by the in situ 
method, suggests that in the latter process a higher 
degree of polymerization has occurred so as to 
impart treated 
leather. 


greater abrasion resistance to 


2. Water Absorption 


Curves for water absorption are shown in 
figures 5, 6,7, and 8. The resins in the specimens 
shown in figure 5 were polymerized in situ. The 
curves for polybutyl methacrylate and copolymer 
of butyl methacrylate and ethyl acrylate prac- 
tically coincide. The relative position of the 
curves for these two resin treatments of vegetable 
clearly establishes the superior 


erust” leather 
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Ficure 5. Effect of resins (polymerized in situ) indicated 


on water absorption of vegetable-tanned crust leather 


Mean values of four specimens for each treatment plotted. @, Controls; 
@ ethyl acrylate (40.6%); X, butyl methacrylate (43.3%); D, ethyl acrylate- 
neatsfoot oil (21.0% KR, 22.0% ; ©, butyl methaerylate-ethyl acrylate 


copolymer) (46.0%). 
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Figure 6. Effect of solution-polymer resins on water 
absorption of vegetable-tanned crust leather. 


Mean values of four specimens for each treatment plotted. , Controls; 
@. polyethy!l acrylate (6.8%); DB, polyethyl acrylate (15.9%); @, polybuty! 
methacrylate (8.3°; 
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Figure 7. Effect of different resins on water abscrption of 
vegetable-tanned crust and factory finished sole leather. 
Mean values of two specimens for each treatment plotted , Vegetable 


rust (controls); D, vegetable crust (ecpolymer, 10°); A, factory finish (con- 
trols); @, vegetable crust (Thiokol, 34%); @, factory finish (Thiokol, 18%) 
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Figure 8. Effect of different resins on water absorption of 
degreased chrome retan upper leather. 


Mean values of two specimens for each treatment plotted , Degreased 
chrome retan (controls), §§, degreased chrome retan (Thiokol, 50%): & 
degreased chrome retan (polyethylacrylate, 50% 


resistance to water absorption by polybutyl 
methacrylate even when copolymerized with ethyl 
acrylate in equal amounts. The ethyl acrylate 
and ethyl acrylate-neatsfoot oil treatments have 
curves that are similar in general shape. However, 
the water absorption curve for ethyl acrylate is 
generally better than for the resin-oil treatment, 
and especially better for the ':-hour period. 

Figure 6 illustrates the effect of treating vege- 
table-tanned crust leather with solution-polymers. 
It will be noted that polybutyl methacrylate is 
about five times as effective in decreasing resist- 
ance to water absorption as polyethyl acrylate for 
approximately the same resin contents. Com- 
pared to the curves in figure 4, it will be observed 
that the reductions in water absorption are much 
less, as are the amounts of resin in the treated 
specimens. The ':-hr absorption value for ethyl 
acrvlate-neatsfoot oil (12 percent of resin and 22 
percent of oil) in figure 4 is the highest of the 
treated specimens, and yet it is only about two- 
thirds the value of that obtained for 16 percent 
polyethyl acrylate solution-polymer. The more 
favorable results may be partially due to the 
neatsfoot oil or to a presumably higher degree of 
polymerization by the in situ process. 

Figure 7 shows the effects of different solution- 
polymers and Thiokol LP-2 on vegetable-tanned 
crust leather and factory-finished sole leather. 
The ‘:-hr absorption values for vegetable-tanned 
crust leather treated with n-butyl methacrylate- 
ethyl acrylate solution-copolymer and Thiokol 
LP-2 are shown to be about the same. However, 
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as time increases to 24 hours, the specimens 
treated with Thiokol LP-2 absorbed only about 
half as much as those treated with the acrylate 
solution-copolymer. A substantial part of this 
difference in absorptions for longer periods may 
be due to the fact that there was 34 percent of 
Thiokol LP-2 in the leather but only 10 percent 
of the copolymer. From a service viewpoint, the 
.-hr absorption values are probably more signifi- 
cant than the 24-hr values. Factory-finished sole 
leather treated with Thiokol LP-2 has the most 
favorable absorption curve. 

Figure 8 shows the effect of different resin treat- 
ments on degreased chrome-retanned upper leath- 
er. The ‘s-hr absorption value for polyethyl 
acrylate is approximately two-thirds that of spec- 
imens treated with Thiokol LP-2. There was 
approximately 50 percent of resin in both sets of 
treated specimens, so that direct comparisons as 


TABLE 7. 


Type of leather Preat ment 





the two treatments are 


to the effectiveness of 
valid. 


3. Water-Vapor Permeability 


Water-vapor permeability data on vegetable- 
tanned crust leather (table 7) show that 8.3 per- 
cent of polybutyl methacrylate and 15.6 percent 
of polyethyl acrylate both reduce the permeability 
from 567 to about 325 (g/m*)/24 hr. Seven per- 
cent of polyethyl acrylate reduces permeability by 
an intermediate amount; i. e., 480 (g/m*)/24 hr. 
The last three items of table 7 show degreased 
chrome retanned upper leather treated with 
polyethyl acrylate and Thiokol LP-2; both de- 
crease water vapor permeability by about the 
same amount, i. e., to about 400 (g m*)/24 hr. 
This value is considered to be ample for shoe 
upper leather. 


Physical tests on treated and untreated leathers 


Mean water absorp- Mean 
Mean 
Mean tion of four speci- water relative 
resin con- mens for Vapor 
stiffness 
tent of permea- 
factor of 
four bility of fans 
specimens 14 hr 24 hr four specimens 


specimens 


Percent Percent Percent (a/m®)/24 hr I 
Vegetable crust Ethyl acrylate * w.6 8.2 22.9 
Do Ethyl acrylate-Neatsfoot * 12.0 18.0 22.2 
Do Ethyl acrylate-Buty! methacrylate * 16.0 3.4 11.2 
Do Butyl methacrylate * 43.3 5 11.2 
Do None 00 71.7 72.8 
Do Ethyl acrylate solution-polymer 1.9 31.2 43.8 316 14.5 
Do do 6.58 62.5 74.3 sy) 20.1 
Do Buty! methacrylate solution-polymer 8.3 26.0 0 331 45.7 
Do None 0.0 67.1 76.6 7 1s. 1 
Factory finished vegetable-tanned do 0 b 25.7 35.3 
Do Thiokol LP-2 solution 617.6 b4.1 » ISS 
Vegetable crust None oo 4.6 b 69.3 b7 
Do Thiokol LP-2 solution b 34.1 bO.6 » 26.4 b 23.2 
Do Ethyl! acrylate-Butyl methacrylate copolymer solu- 610.4 44 » 48.6 
tion-polymer. 
Degreased chrome retan upper None 0.0 » 22.9 bOL4 b 1011 64.9 
Do Thiokol LP-2 solution b M1 » 16.7 649.3 » 408 7 
Do Ethyl! acrylate solution-polymer b 49.5 9.5 b 42.1 » 34 1.3 


*Polymerized in situ, 
bMean of two specimens. 


4. Relative Stiffness 


The relative stiffness factors of several leathers 
and treatments are shown in table 8. It will be 
noted that 16-percent polyethyl acrylate has prac- 
tically no influence on stiffness of vegetable-tanned 
crust leather, whereas 8 percent of polybutyl 
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methacrylate more than doubles the stiffness fac- 
tor. Table 7 shows that treatment with Thiokol 
LP 2 triples the stiffness factor of vegetable- 
tanned crust leather. However, degreased chrome 
retan upper leather shows no appreciable change 
in stiffness when treated with soft acrylate or 
Thiokol LP. 2 resins. 
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V. Summary 


1. A laboratory procedure is described for poly- 
merizing acrylate monomers in situ. By use of 
small polymerization chambers, as much as 50- 
percent resin can be obtained in vegetable-tanned 
crust leather. 

2. Specially prepared solution polymers have 
been made to impregnate leather. Up to 25-per- 
cent resin has been obtained in vegetable-tanned 
crust, and 50 percent in chrome retanned upper 
leather by this procedure. 

3. Treatment with Thiokol LP-2 can be con- 
trolled to give a very stiff sole or to produce a soft 
shoe upper leather. 

4. Abrasion resist..ace of vegetable-tanned crust 
leather is improved about 75 percent by n-butyl 
methacrylate polymerized in situ. Treatment 
with Thiokol LP-2 improves abrasion resistance 
by about 30 percent, whereas solution-copolymers 
of n-butyl methacrylate and ethyl acrylate con- 
tributed no improvement. 

5. Half-hour water absorption data on vegetable- 
tanned crust leather indicate that treatment with 
n-butyl methacrylate, polymerized in situ, shows 
95-percent reduction over untreated leather. On 
the same basis, solution-copolymers of n-butyl 
methacrylate and ethyl acrylate as well as Thiokol 
LP-2 are all about equally effective. For de- 
greased chrome-retan leather, polvethyl acrylate 
solution polymer reduces }.-hr water absorption 
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about one-half, and Thiokol LP-2 reduces ‘:-hr 
water absorption about one-third. 

6. Water vapor permeability of degreased 
chrome-retan upper leather is reduced by 60 per- 
cent. However, the reduced value of about 400 
g/m?/24 hr is considered ample for comfort. 

7. Stiffness of vegetable-tanned crust leather 
treated with n-butyl methacrylate polymerized 
in situ is not appreciably different from untreated 
leather. Vegetable-tanned crust leather treated 
with Thiokol LP-2 has a stiffness factor somewhat 
more than double that of untreated leather. De- 
greased chrome-retan upper leather shows no ap- 
preciable change in flexibility when treated with 
soft acrylate or Thiokol LP-2 resin. However, a 
different cure for Thiokol LP-2 in treated chrome 
retan upper leather reduces the stiffness slightly, 
as does polyethyl acrylate solution-polymer. 
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Noise Spectrum of a Diode With a Retarding Field’ 


By Jacob J. Freeman 


A general expression for the noise-power spectrum generated by the random emission 


of electrons of arbitrary trajectories within a waveguide is obtained. 


A numerical solution 


of the potential distribution within a plane diode is effected, and a family of curves from 


which the potential distribution may be obtained is graphed. 


tween several actual potential distributions and Von Laue’s solution is given. 


A graphical comparison be- 


The above 


results are utilized to derive the equivalent mean-square fluctuation current due to the space 
charge within a diode for two cases of potential distribution, namely, the linear distribution, 


and that distribution which occurs at the neighborhood of the beginning of the retarding 





field. For the ease of the linear potential distribution, the equivalent noise temperature of 


the diode conductance is shown to be equal to the cathode temperature. 


I. Introduction 


A considerable increase at high frequencies in 
both the shot noise and conductance of a nega- 
tively biased diode over the theoretical values 
computed on the basis of zero transit angle has 
heen reported in the recent literature by several 
investigators. C. N. Smyth [1]? observed ex- 
cessive damping and noise from a negatively 
biased diode at 3,300 Me/s, and suggested it was 
due to currents induced by space charge. A. Van 
der Ziel [2] and A. Versnel [3] reported preliminary 
measurements on the equivalent noise current, 
and on the conductance of a negatively biased 
diode at wavelengths of 5.8 and 7.25 m, over a 
range of anode voltages. 

In this paper the noise-power spectrum of a 
diode with a retarding field * will be investigated 
theoretically, taking into account only the noise 
power generated by those electrons with insuffi- 
cient energy to reach the plate, in order to be able 
io evaluate the equivalent mean square current 


fluctuations due to space charge alone.’ Because 


A dissertation submitted to the faculty of the Graduate School of Arts 
nd Sciences of the Catholic University of America in partial fulfillment of 
he requirements for the degree of Doctor of Philosophy 

? Figures in brackets indicate the literature references at the end of this 
maper, 
) By a retarding field is meant one in which the slope of the potential curve 
slways negative. 
‘ For brevity, the term “‘space charge”’ will henceforth be used to designate 
© aggregate of electrons that have insufficient energy to reach the plate. 
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of the mathematical intractability of the potential 
distribution, the analysis is carried out only for 
the case of large saturation current, when the 
slope of the potential curve at the anode is small, 
and the case of small saturation current, when the 
potential distribution is essentially linear, For 
the latter case, the equivalent temperature of the 
conductance is shown to be equal to the cathode 
temperature. 

Coincidentally, the noise power spectrum gen- 
erated by electrons of preseribed trajectory within 
a waveguide of arbitrary cross section is developed, 
partly because the results can be immediately 
applied to the specific problem of a diode across 
a transmission line, and partly because the results 
are necessary for intended investigations of noise 
from transit-time devices within waveguides. 


II. Noise Spectrum 


In this section ° will be derived an expression 
for the noise power generated by random emission 
within a waveguide, or a transmission line. The 
latter is to be regarded as a special case of the 
former, in which the particular mode propagated 
is the TEM mode, and the word “guide” will be 
used to designate beth. Although in this paper 
only the particular case of a retarding-field diode 


§ MKS units are employed throughout this paper. 
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connected to a transmission line will be consid- 
ered in detail, the noise spectrum of a diode 
inserted within a eylindrical guide of arbitrary 
cross section will be treated, since the generaliza- 
tion involves little further complication, and it 
seems desirable to have an explicit treatment of the 
problem for purposes of future reference (e. g., 
computation of the noise power generated by 
electronic transit-time devices within waveguides). 
It is assumed that the cathode and anode surfaces 
of the diode coincide with parts of the guide sur- 
face, so that there are no discontinuities due to 
the electrode structure. In other words, elec- 
trons emitted from a given portion of the guide 
surface describe a trajectory that is prescribed 
for electrons with given initial velocity by a given 
direct-current potential and impinge on another 
portion of the guide surface. It is assumed that 
the given direct-current potential has no minimum, 
so that the results of this section will apply either 
to the case of a temperature-saturated diode, or 
one with a retarding field, in which the anode 
voltage is sufficiently below the cathode, so that 
the slope of the potential curve at the anode is 
negative. If no direct-current potential minimum 
exists, the perturbations of the trajectory of elec- 
trons with a given initial velocity due to small 
fluctuations in the space-charge potential will be 
small, and it is assumed that they can be neg- 
lected. If a potential minimum exists (fig. 1), 
there is a discontinuous change in trajectory for 
electrons having just enough initial energy to 
reach the potential minimum, so that a small 











CATHODE ANODE 
= 
> 
x —. 
Fiaure 1. Voltage distribution of diode with minimum. 
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fluctuation of potential may cause a large fluctua- 
tion of trajectory for this class of electrons. 

Accordingly, then, it may reasonably be assumed 
that each electron with a given initial velocity 
traverses a fixed trajectory independently of all 
other electrons; thus, one may regard the result- 
ing electromagnetic pulses generated by each 
electron as random events. 

The guide is considered to have its axis in 
the z-direction, and to be terminated at either 
end by a matched load, with the diode in the 
The immediate problem is to compute 
the average noise energy absorbed by each load, 
due to electrons having an arbitrary trajectory. 
Only the case where the range in frequency is such 
that the guide supports only the dominant mode 
will be considered. 

Then if E=E, exp(jot) represents the electric 
field of the dominant mode generated by a current 
density J=J) exp(jof), it is shown in appendix | 
that 


center. 


E, -= bir) | Jir’)-b(r exp —Jja 2’—2z \dr’. 
. (1 


Here w is the magnetic permeability, and dr’ is an 
element of volume occupied by the current density. 
w=2rf is the angular frequency, a the phase con- 
stant for the dominant mode, and 6 is the corre- 
sponding normal mode [4]. Some properties of b 
are briefly discussed in appendix 1. r(z,y,2) and 
r’(x’,y’,2’) are the coordinate vectors of the field 
and source points, respectively. If the current 
density J(r’) has negligible cross section, J(r’) 
dr’=I(r’)ds’, where I(r’) is the current flowing 
at (2’,y’,2"), and ds’ is the incremental displace- 
ment along the path of current flow, s’, corre- 
sponding to dr’. Equation | then becomes 


Mw F , , Pian ‘als. al os, o’ 

E\(r) — pe bir) | I(r’)ds’ -b(r’ \exp[ja(z’—2)],2>2’. 

(2) 

When the current /(x’,y’,2’) consists of an electron 
of charge e following the trajectory s’, 


I=¢6(t—t’) = | exp[jw(t—t’)]dw, (3 


where ¢’ represents the time at which the electron 
is at (2’, y’, 2’), and 6(t—t’)=0, for t#t’, and 
Si(t—t)dt=1. Accordingly, the field generated 
by an electron is 
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‘on 
nd 
ted 


‘ch 





. ueb Pe dw ; ‘ a @ 
E(t) =— > | og XP [j(wt—az)] ~ 
"to ; 
v(t’)-6 exp [j(a2’—et’)|dt’, (4) 


where the variable of integration is changed from 
tot’, v(t’) =ds'/dt’, and ¢, is the transit time of 
the electron. 


Let 
t . 
8(w) ( v(t’)-b exp [j(az’—at’)|dt’. (5) 
Then, 
™ eb? w , P ; 
E(t ae = {- ~ xp [j(wt—az)|s(w)dw. (6) 


Since H=j/uw curl E, one finds similarly that 


» -2i ‘exp[j(wt—az)]8(w)dw 
[ ~aecone b}. 


a 


(7) 


Since it is desired that E(t) be real, @ is defined 
as an odd function of w, so that the imaginary 
parts of the integrals in eq. 6 and 7 vanish, since 
they are odd functions of w. 

The instantaneous power flow due to the pas- 
sage of a single electron, is, by Poynting’s 


Theorem, 


Pit EE H-enda. (S$ 


~S 


where nis a unit normal to the surface S’. 

If one utilizes the orthogonality properties of the 
normal modes, as stated in eq. 94 and 95 in 
appendix 1, namely, 


(_ ee 


and 


[b,x (Vrb,]-kda=0, 


it is found on substituting eq 6 and 7 in eq 8, that 


pe = w . ? 
Pit —J | exp (3(wt—az))s(w)d 
162° Jo @ P ste az) )O\e0 0 
exp [)(wf—az)]s(w)dw. (9) 


Since the pulses of energy are random, or statis- 
tically independent, the average power flow, by 
’ , ry - ° , 
ampbell’s Theorem [5], is 


{ 
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P=P(t)=n Pi(t)dt, 


. 


where \ is the average number of electrons emitted 
per second, and Xe=/, the direct current. 


= hew [{° ‘ 
I am) | 


| 8(w) exp [j(wt—az)]dwdt. (10) 


w . 
exp [j(wt—az)|s(w)dw 
a 


By Parseval’s Theorem [6], 


| A (Of Ftdt Fi(f) F(t fodf, (11) 


where f; and f, are the Fourier transforms of F, 


and F),. 


Substituting eq 11 in eq 10, 


S(@)S w)dw. 


P he“ u ( w 


St a 


Or, since s(—w)=s*(w) where s* is the complex 


conjugate of s, and since s*s is even in a, 
ple ( w 
4m Jo 


S(q@) 2dw. | 12) 
a 


If one defines wu/a=Z, the wave impedance, then 


el 


4dr, 


Z(w)'8(w) *dw. 13) 





We consider now a plane parallel transmission 
line, with planes a distance @ apart, and of width 
d (fig. 2). For this ease, b=i/ yad, and a=wy ue, 
where iis a unit vector in the z-direction, and € is 
the specific inductive capacity of free space. 
Hence, 


— el oo . 
P raya 8(w) *dw, (14 


J 
where 
( v(t’)-iexp j(w y ue =’ — wt’) 
S(w) 
J0 yad 


If the trajectory is considered to be in the s,y, 
plane, and the linear dimensions of the cathode 
small compared to the wavelength corresponding 
to any frequency in the range considered, then 
one may take 2’=0 and let the origin of coordi- 
nates be at the cathode. Equation 15 becomes 
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*’p(t’)-d exp (—jot’)dt’ 


(16) 
J yad : 


S(qw) 


In order to compute the power spectrum from 
a diode with a retarding field the trajectories must 
be computed, and this requires that the potential 
distribution for the case of a retarding field be 


determined. This is indicated in the next section. 











Ye 


Figure 2. Schematic of parallel-plare transmission line. 


III. Potential Within a Plane Parallel 
Dicde Subject to a Retarding Field 


1. Rigorous Solution 


The potential distribution within a conducting 
diode has been thoroughly investigated by Fry [7] 
and Langmuir [8], among others. However, Fry 
limited his investigations to the case where the 
potential has a minimum. Von Laue [9] has 
exhaustively treated the potential distribution of 
a diode for the case of a retarding field, where it is 
assumed that the plate current, 7, makes a negli- 
gible contribution to the potential in comparison 
to the contribution made by those electrons that 
never reach the anode. In other words, Von Laue 
assumes that the electron gas between the two 
electrodes is in statistical equilibrium, and that 
the space-charge distribution is governed by the 
familiar relation, 

‘ V 

p py CXP ( I T ). 
where p is the space charge per unit volume. 
Since the rigorous potential distribution for the 
case of a retarding field has never been investi- 
gated, it was believed desirable that this be done, 
in order to check the accuracy and range of Von 
Laue’s solution. 

The separation, a, between cathode and anode 
(fig. 3) is considered to be sufficiently small with 
respect to the area of the electrodes so that edge 
effects may be potential 
distribution is governed by 


neglected, and the 
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Py = 17 
dr ry stil 


where V is the potential, p the charge density, and 
¢=1/367.10~° farad per meter. 


The anode is assumed to be at a voltage V, 
sufficiently negative with respect to the cathode 


CATHODE ANODE 














Figure 3. Voltage distribution cf diede with retarding field, 


so that the slope of the potential V(sz) is always 
negative. 

The following analysis follows Fry’s analysis 
of the a-space precisely, and is outlined here for 
completeness. The only difference, of course, is 
that here the constant of integration is not zero, 
as it is when one has a potential minimum. If 
ry represents the initial velocity of an electron, 
and ¢ is its velocity at position v, 


| , » I 
=» Me? + eV=SMri; IS) 
where ¢ 1.60 10°" coulomb, is the charge of 


the electron, and 1/=9.10~"', is the mass of the 
electron in kilograms. 

At a given position v, those electrons having 
initial V2Mi<tiV  nevei 
reach s; those having initial velocities such that 
eV< 1 /2Mi<eV, twice, going and 


velocities such that 
pass by tr 
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Those having initial velocitis such 
-eV, are captured by the plate and 


Accordingly, the charge 


returning. 
that 1/2.\Jes 
pass the point only once. 
density is given by 


n(0) "yr n(v) 


diy + 2¢ 
v =o 
J yrVia 


where n(?) is the number of electrons emitted 
per second per square meter with initial velocities 
M. Substituting 


p(r) =¢ diy, (19) 


J yr 


between ¢p» and v%+deo, and r= 2 
eq 19 and 18 in eq 17, 


. v -~ f- l (Uo) dro 42 ye 
wear OS 
Integrating eq 20 between the limits 2 and a, we get 
(<)- rind ) 
dr /, dr 
2M ce : 
[ | n(t) yeR—rvV, doy 
€ wo yrt 


rv dt. 


n(V) dv 


= > 1 (20) 
ye—rl 


N() Vl 


yr hi . 
n(%) ye—rV dry + (21) 


« yrl 


Since the initial velocities are distributed 


according to the Maxwell-Boltzmann distribution, 


- exp ( _ =e) 


in which N is the total number of electrons emitted 
per square meter per second, & is Boltzmann's 
constant, k=1.38 10°", and T is the absolute 
temperature. 


a Ne . (23) 
Where /) is the saturation current density. If one 
changes to reduced variables, 


and 


where 
y2rM el, 
€ (kT)32 
and /, is the plate current, then eq 21 reduces, 
when n(v) is replaced by its value, (eq 22), to the 
following: 
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815197—49 6 


4] (3 a ” S 9 V7). 
(‘r) (jn). +6 (1+ erf yn) (14243) 


Here erf 9 is the error function, 


9 


erf x 


| exp(—f)dt, 
0 


V We 


and (dn/dt)2 is the value of (dy/d&)* at ¢& . that 
is, at the anode, z=a. 
For brevity, let 


fin) = (ert vnyer—(1 + > va): 


\r 
and 


(i) i 


26 may be written as 


dn\ sia 
(ut m +-f(m). 


Accordingly, eq 


Since ¢ is always negative, 


_—_— l" dn a 
J0 vf(n) -m- 


The value of m is determined by 
condition at the cathode. 


the boundary 


At the cathode, 1=0,é ayi=hn=eVUkT=m, 
so that 

£=- \" dn 

; Jo yf(n) +m 


fixes the value of m?. 


Since for large values of n, {(n)—-2e", the upper 
limit of the integral in eq 32 may be replaced by 


infinity, so that one has, approximately, 


a dy 

~ |, (f(n) + m?)b 

Since &; ts proportional to 1}, it may be inferred 
from the above equation that a large value of 
m® implies a small value of plate current, and that 
a small value of m? implies a relatively large value 
of plate current. 

North’s paper [10] tabulates /(n) (his ¢e()), and 
tabulates the solution of eq 14, for the case 
m=, (his &). Using the computed values of 
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t(n) from North’s paper, & was determined as a 


function of » by numerical integration, for selected 
values of m?, and the results are tabulated in 
table 1 for OSy»=20, and are graphed in figure 4 
for O=nX10. To determine €— for »>9, one 
observes that for the tabulated values of m’, 
f{(n)+m?~2e, for 7 >9. Accordingly for »>9, 
eq 14 yields 


E(n) =£(9) + y2(e-**—e ™), (33) 


from which &() may be calculated for larger 
values of » than are tabulated in table 1. 

To determine the potential distribution of a 
given tube, one first determines & and » from the 
given values of cathode temperature, plate cur- 
rent, or saturation current, plate voltage, and 
plate-cathode separation, a. Then, from the 


Values of & [° j dz 
vy (x) + 


TABLE 1. 





” m=0-— m=0.05—-F m=0 — m=0.2-E m 
“ o “ “ i 
ool 0. 0426 0. 0809 0. 0221 0 
a2 osi4 oeol (437 
Os 171 Ne) 652 
wt 1508 1147 OSSS 
Os 0. 4281 Isl4 . 1402 157 
we 2108 1648 1252 
07 2586 ISS4 . 1443 
(is 2650 2011 1628 
wn vo ed 2337 TSOs 
ih sl ‘M4 2542 . 19M4 
A S170 a) re | a SS 
w UTS wT 57s ASs7 
i 1. 10st “VIS HUT A006 
" 1. 2173 WT sol Us 
ro Lolo LL sy Ol Pei 
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Ww 2. 5444 2. 2402 2. 1352 2. O82 1 
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a 2. as 2. 2406 2. 1446 2.0176 1 





family of £, » curves plotted in figure 4, one deter- 
mines the particular curve of the family that 
passes through the point &, m. From this curve 
ne). (n< mé<§&), the potential Vo (2) within 
the diode may be determined from 


Viz) "1-7 
Vi mn 
2 &-—§ 
a &, 


according to the definitions of — and » (eq 24 and 


(34) 


; (35) 


2. Approximate Analytic Expressions for the 
Potential 


In order to obtain an analytic approximation 
for the potential, one must examine approxima- 
tions to f (m) that will make the equation 


as a function of n, for selected values of m? 
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Figure 4. 


(nm dx 
-t=| 


o vf(2) +m?’ 


integrable in terms of tabulated functions. For 
n large, f(n) ~2e", and as 7-0, f(y). 
For T=1,000° K, LT 11.7, and »=11.7(V,—V). 


For negative plate voltages of one volt or more, 
y will be larger than, say 2, over most of the volt- 
age range, so that the approximation /(n)~2 
(e’—1) will be fairly good over most of the range 
This is illustrated in table 2, where f(n) 
tabulated for 0<n<9. For 


of ». 
and 2(e’—1) are 


TABLE 2. A comparison of f(n) and 2(e"—1) 


n fin Je 1 
0 0 0 
| 0.2 0. 204 O44 
4 716 as 
6 1.27 1. 64 
& 1. us 2.46 
1 2. 88 3.44 
2 11.8 12.78 
; 36.9 38.18 
4 1. 6 107 
5 208 205 
6 sus sod 
7 2, 190 2, 191 
S 5, 060 5, 960 
oo) 16, 200 16, 204 
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large values of 9 the approximation is quite good, 
but becomes poor for 0<y<2. On the other 
hand, the error in — due to the approximation is 
the addition of m’, especially for 
Accordingly, let 


reduced by 
larger values of m’. 


dd Lv de 
=2 | 4 | a » (36) 
e Oye T 


0 y 2e7-+ m*—2 V2. 


where 


Equation 36 is a solution of 


dy 
a4 e", (38) 
de 


which is Von Laue’s expression [9], inferred from 
the assumption that equilibrium exists. It is 
equivalent to 
‘D) 
a - tee 
d*\ V kT ir 
> — € ’ 

dz’ € 
which in turn implies that the space charge 
distribution is 


eV 
P= Pot vail (40 
where 
2eM 
a=, V [T (41) 


If one integrates the rigorous expression for p in 
eq 19, one gets 


wv (14008 «[%=")) 
/2eM -r( terf V kT 


p lo 4 iT € 5) ° (42) 
Since 
orf Ver ~V)=erf Je~t, 


for 7 >2, it is seen that the approximate expres- 
sion for p approaches the rigorous expression 
fairly closely for »>2, and accordingly, the 
approximate expression for the potential must 
also approach the true potential in this region. 
For —1SA°<0, (OSm’?<2), the integral of 


eq 36 is given by 
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2 | om? 
—£ ¥ [cos (yi- 9 € ve) 
/ m - 
vio 2 


/ 2 
cos~ \! -"9 | (43) 


©. (m*=32), 


For A“ 
em), (44) 


g—§,= ¥ 2Ae-**- 


For A?>0, (m?>2), 


2 : im? 
* “ 1 haaiess 
~ [ inh Vv » 1 
sinh (e-v v3) (45) 


The value of m? in the above equations is fixed 
by the boundary conditions at the cathode, 
namely, n= When &=&,. For large values of 


m, eq 45 gives 


wer 
- 


m —} (46) 
and the resulting potential distribution is, of 
course, the same linear funetion of distance one 
gets in the absence of space charge. 

In order to compare eq 43, 44, and 45, which are 
Von Laue’s solutions, with the rigorous solutions 
of eq 32 obtained numerically, values of the two 
solutions are tabulated in table 3, for 7,10, and 
for m?=0, and m?=3. For m?=3, the agreement 
between the two solutions is excellent, and even 
for m?=0, the worst possible case for a given value 
of m, the agreement is fair. 

Unfortunately, the expressions for the potential 
obtained from eq 43 and 45 do not permit the time 
of transit of an electron as a function of its position 
to be calculated explicitly, which is necessary for 
an evaluation of the noise power spectrum. 

However, the potential distribution given by eq 
44 approximates fairly well the shape of the po- 
tential curve for small m?, that is, 0<m*< 20 
(large /,). In figure 5 graphs are plotted of the 
true reduced potential distribution, V/V, 
(m—n) m as a function of the reduced distance 
r/a=(&,—&)/&, for m?=0, m*= 10, and m?=20, for 
m —10. Also, the relative potential distribution 
obtained from eq 44 is plotted and is seen to fall 
about midway between the curves m?=0 and 
m’=—20. Accordingly, for the case of large plate 
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Figure 5. Comparison between approximate and rigorous 
potential distributions. 


A, m?=0; B, approximate formula; C, m?=10; D, m?=20; E, m=—m/g:. 


currents, one may use eq 44 for obtaining an ap- 
proximate potential distribution and thence cal- 
culating the noise spectrum. 


TABLE 3. 
VeVi =1l—am, as a function of the reduced distance, 


A comparison of the reduced potential, 


xa=1—£&, obtained from the rigorous solution (eq 14), 
and from Von Laue’s appreximation (eq 43, 45) 


mi=0 m=3 
. . Approxi- , Approxi- 
Rigorous mate Rigorous nate 
l—v 1—é/§, 1—é/é 1—t/é, 1—t/&, 
i- 
i) 0 0 0 i) 
| ol 0. 0025 0.0015 0.0045 0. 0048 
2 LLied 0081 O12 O25 
3 0130 0138 0247 0264 
| 4 O24 O24 O70 . 0466 
5 042 470 O7S89 O84 
| .6 0719 0814 . 135 . 139 
we . 1221 1385 227 . 23 
s a ss | 379 386 
4 371 412 fay . 0 
95 522 wie 79 ToS 
us ory ra) 1 v6 915 
Lo 1000 1.000 1.000 1.000 





For the case of small plate currents, and large 
negative anode voltages, the reduced potential 
distribution approaches a straight line of unit 
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slope, also graphed in figure 5 under the label 
m=—m/i. Accordingly, the following section 
proceeds with the calculation of the noise-power 
spectrum for the two cases, large plate current 
(m? small) and small plate current (m? large). 


IV. Noise Power Spectrum for a Diode 


1. Case 1: Noise Power Spectrum for Large Plate 
Current (m? small) 


If one substitutes the values of — and yn from eq 34 
and 35 into eq 44, one gets for the potential 
energy of an electron, U, 


U=Ve=2kT In(1+6n), (47) 


where 
— at ——s (48) 


In this paper only the contribution to the noise 
spectrum made by those electrons that never 
reach the plate, but which return to the cathode 
will be considered, since the effect of those elee- 
trons which constitute the plate current has been 
treated by Rack [11]. Accordingly, one considers 
only electrons whose initial kinetic energy, £, is 
less than eV}. 

For an electron emitted at t=0, 
which it has position sz is given by 


tad dx 
{= —_— — . ‘4é 
J, ay 2/M \ E-l (x) (49) 


the time at 


If one substitutes for ("(x) its value given by eq 47 
and integrates 


de /xM exp(E/2kT - 
‘WVkT 2% erty/ oT 


ety ste In( l — bx) | (50) 


for OStSt)/2, where fy is the transit time of the 
electron. 


, —jut, 
21 exp i $ — 
&(w) — ( sin 
yad 0 
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wl erf J _- In (1+ 62) 


Sat eu 7(7) «4. [£ (51) 
OV ET erty op 7” 


so that 
te erty) ofp —In(1 + br) 
erty or 
for OStSt,/2. 
According to eq 14 and 16, the noise power, 


Px, absorbed by each load in the frequency range 
Af, due to electrons with initial energy F, is 


> ely aw paw ne 
Ps: 2 tie on (53) 
where 
ten. "YY — ’ 
ata fr" texp( jot)dt (54) 
0 yad 


and where /, is that portion of the saturation 
current composed of electrons of energy F. 

Since t(z) and ¢t’(x), the first and second times 
at which the electron passes the point ., respec- 
tively, are related by 


t’ - é—t, 


and since 
—r(t), (55) 


v(t—t) 
w=} it) dt 
8(w) Vea , v(t)exp(— jot) 


1,/2 
{° r(tyexp[jw(t— i]t. (56) 
0 
Letting 


¢,/2 
8; (w) i.” r(t)exp(—jot)dt | exp(—jwt)dr 
« 0 Jf 
(57) 
where z,, is the maximum value of s, 


8(w) Jed [s:(w) —exp (— jwty)s;(—w)]. (58) 


Substituting eq 52 in eq 57 and 58, 


dr. (59) 


, BE 
2 erf vor 








To simplify the above somewhat, the substitution 








ya) ayn (1+ 62), 


is made, and one obtains 


, Ek t 
41 exp (5p 7- je 2") ev f2 
sald bad ‘’ 
Let 
er 
A( (24 er) Tig ye~” sin | — Ae d 
{2 2 erf y RT 
(62) 
Then 


— E—eV,\ go bk sae 
8(w)|"=—9 exp ( kT ') (arr) (63) 


where the value for 6 given by eq 48 has been sub- 
stituted. 

Since the energy distribution of the electrons is 
Maxwellian, 


Te a exp (—E/kT)dE, 


(64) 
and since the resultant noise power due to un- 
correlated noise sources is simply the sum of the 
powers of the individual sources [5], the resultant 
noise power spectrum due to all the electrons that 
return to the cathode is 


ly 
A (s/t) air) ES (65) 


where P represents the noise power absorbed by 
matched load in the interval Af=Aw/2r. 


Selpa /y 


™ ray |, 


each 


* VE/2T 


A(, £/2kT)= | ue~ sin Ye exp (“5 oy 4. ') orf u | du. 


0 


Integrating the above equation by parts, 


A(y k/2k7 


ir 9 
ye * sin 


(60) 


wt, erf y | 


dy. (61) 
lo erf 


van 


The diode shunted by two matched transmission 
lines is now replaced by its equivalent circuit, 
namely, a diode shunted by impedance 7,/2. The 
equivalent mean square fluctuation current, 7? is 
defined by 

2P=7,/2%, (66) 
where Zo= yu/ea/d is the characteristic imped- 
ance of the line. If the impedance of the diode is 
large compared to 7, which is assumed to be true, 
? represents the equivalent mean square fluctua- 
tion current that would flow through a diode short- 
circuited at frequency f. 


From eq 65 and 66, 


parr Vidas, (67) 

where 
F(V,,f) i “AX \VEPRT)E, (68) 

Since 


_ who E—e\ 
Qerf vE kT" rM /4kT exp ( kT 
from eq 48 and 51, letting 


c=dw yrM/4kT, (69) 


(70) 


sc [ (vEur E-eV ) 
a (Eel ‘ [, cos , ¢ exp aL? ‘| erf u ' du— 


vE/2kT cos a Fn ert YE2KT! | 








One observes that A( y£/2k7)=0, for c=0 and 
c=, Since the expression within brackets is 


always less than 2, £/2kT in absolute value, 


Vary BE /2kT 
i — ¢ V; 


k 
CeXP] opT 


ie ) my a , ak T 
, AWvERET dE< 2c2exp(—eV{/kT)’ 


(73) 


A(, £/2kT)< 


and 


F(V;,f) 


e 


which gives an upper limit to F(V,,f) for very high 

frequencies. 
€ 

A(VE/2kT) = 

Since A®(yE/2k7T) approaches zero as kT 

for kT small, and increases exponentially as 

exp(E kT) for E/kT>1, and since one is interested 

My Vi 
in A°*dE for an upper limit F/k T=eVi/kT > 10, 


JV 


the contribution to the integral from large values 
of the argument far overshadows the contribution 
from small values of the argument. Hence, we 
may replace the value of A(y£/2kT) by its 
asvpmtotic value for large values of the argument, 


c . wai ( es f ) 


\ 


A(yE/2kT) =~ (77) 
Physically, this means that despite their pre- 
ponderance in numbers, the low-energy electrons, 
because of their small excursions toward the plate 
times, contribute 
the low-frequency 


transit 
toward 


and consequent short 


practically nothing 
noise spectrum. 
Accordingly, 


"ely o T'n2 
; dE a ’ 


. 


7 "eV; 
IyeZ, - 


P=3IT 


(4/eV,)? 


where 


. = = Vi 
c exp/( ah: P)Fespieraen 
9 


E, , : *\2 
“er feos aoyE sinoykl 


Y 


Ve 


(a) Low-Frequency Approximation 


For a cathode temperature of about 1,000°K, 
and with a, the cathode-anode spacing expressed 
in millimeters as a,,, one has 


c< 5X 10 fay. (74) 
Since ESeVj, one may approximate the sine of the 
argument in eq 71 by the argument for frequen- 
cies such that ¢<0.3. 
rE 
ToL S—eV,] (VT ; 
A(,yF/2kT)=c¢ exp _ car | ue“ erf udu. 
2k 7 J0 
(75) 


Or, integrating by parts, 


erf y E/2kT—erf VE /2k | 


and 
Pel pa*w? a4 Af. 
(The above formula holds only for 
C=dwyrM/AkT<0.3). 
2. Case 2: m* Large (J, Small) 
For this case, the potential energy of the elec- 
tron is 


eV iz 
a 


U(r) 


and the time, f, an clectron reaches a point sr, is 


if \ 1h" | O=x1<1,/2. 


one 


given by 


Following the same procedure as in case 1, 
finds for the noise power absorbed by each load 
in the frequency interval Af, 


oyk ek ) dEAY, 


0 wy2M aleV,. 


The equivalent mean square fluctuation current through the diode is 


2] 


2 er vy | 
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E —— - 
(ier feos oyF sin oy EV 


—_ dEAf. 
o \ BE ok 5 I Af} 








For low enough frequencies, such that 
oyeV,<0.3, the expression in brackets approxi- 


mately equals ek 9. Since also, 


«Vy -2 
( FerdE=b(kT)‘, for eV\/kT >10, (81) 





en *¢ long 


(Ser 


For large values of Vv, or for high frequencies, 
the upper limit in eq 80 may be taken as infinity 
with negligible error. Accordingly, eq 80 may be 
rewritten as 


(82) 





telat 


In appendix 2 it is shown that 


o\k ek ) 


os oy _ 


cos oyF sin oyEU dE- Cet a4 owe «F(5, 


sin o \ EV 





where 
(85) 


and F is Kummer’s function |12], also known as 
the confluent hypergeometric function. 
Accordingly, the expression P 


for ¢ be 
written, 


may 


ee ( 0 ’ ‘ 
=. ~ +2w)e-"F (5 w)— 1 haf. (S86) 


(i: ny 4 


V. Equivalent Temperature of the Con- 
ductance Due to Space-Charge 


The shunt conductance, g, of the space-charge 
within a diode has recently been calculated by 
Begovich [13], who finds® 

a 


(a) 


. . z FR : 
where w is defined above (eq 85) and F (5:5: w) is 


ey : 2w)e-" 5H —1 . (87) 


the confluent hypergeometric function [12]. 


If the above value of g, and the value of 7? 
given by eq 86 are substituted in the Nyquist 
formula [24], 


P= 4kT gd f, 


where 7, is the equivalent temperature of the 
conductance, it is immediately found that 


* The writer is indebted to N. Begovich for a personal communication 
amplifying the material contained in reference [13]. 


dEAf. (83 

oyE e@E ) ~- 
; w)-1}, (84) 

Zee 3; (88) 


in agreement with the experimental results ob- 
tained by Van der Ziel [3]. 


VI. Summary 


1. A general expression for the noise power 
spectrum generated by the random emission of 
electrons of arbitrary trajectories within a wave- 
guide is developed in section I, 

A numerical solution of the potential dis- 
tribution within a plane diode is effected in see- 
tion II, and a family of curves from which the 
potential distribution may be obtained is graphed 
in figure 4. <A graphical comparison between 
several actual potential distributions and_ their 
approximate representation is given in figure 5. 

3. The results of sections I and II are utilized 
to derive the equivalent mean square fluctuation 
current due to the space charge within a diode 
for two types of a retarding potential distribu- 
tion, namely, the linear distribution, and those 
distributions whose slope at the plate is small. 

4. Utilizing the results of section III, and a 
recently computed value of the conductance due 
to the space charge within a diode. the equivalent 
noise temperature is shown to be equal to the 
cathode temperature. 


The author is greatly indebted to K. F. Herzfeld 
for his many helpful discussions and generous 
guidance during the course of this work. Also, 
the writer is extremely grateful to H. Lyons, chief 
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of the Microwave Standards Section, National 
Bureau of Standards, for his wholehearted coopera- 
tion and valuable criticisms. 
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VIll. Appendix 1 


Let x,y be generalized coordinates suitable for 
describing the guide surface S, and z the usual 
cartesian coordinate which is parallel to the axis 
of the guide. The normal modes, 6, of electric 
type for a cylindrical guide are a doubly infinite 
set obtained by two independent processes [15,4]. 


6,=VU:x<k, (89) 
where k is a unit vector in the z-direction, and 
where 

VU dry) + BU (ry) =0, (90) 


and n-7U';=0 on S, the surface of the guide. n 
is normal to S. 


b= BrlrwkFjaVe(r,y), (91)? 
’ The rj should not be confused with the letter ¢ used in section I to represent 
elocity. 
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where Y*v,+ 8?v,=—0, and v,=0 on S, except for the 


TEM ease for which 8,=0, and 


nx Ve.=—0o0n S. (92) 


(The 8, are, in general, different for the (7°, and 
v;.) 


j 


If the WU’, and v, are normalized, that is, 
Selda= ferda=1, (93) 


where the integral is taken over the transverse 
plane S’ of the guide, then the following orthogonal 
properties may be demonstrated [15,4]: 


Seb;- b,da=5),, (94) 


where 6,,=0, for 7#h, and 6,,—1, for j=h. 


Je (6; curl b,)-kda=0 for all j and h. (95) 


The phase constant a corresponding to each b, 
is given by 
a’ = pew” — B?, (96) 


and L,=6, exp (+ jaz) satisfies the vector wave 
equation 
curl curl Lj= yew*L,. (97) 


Utilizing the above properties of the normal 
modes, the field within a guide excited by an 
arbitrary current distribution will be determined 
by a method due to Bethe [4]. 

Using Gauss’ Theorem, if A and D are two 
arbitrary vectors, 


[14 <curl D— Dx curl A] - nda= 
| (D- 9X VXA-A- TX VX Dd, (98) 


where n is the outward drawn normal to S, en- 
closing the volume V whose element is dr. If one 
lets A= E, the electric field which satisfies 


curl curl E— pew*E= —joyd, (99) 
where J is the current density, and lets D be a 
solution of 
VX<V7X< D—pvew’*D=—0, 


then eq 98 becomes 


| 2 VX D—DX 7 XB) -nda=~yjw | J+ Dar. 
(100) 
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Consider a current element within the guide at 
z=2’, and take matched loads at the two planes, 
z=0 and z=2,. Expanding the field - in terms 
of the normal modes, let 


E=S0A,b, exp (—ja;z)z>2’ ) 
: (101) 
E=SOC,/6, exp (+ja;z)z<2’ ) 
Since E is continuous at z=’, and since the b, 


are orthogonal, 


A, exp (—Jjaj2’)=C, exp (jaj2’). (102) 

If one substitutes the value of E given by eq 101 
in eq 100, and lets D= b, exp(—ja,z), and inte- 
grates the left side of eq 100 over the plane z=0, 


one gets 


= I p>» C,[b,> < (b, exp(—ja;z)) 
—b, <7 x (6, exp + ja,2z)|-kda, 
or, 
—SOS [Or K7xb —ja,b, K(k <b, 
h 
—b ‘ V . b,—jayb; . (k . b,)|-kda. 


Because of the orthogonality relation (eq 96), 
the above simplifies to 


“ie C, [—ja; fb, -b,da—ja, fb,-b,da=2Cja,j. 
The surface integral over the plane z= 2, vanishes, 


and one has, from eq 100, 


C, Qa, ef J-b,exp(—ja;2’)dr. (103) 
Substituting eq 103 and eq 102 in eq 101, one 
finally gets 

E — Hsp, POs ——a =i) g 
- j a 
(104) 
IX. Appendix 2 
Let 
] | kx ‘T (cos oy E_ sin oy EV GE, (104) 

Jo l oyE ek )S 

feos r_ sin zr)? ows , (2x) *-? (7—1) . 
— 2 -— . ; ’ 0: 
a 2 | 1) (ay! Q+1) ny 
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which may be verified by performing the square, 
expressing the result in trigonometric functions of 
2x, and then substituting for each trigonometric 
function its corresponding infinite and 
combining. 

If we substitute, eq 105 in eq 104, we obtain 


series 


(j—1)¥(7+2) 
(27)!(9+-1) 
(106) 


T= (kT) > (—1)/2—""(eXk TT)" 
j=l 


where we have utilized the formula for the gamma 
function, 


, (ii) 
| zr" l¢ “ay ¥ ‘ 
e 0 


a" (107) 
Using the duplication formula [16] for the y- 
function, 
: ‘ . 2°) : i. a 
(27)!=2py(29) OA bei Nx(j+5) 
\r « 


eq 106 simplifies to 


(kT 4 (3\ 2 (—@kT) | L\S (—okT)! 
=. (14D) a) 
2 IH /.,3 
(i+3) W(i+3) 
(108) 
From MacRobert’s definition of the confluent 


hypergeometric function, it follows that 


¥(p) w"y(a+n) 


(109 
y(a) nao ¥(N-+ 1) ¥(p ) 


F(a,p,w) 


Accordingly, eq A 109 may be written 


OD Aide) +a(1—) 


where w=o0°7kT 


[= (110) 


However, since 


F(1,o—1,—w)=1 — wF(1,p,—w), (111) 
I sr a4 20) F(1,5,—w)—1)- (112) 
2 { 2 j 
Finally since 
F (a,p,w) =e" F (p—a,p,—w), (113) 
ie 
I ee) [1+ 2w)e"F( 5-5) (114) 


Wasuineton, October 6, 1948. 
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